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PREFACE 

The  work  reported  herein  was  conducted  by  the  Arnold  Engineering 
Development  Center  (AEDC),  Air  Force  Systems  Command  (AFSC).  The  results 
presented  were  obtained  by  ARO,  Inc.,  AEDC  Division  (a  Sverdrup  Corpora- 
tion Company),  operating  contractor  for  the  AEDC,  AFSC,  Arnold  Air  Force 
Station,  Tennessee.  Elton  R.  Thompson  was  the  Air  Force  project  manager. 
The  work  was  done  under  ARO  Project  No.  E32A-P0A,  and  the  manuscript  was 
submitted  for  publication  on  October  5,  1978. 

The  author  of  this  report  is  a Computer  Applications  Analyst  in  the 
Computational  Services  Branch,  Facility  Support  Department,  ARO,  Inc. 
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1.0  INTRODUCTION 

Theoretical  calculations  of  exhaust  nozzle  performance  are  often 
required  to  aid  in  the  evaluation  of  engines  tested  at  the  Arnold 
Engineering  Development  Center  (AEDC).  In  the  past,  most  propulsion 
nozzles  have  been  axisyimuetric  and  well-established  computer  programs 
are  available  for  computing  the  flow  field  in  such  nozzles.  However, 
future  tests  at  AEDC  will  involve  engines  with  three-dimensional  (3-D) 
exhaust  nozzles,  and  no  programs  have  been  available  for  computing  the 
flow  field  in  3-D  nozzles.  Consequently,  development  of  3-D  computer 
programs  for  both  the  transonic  and  supersonic  portions  of  the  flow  was 
initiated.  This  report  describes  the  supersonic  computer  program,  which 
is  based  on  the  method  of  characteristics  (MOC).---Of  course,  the  compu- 
ter program  described  herein  is  not  limited  to  computation  of  propulsion 
nozzle  performance,  but  is  also  applicable  to  many  other  3-D  supersonic 
internal  flows. 

Except  for  the  bound ary- layer  region  near  the  wall,  supersonic 
nozzle  flow  can  be  computed  with  adequate  accuracy  by  assuming  the  fluid 
to  be  inviscid  and  adiabatic.  In  addition,  most  nozzle  flows  do  not 
contain  strong  shock  waves,  so  the  nozzle  performance  can  be  adequately 
predicted  by  assuming  the  flow  to  be  shock-free.  However,  because  of 
flow  phenomena  upstream  of  the  nozzle  entrance,  the  flow  in  propulsion 
nozzles  is  often  significantly  rotational.  The  rotatlonallty  of  the 
flow  entering  the  supersonic  region  persists  throughout  the  flow  field. 

In  this  study,  the  rotational  MOC  was  chosen  as  the  basis  of  the 
numerical  analysis.  The  mathematical  theory  of  characteristics  is 
well  established,  and  many  computer  programs  have  been  based  on  the 
method,  particularly  for  planar  and  axisyrametric  flow.  Several  programs 
have  been  developed  to  compute  the  flow  over  3-D  bodies  (e.g.,  Refs.  1 
and  2)  and  special  types  of  3-D  internal  flow  have  been  solved  with  the 
method  (Ref.  3).  However,  few  attempts  have  been  made  to  compute 
general  3-D  supersonic  internal  flows  with  the  MOC. 


5 


2.0  GOVERNING  EQUATIONS 


In  this  section,  the  governing  flow  equations  and  the  resulting 
characteristic  equations  are  presented.  These  equations  are  well  docu- 
mented in  the  literature  (e.g.,  see  Ref.  1);  therefore,  no  development 
will  given. 


2.1  EQUATIONS  OF  MOTION 


The  steady,  inviscid,  3-D  flow  equations  for  an  ideal  gas  are: 


Continuity: 


Momentum: 


Energy : 


+ 

+ w k. 
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2.2  CHARACTERISTIC  EQUATIONS 

The  equations  of  motion  give  rise  to  two  sets  of  characteristic 
surfaces.  The  defining  equations  are 

2 

(ufx  + vf  + wfz)  =0 

and 

(ugx  + '•'gy  + ws/  - a.2  (sx  - Sy  + g2)  = 0 CB) 

where  f (x,  y,  z)  = 0 and  g (x,  y,  z)  = 0 are  the  characteristic  surfaces. 
The  first  surface  is  composed  of  streamlines  and  the  second  is  the  Mach 
conoid.  The  equation  of  a ray,  or  bicharacteristic,  of  the  Mach  conoid 
may  be  expressed  as 


dx  = (cos  /3  sin  9 ~ sin  jS  cos  8 cos  5)  dL 


(9) 


dy  = ^cos  |8  cos  8 sin  ip  — sin  (sin  8 sin  ip  cos  8 - cos  \!i  sin  sjjdL 


(10) 


dz  = os  ,S  cos  8 cos  i (r  - sin  /3  (sin  8 cos  dj  cos  8 +■  sin  d/  sin  dL  (11) 

where  3 is  the  Mach  angle,  dL  is  the  distance  along  the  bicharacteristic, 
(5  is  a parametric  angle,  and  0 and  i|/  are  related  to  the  velocity  vector 
by 


u = q sin  9 


(12) 


v = q cos  8 sin  \p 


03) 
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w = q cos  0 cos 


(14) 


The  parametric  angle  6 lies  in  a plane  normal  to  the  velocity 
vector  and  is  measured  from  the  plane  containing  q and  x.  The  rela- 
tionships between  the  variables  are  shown  in  Fig.  1. 


X 


Figure  1.  Coordinate  system. 


2.3  COMPATIBILITY  EQUATIONS 


The  compatibility  equations  are  determined  by  the  flow  equations 
and  the  requirement  that  the  derivatives  in  a direction  normal  to  the 
characteristic  surface  disappear.  The  compatibility  equation  which 
applies  along  the  Mach  conoid  is 
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°ot  ft  gE  + cos  S — + cos  e Sin  5 itt 

pq2  dL  dl.  dL  (15) 

- sin  p ( cos  & cos  5 OiL  - sin  S = 0 
\ <9N  dN/ 

3 3 

where  — and  are  derivatives  along  and  normal  to  the  bicharacteristic. 


The  compatibility  equations  along  a streamline  are 


RtlT  = 1 di. 

y- 1 p 


<16) 


and 


1 dp  =■  - q <lq  (17) 

3.0  NUMERICAL  PROCEDURE 

The  numerical  procedure  is  similar  to  the  one  developed  by  Strom 
(Ref.  1)  for  external  flow.  The  method  traces  streamlines  from  a known 
reference  plane  to  the  next  reference  plane  as  illustrated  in  Fig.  2 for 
a field  point.  The  flow  properties  at  approximately  equally  spaced 
points  on  an  initial  plane  normal  to  the.z  axis  are  assumed  known.  The 
reference  planes  throughout  the  flow  field  are  assumed  to  be  normal  to 
the  z axis. 

3.1  METHOD  OF  SOLUTION 

A new  reference  plane  is  located  a distance  dz  from  the  initial 
plane.  This  distance  must  be  determined  such  that  the  Courant-Friedrichs- 
Lewy  (C-F-L)  stability  conditions  are  satisfied  (Refs.  1 and  4).  The 
C-F-L  stability  conditions  are  satisfied  if  dz  is  smaller  than  the 
minimum  intersection  distance  of  the  Mach  conoids  from  the  initial 
points. 
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• Known  Points  in  Initial  Plane 
o Point  in  New  Reference  Plane  (P2) 
x Points  in  Initial  Plane  on  Bicharacteristics  from  P2 
— Streamline 
- — Bicharacteristic 

Domain  of  Dependence  of  Difference  Network 

Domain  of  Dependence  of  Differential  Equations 

For  the  C-F-L  stability  conditions  to  be  satisfied,  the  domain  of 
dependence  of  the  differential  equations  must  be  contained  within 
the  domain  of  dependence  of  the  difference  network. 

Figure  2.  Field  point  network. 


The  intersection  of  a streamline  from  a point  in  the  initial 
plane  with  the  new  reference  plane  locates  point  P ^ . Four  bicharac- 
teristics, 90  deg  apart,  from  P 2 to  the  Initial,  plane  yield  points  P.^, 

P4»  P5*  and  V 

The  properties  at  and  eight  surrounding  neighbors  in  the  initial 
plane  are  surface  fitted.  The  flow  conditions  at  points  Pj  - Pfi  are 
evaluated  from  the  surface  fit.  The  compatibility  equations  are  applied 
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along  the  bicharacteristics  and  the  streamline  to  obtain  the  flow  condi- 
tions at  P^'  This  procedure  is  iterated  until  the  flow  conditions  at 
converge.  When  all  points  on  the  new  plane  are  determined,  the  new 
plane  becomes  the  initial  plane  for  the  next  calculation. 

3.2  WORKING  EQUATIONS 

Several  of  the  computational  steps  are  the  same  for  both  field 
points  and  body  points.  These  common  phases  of  the  procedure  will  be 
discussed  first.  Then  the  computational  steps  required  for  field  and 
body  points  will  be  presented  individually. 

3.2.1  Common  Procedures 

The  surface  fit  of  the  flow  properties  at  and  eight  neighboring 
points  in  the  initial  plane  is  made  once  in  the  iteration  procedure  to 
obtain  and  ^ts  properites.  The  fit  chosen  is  a spline  surface  fit 

(Ref.  5)  which  fits  a variable  w as  a function  of  x and  y: 

9 

W (x,y)  = aQ  + Uj  x + a2y  + ^bjrflnr?  (18) 

i^l 

where 

ri  = (x  ~xi)  + (Y  -Vi)2 

Five  variables  are  fit:  p,  p,  q,  9,  and  il>.  The  surface  fit  is  used  to 

g 

evaluate  the  variables  at  P,  - P,.  The  partial  derivaties  -r—  and 
g 3 6 9x 

are  also  evaluated  using  this  fit: 

^ = &l  + 2j^b.  (l + ]nrf)(K-xi)  09) 
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dw‘ 

— = n„  + 

dy  2 


2Sbi  6 + ln  rf)(y  - >;) 


(20) 


i=i 


The  derivatives  and  used  in  the  compatibility  equation  along 

a bicharacteristic  must  be  obtained.  The  procedure  for  evaluating  them 

3 0 

is  the  same;  therefore,  only  the  procedure  to  obtain  — will  be  given. 
.30 

Write  -rrr  as 
3N 


fdO\  _ (<W\  /<?x\  +/dd\  /<9y\  + (dd\  /Sz\ 

\dN/  " \<9x/.  \dW  . W/  i \dW  . wv  ; W*7( 

' l • 1 I 1 1 


(21) 


3030  30 

The  derivatives  t—  and  - — are  obtained  from  the  surface  fit.  Then,  — 

Sx  3y  3z 

is  obtained  by  writing  the  equation  as  a difference  equation  and  solving 

J " (|),  C22) 


The  derivatives  normal  to  the  bicharacteristic  are  obtained  from  the 
coordinate  transformation  as 


(— \ = -costfsinS; 

VdN/.  1 

(23) 

(dy\  m 
\<?N/ . 

I 

sin  Q sin  <A  sin  Sj  + cos  cos  Sj 

(24) 

= 

sin  6 cos  t !j  sin  5j  - sin  tfj  cos  Sj 

(25) 

3 0 

The  final  form  of  — becomes 

3N 


- (M\(ch\  * (jW\  fdy\ 

" \dy  ).  \dK). 

- [*■♦("),  (,«-x')t(f)1  (,2-y‘)] 


(26) 
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For  the  first  Iteration,  the  flow  properties  at  P ^ are  assumed  to 
be  the  same  as  at  . Thereafter,  the  last  computed  values  are  used. 

In  the  sections  on  the  field  point  and  body  point,  the  average  values  of 
flow  properties  are  the  averages  of  the  two  points  involved  in  the 
calculation  (for  example,  9 = (8^  + 

3.2.2  Field  Point  Routine 

Given  an  initial  field  point  and  eight  neighboring  points  in  the 
Initial  plane  (Fig.  2),  spline  surface  fits  of  p,  p,  q,  0,  and  ip  as 
functions  of  x and  y are  made.  The  iteration  for  the  point  and  its 
flow  properites  consists  of  the  following  steps. 

Step  1, 

The  intersection  of  the  streamline  from  P^  with 
plane  locates  the  new  field  point  P2. 

x2  = Xj  - sin  6 dz/(cos  0 cos  t/i) 

y 2 = + sin  0 dz/cos  $ (28) 

z2  = Z,  +■  dz  (29) 

Initially  6 and  ip  are  the  values  at  P^.  Successive  Iterations  use  the 
average  of  9 and  if)  at  P^  and 

Step  2 

Four  bicharacteristics  extended  from  P2  to  the  initial  plane  give 
the  base  points  P^,  P^,  P^,  and  P^,  The  flow  properties  at  are 


the  new  reference 


(.27) 
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assumed  to  be  the  same  as  those  at  for  the  first  iteration.  There- 
after, the  flow  properties  at  P2  are  the  last  values  computed.  The 

bicharacteristics  are  located  at  the  parametric  angles  5 » 0,  tt, 

• 3ir  L 

and  — The  equations  for  the  base  points  are 

x;  = x g (cos  /3sin  6 - sin  jS  cos  8 cos  fij)  dL,  (30) 

Y j = y%  -(cos  /3  cos  6 sin  if/  - sin  /9  (sin  ff  sin  if/  cos5;  - cos  i/r  sin  sSjdLj  (31  ) 


(32) 


dLj  = dz/|cos  /3  cos  Q cos  if/  - sin  /8  (sin  8 cos  if/  cos  5j  -r 


sin  if/ 


>0] 


(33) 


For  the  first  iteration,  the  angles  &,  0,  (fi,  and  6 are  those  at  P2< 

Successive  iterations  use  the  average  of  the  values  at  P„  and  P . . 

2 1 

Step  3 

The  compatibility  equations  along  the  bicharacteristics  are  solved 
for  p,  0,  and  ip  at  P2>  Only  three  of  the  bicharacteristics  are  required; 
however,  to  improve  accuracy,  four  solutions  are  obtained  using  three 
bicharacteristics  at  a time.  The  results  of  the  four  solutions  are 
averaged  to  obtain  the  values  of  p,  0,  and  iJj.  The  compatibility  equation 
in  difference  form  is 

Y (P2  ~ P>)  + cos  Si  @2  - + cos  6i  sin  5i  (^2~:^i) 

Pi  **  go  (34) 

+ sin0.  (cos  8-  cosSj^).  - sin  5^)1  dL.  =0 
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The  equation  contains  the  three  unknows  and  ^2'  T**e  compati- 

bility  equations  using  P^»  P^,  and  are  solved  simultaneously  for 

P21*  021*  and  ^21  “ simllarily»  <p4>  p5»  p6>»  (p5»  pg»  p3)»  and  (p6*  p3> 
P^)  are  used  to  obtain  values  of  p,  0,  and  ty.  Thus,  * (P21  + P22  + 

P23  + P24^^*  *n  tlie  compatibility  equation,  p^,  p^,  q^,  and  6^  are 
the  average  of  the  values  at  P^  and  P^. 

The  variables  p^,  T^,  and  qj  are  obtained  from  the  compatibility 
equations  along  the  streamline.  The  flow  along  a streamline  without 
shocks  is  isentropic.  Therefore,  the  equations  integrate  to  the 
following  forms  for  an  ideal  gas: 


y-l 


12  - V2  FT  <T.»  - T2> 

where  and  p are  the  stagnation  values  for  the  streamline. 

Step  4 

The  preceding  steps  are  repeated  until  the  flow  conditions  at 
converge.  If  there  is  no  convergence  after  25  iterations,  the  results 
are  checked  with  a reduced  convergence  criterion.  If  this  is  satisfied, 
the  point  is  accepted  and  a message  printed.  Otherwise,  the  program 
terminates. 


(36) 

(37) 
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3.2.3  Body  Point  Routine 

Given  the  properties  of  an  initial  body  point  and  those  of  eight 
neighboring  points  on  the  initial  plane  (Fig.  3),  spline  surface  fits 
are  made.  The  body  surface  is  assumed  to  be  given  by 

B (x,  y,  z)  = 0 (.38) 

where  B is  either  a known  function  or  a surface  fitting  element  in  the 
region  of  the  body  point.  The  iteration  for  the  new  body  point  P^t  in 
the  new  reference  plane,  and  its  flow  properties  consists  of  the 
following  steps: 

Step  1 

The  new  body  point  P£  is  located  at  the  intersection  of  a plane 
through  P^,  defined  by  the  body  unit  normal  and  the  unit  velocity  vector 
tangent  to  the  body  at  P^,  with  the  body  surface  at  the  new  reference 
plane.  This  requires  the  simultaneous  solution  of  the  following  equations 

B (x,  y,  z)  = 0 (39) 


and 

cos  0 cos  41  — cos  6 sin  ^x 
- 112  cos  6 cos 

where  n^,  n^,  and  n^  are  unit  normals  to  the  body  surface.  The  first 
iteration  uses  the  values  of  &,  i|j,  n^ , n^,  at  P^.  Successive 

iterations  use  the  averages  of  the  properties  at  and  P£. 


xj)  + (n  j sin  6 

y-j'J  = sin  6 - nj  cos  6 sin  dz 


(40) 
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Step  2 


This  step  is  similar  to  Step  2 of  the  field  point  routine.  However, 
only  three  bicharacteristics  are  used  for  the  body  point  calculations. 

The  parametric  angles  used  are 


S*  = 


sin  t)  cos 


* 


"2 


6 - 


(41) 


S 


4 


53- 


''2' 


(42) 


8 


'2 


(43) 


Step  3 

The  compatibility  equations  along  two  bicharacteristics  are  solved 
simultaneously  with  the  condition  for  flow  tangency  to  the  surface  at 
P2.  The  equation  for  flow  tangency,  which  provides  an  additional  rela- 
tionship between  0^  and  may  be  expressed  as 

n ^2  $2  cos  VJ2  + n 2 2 sin  ^2  ^ n32  tos  ^2  a*n  ^2  = 0 (44 

where  n^2»  n22*  an<*  n32  are  un*t  normals  to  fche  body  at  P^.  A Newton- 
Raphson  routine  is  used  to  solve  this  system  of  equations.  The  base 
points  and  P^  are  used  to  obtain  P2^,  021*  an<*  ^21^  P3  ant*  ^5  are 
used  to  obtain  p22>  ®22’  ant^  ^22’  ^bese  two  results  are  averaged  to 
obtain  the  values  at  p2  (p2  = (p21  + p22)/2). 

The  compatibility  equations  along  the  streamline  are  solved  the 
same  way  as  for  the  body  point  yielding  p2,  q2»  and  T2* 
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Step  4 


The  Iteration  procedure  for  the  body  point  is  the  same  as  the 
procedure  used  for  the  field  point. 


• Known  Points  In  Initial  Plane 
o Point  in  New  Reference  Plane  and  Body  Surface  (Pg) 
x Points  in  Initial  Plane  on  Bicharacteristics  from  P2 

Streamline 

Bicharacteristic 

Domain  of  Dependence  of  Difference  Network 

Domain  of  Dependence  of  Differential  Equations 


Figure  3.  Body  point  network. 


3,2.4  Computer  Program 


The  preceding  equations  were  programmed  in  FORTRAN  IV  language  for 
solution  on  an  IBM  370/165  computer.  A description  and  listing  of  the 
program  is  given  in  Appendix  A,  and  an  example  problem  is  given  in 
Appendix  B.  The  program  requires  a storage  capacity  of  approximately 
192,000  bytes. 
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4.0  RESULTS  AND  DISCUSSION 

The  present  3-D  MOC  program  was  evaluated  by  computing  the  flow  in 
a typical  axisymmetric  nozzle.  With  y ".1.24  and  uniform  flow  in  the 
throat,  the  nozzle  produces  an  exit  Mach  number  of  4.1.  The  3-D  MOC 
computations  were  done  in  three  different  ways  (see  4) . 

([and 

z Axis  - Case  1 

z Axis  - Case  2 


i Axis  - Case  3 

Figure  4.  Geometry  for  test  cases. 

Case  1 - The  z axis  was  aligned  with  the  geometric  centerline 
of  the  nozzle. 

Case  2 - The  z axis  was  parallel  to  the  nozzle...  centerline  but 

was  displaced  by  a distance  of  10  nozzle  throat  radii; 
thus,  the  z axis  for  the  flow  field  is  completely  outside 
the  nozzle. 


Final 


19 


AEDC-TR-78-6B 


Case  3 - The  z axis  was  rotated  5 deg  relative  to  the  nozzle 

centerline.  In  this  case,  the  flow  in  the  starting  plane 
is  nonuniform  and  was  obtained  by  interpolation  from  the 
flow  field  predicted  for  case  1 , 

The  predicted  wall  Mach  number  distribution  for  case  1 is  shown  in 
Fig.  5 along  with  the  predictions  of  the  Lockheed  axisymmetric  MOC 
computer  program  (Ref.  6)  which  is  well  verified  and  widely  used.  The 
results  from  the  two  programs  are  in  good  agreement.  At  the  nozzle  exit 
plane,  the  circumferential  variation  of  wall  Mach  number  calculated  by 
the  3-D  MOC  program  is  much  less  than  one  percent. 


Figure  5.  Results  for  z axis  aligned  with  nozzle  centerline. 

The  results  of  the  three  3-D  computations  are  compared  in  Fig,  6 
where  it  is  seen  that  the  three  predicted  wall  Mach  number  distributions 
are  essentially  identical. 

These  computations  of  axisymmetric  nozzle  flow,  done  the  "hard  way" 
with  the  3-D  MOC  program,  indicate  that  the  program  is  acceptably  accu- 
rate and  is  capable  of  computing  the  supersonic  flow  in  complex  configu- 
rations . 
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Figure  6.  Results  of  3-D  MOC  program. 


In  each  of  the  test  cases,  346  streamlines  were  traced  throughout 
the  flow  field.  For  cases  1 and  2,  78  reference  planes  were  computed 
which  required  a CPU  time  of  approximately  25  min.  Fewer  reference 
planes  were  computed  for  case  3,  with  a corresponding  reduction  in  CPU 
time. 


5.0  CONCLUDING  REMARKS 

The  results  obtained  for  an  axisymmetric  nozzle  with  the  3r-D  MOC 
computer  program  are  acceptably  accurate.  However,  additional  verifi- 
cation of  the  program  is  clearly  required.  The  predictions  of  the  code 
should  be  compared  with  experimental  results  for  3-D  supersonic  flows, 
whenever  adequately  detailed  experimental  results  become  available. 

Two  extensions  of  the  3-D  MOC  computer  program  are  recommended. 
First,  the  program  should  be  modified  to  include  a constant-pressure 
free  boundary  condition  for  part  of  the  flow  field.  This  modification 
would  permit  the  computation  of  the  exhaust  plume  from  a 3-D  nozzle. 

The  second  recommended  modification  is  the  inclusion  of  shock  formation, 
both  within  the  nozzle  and  in  the  exhaust  plume. 
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APPENDIX  A 

DESCRIPTION  OF  COMPUTER  PROGRAM 


DESCRIPTION  OF  ROUTINES 


NAME 


FUNCTION 


/ 


MAIN  CALLS  CARDIN,  INGEOM,  INPUT,  AND  THREED 

BOYFIT  CURVE  FITS  BODY  IN  AREA  OF  INTEREST 
BODY  BODY  POINT  ROUTINE 

CARDIN  READS  AND  PRINTS  INPUT  CARD  IMAGES 

CUT  INTERPOLATES  THE  BODY  GEOMETRY  IN  CURRENT  PLANE 

DEi-TAF  COMPUTES  THE  B (CHARACTER  1ST  IC  ANGLE  DELTA 

01  ST  COMPUTES  MIN,  DI ST.  OF  CHARACTERISTIC  INTERSECTIONS 

FIELD  FIELD  POINT  ROUTINE 

FIND  LOCATES  POINTS  NEAREST  TC  GIVEN  POINT 

FINDBP  LOCATES  BODY  GEOMETRY  POINT  NEAREST  GIVEN  POINT 

FIT  SURFACE  FIT  ROUTINE 

IDENT  WRITES  HEADER  AND  TRAILER  LABELS  ON  PLOTS 

INGEOM  READS  GEOMETRY  INPUT 

INPUT  REAOS  STARTING  PLANE  INPUT 

NEIGH  LOCATES  (OR  READS)  8 NEIGHBORS  OF  EACH  INPUT  POINT 

NEW? AP  NEWTON-RAPHSCN  ROUTINE 

NORMAL  COMPUTES  THE  NORMAL  VECTOR  TO  THE  BODY  SURFACE 

OUTPUT  CONTROLS  PRINTED#  PLOTTED,  AND  TAPE  OUTPUT 
SI  MO  SOLVES  SYSTEM  OF  SIMULTANEOUS  LINEAR  EQUATIONS 

SOLVBP  SOLVES  FDR  THE  LOCATION  OF  A BODY  POINT 
SORT  SORT  ROUTINE 

THREED  CONTROLS  THE  3-D  CALCULATIONS 
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INPUT  NOMENCLATURE 


G RATIO  OF  SPECIFIC  HEATS 

HEAD  PROBLEM  DESCRIPTION 

(OBUG(I)  SPECIAL  PRINT  FLAGS  =0  NO  PR  I NT  I *1  PRINT 
IDBUG  ( 1)  PRINT  INPUT  GEOMETRY 
IDBUGl  2)  PRINT  FINAL  GECMETRY 
IOBUGC  3)  PRINT  INITIAL  PLANE  INPUT 
I DBUG ( 41  PRINT  FINAL  INITIAL  PLANE  DATA 
I D6UG < SI  PRINT  NEIGHBORS 
IDBUG(ll)  PRINT  NORMAL  CALCULATION 
IMV  INITIAL  PLANE  DATA  TYPE 

INV-0  Q THETA  PS  I 

I MV=1  M THETA  PS  I 

IMV-2  to  U V 

INEIGH  NEIGHBOR  CONTROL 

INEIGH=0  NEAREST  POINTS  ARE  NEIGHBORS 

I NE IGH=t  SPECIAL  CCMPUTATICN  OF  NEIGHBORS 

INEIGH-2  READ  NEIGHBORS  FROM  CARDS 

IF  ITYPE=0  AND  INEIGH=0*  INEIGH  IS  SET  TO  I 
I PL OT  *0  NC  PLOTS:  *1  PLCTS 

IT  III  COORDINATE  SYSTEM  FOR  BODY  POINTS  AT  THIS  STATION 

IT{IJ=0  RECTANGULAR  COORDINATES 
IT(I)=1  PCLAR  CCCRDI NATES 

I T( I 1*2  PCLAR  CCCRDI NATES  - AXISYMMEHIC 

I TYPE  COORDINATE  SYSTEM  FOR  STARTING  PLANE  INPUT 

I TYPE  =0  PCLAR  CCCRDI  NATES  - AXISYMMETR1C 
ITVPE  = 1 PCLAR  COORDINATES 
ITVFE=2  RECTANGULAR  COORDINATES 
IT  11  FLOto  FIELD  TAPE  OUTPUT  IT  11=0  NG  TAPES  IT  1 1=1  TAPE 

7 T 1 2 BODY  DATA  TAPE  OUTPUT  IT12=0  NC  TAPES  I T 1 2=  1 TAPE 

I\f{ I ) NO.  CF  BODY  POINTS  AT  THIS  STATION 

JMAX  MAXIMUM  NO.  OF  PLANES  TO  BE  CCMPUTEO 

JPT  PRINT  EVERY  JPT  PLANES 

N9T  =0  INPUT  FOR  FIELD  POINT J =1  INPUT  FOR  BODY  POINT 

NEI J.Kl  THE  J-TK  NEIGHBOR  OF  THE  K-TH  POINT 
NPTS  TOTAL  NO.  OF  POINTS  IN  STARTING  PLANE 

IF  ITYPE=0  NO.  OF  POINTS  ON  A RAY 
NR  AYS  NO.  OF  RAYS  IF  ITYPE=0J  OTHERWISE  NOT  USED 

NF  NO.  OF  STATIONS  OF  BODY  GECMETRY  INPUT 
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INPUT 


PC  I J 
PSK  I J 
P.'0<  I > 
PTO 
Q ( I 1 
RO 

THETM 
TTO 
X(  I ) 

XCC  J*  I 
Y { I » 

YC< J.I 
ZC(  I) 
ZMAX 

zo 


NOMENCLATURE  (CONT*» 


PRESSURE 

ANGLE  SHOWN  IN  FIG*  1 i OR  V IF  IMV=2 

TOTAL  PRESSURE  AT  POINT  CIF  0*0*  SET  TO  PTOJ 

TOTAL  PRESSURE  FOR  IRROTATICNAL  FLOW 

VELOCITY  IF  IMV=0;  MACH  NUMBER  IF  IMV=1;  W IF  IMV=2 
GAS  CONSTANT 

I)  ANGLE  SHOWN  IN  FIG*  1 ♦ OR  U IF  IMV=2 
TOTAL  TEMPERATURE 
POINT  IN  STARTING  PLANE 

X IF  ITYPE  = 2 * R IF  ITVPE=0  OR  I 
> BODY  POINT  X IF  IT1I)=0;  R IF  IT(t)=l  OP  2 
POINT  IN  STARTING  PLANE 

Y IF  I TYPE =2  1 ANGLE  IF  ITYPE=0  OR  1 
I BODY  POINT  Y IF  IT(I>=0;  ANGLE  IF  IT  II  1=1  OR  2 
Z LOCATION  OF  BOOY  STATION 
MAXIMUM  Z LOCATION  FOR  CALCULATIONS 
Z LOCATION  OF  STARTING  PLANE 
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INPUT  CAROS  IN  ORDER  READ 


CARD 

FORMAT 

CONTENTS 

1 

( 20A 4 ) 

(HEAD( I ) .1 *1 .20 ) 

2 

(8011) 

(IDBUG(I) ,1=1,80) 

3 

(IS) 

NS 

4 

5 

(215.E10.0) 

( 8E  1 0 « 0 ) 

READ 

IV(I  > , I T(  I ) ,ZC( I I 
( XCC  J,I ) «YC( J • I ) ,J=1 .IV(I>) 

CAROS  A AND  5 NS  TIMES 

€ 

(4EI0.0) 

G.PTO.TTO.RO 

7 

( 2E  10. 0.1  01  5) 

ZO  *ZMAX , I TYPE  « I MV*  NPTS • NR AYS • JMAX  * 
JPT • I PLOT  * I Tl 1 * IT  12,1 NEIGH 

a 

< 7E 10.0*15) 

X(  I ) ,Y(D,PCI),Q(1),ThETA(1),  PS  t ( II, 
PTG ( I > ,NBP READ  NPTS  CAROS 

9 

(9110) 

READ 

(K,(NE( J,K),J=1,8) .1=1, NPTS1 

CARO  9 ONLY  IF  INEIGH=2 

THREAD 


*****  PRIHARY  OVERLAY 

secondary  overlay 
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C 

THREE  OIMENSIONAL  DUCT  FLOW  SRR00886 

c 

SIMPLIFIED  VERSION  WITH  FOLLOWING  ASSUMPTIONS 

C 

NO  INTERNAL  SHOCKS 

c 

IDEAL 

GAS 

C 

UNIT 

FUNCTION 

c 

5 

CARD  INPUT 

C 

6 

PRINTED  OUTPUT 

c 

10 

PLOTTED  OUTPUT  > 

c 

11 

FLOW  FIELD  OUTPUT 

c 

12 

BODY  OUTPUT 

c 

8 

WORK  DATA  SET 

c 

OVERLAY  STRUCTURE 

c 

OVERLAY 

A 

c • 

INSERT 

CARD1  N.INGEOM.I NPUT. NEIGH, FIND. SORT 

c 

OVERLAY 

A 

c 

INSERT 

THREED, CUT .F IT, DIST ,S IMO.OELTAF 

c 

OVERLAY 

B 

c 

INSERT 

FI  ECO 

c 

OVERLAY 

B 

c 

INSERT 

BODY .FtNDSP.SOLVBP. NORMAL, BDYFIT, NEWRAP 

c 

OVERLAY 

B 

c 

INSERT 

OUTPUT 

COMMON  /TITLE/HEAOI20I 
COMMON /I  DEBUG /I  DBUG  ( 80  > 

COMMON  SSTAG/DUMMYf 1 207 ) • I PLOT .IT11.IT12 
DATA  ISTRP^O/ 

1 FORMAT C 20A4 ) 

2 FORMAT C 801 1 1 
CALL  CARDIN 

20  CONTINUE 

READ!  8.1  «END  =99  I HE  AD 
READ! 8.2> I DBUG 
CALL  I NGEOM 
CALL  INPUT 

IF f I PLOT ,EQ. Q1  GO  TO  21 
IF! ISTRP.NE.OI  GO  TO  21 
I STRP^l 
CALL  IDENTIM 

21  CONTINUE 
CALL  THREED 
GO  TO  20 

99  CONTINUE 

IF1 ISTRP.NE.O)  CALL  IOENTI21 

STOP 

END 


28 


AEDC-TR-78-68 


SUBROUTINE  BDYFITI A .B.C .XB.YB. I * IT) 
IMPLICIT  REAL *8 ( A— H * 0— Z) 

REAL *4  DUMMY! 5051 ) .XB150) *YB(5Q ) 
COMMON  /GEOM/OUMMY*NP 


c 

I T=1 

X=C 

c 

IT=2 

Y=C 

c 

I T=3 

AX+BY=C 

c 

I T=4 

( X— A) **2+< Y—B) **2=C 

I Ml -I— 1 

IF(I.EO.l)  I Ml  = NP 
XlsXB(IMl) 

Y1=YB< I Ml) 

X2=XB< I) 

Y2=YB ( I ) 

IPIssI  + 1 

IF(l.EQ.NP)  I PI =1 

X3=XB( 1P1) 

Y3=YBC I PI ) 

IF  M XI *NE • X2  ) • OR* ( XI  * NE.X3 ) ) GO  TO  21 

IT-1 

A — 1 • 0 

B-0-  0 

C=X2 

RETURN 

21  CONTINUE 

IF« t Y1 .NE.Y2J .OR, IY1.NE.Y3) ) GO  TO  22 

I T=2 

A-0.0 

B = 1 .0 

C-Y2 

RETURN 

22  CONTINUE 

R 1=X1 **2+ Yl **2 
R2=X2**2+Y2**2 
R3=X3**2*Y3**2 
OX2I— 2 • 00* I X2— XI ) 

OY21 =2*00*1 Y2-Y1) 

0X32=2.00*1  X3-X2) 

OY32=2.00*1 Y3-Y2) 

IF10X21*  NE* 0*00)  GO  TO  23 
0X21=2.00*1  X3-X1) 

D X32=— 0X32 

0Y2 1=2. 00*1 Y3-YI ) 

DY32=— D Y32 

H=R2 

R2=R3 

R3=H 

GO  TO  24 

23  CONTINUE 

1F(OX32.EO*O.DO)  GO  TO  24 
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IF ( [ D Y21  /DX21  ) • NE  • ( 0Y32/DX32 1 1 GO  TO  24 
l T = 3 

A=-DY2 1/0X21 
8=1.00 
C=Y1+X1*A 
RETURN 

24  CONTINUE 
IT=4 

A=(R2-R1 1/0X21 
B=R3-R2 
C=0  Y21/DX21 
H=0 Y32 

IF( DX32.EQ. 0.00)  GO  TO  25 

S=B/DX32-A 

H=H/DX32— C 

25  CONTINUE 
8=0  /M 
A=A— B*C 

C=R1-2.00*<  A*X1  *B*Y  1 ) 

C=C4B4*2+A**2 

RETURN 

END 
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SUBROUTINE  BOOYC t ) 

IMPLICIT  REAL«8(A-H.G-Z> 

REAL *4  ZD.XNC 1000) »YN(  1000) • PN( 1000 >• RN( 1 000  >. QNC 1000) • 
•PSNf 1000) .ZE.XC1000) .Yd  000) .PC  1000) .RHO(IOOO). 

*QC 1000) .THETA11000) .XBI50) >Y8 (50  ) .X80C50  ) • YBOC  50  ) • 

, ♦TN( 1000) .PSI(lOOO) 

COMMON  /N  VALU/ZD  <XN »YN  « PN  «RN*  ON.TN. PSN 
COMMON  /VALUE/ZE  .X  *Y «P .RHO • O.THETA.PS l . N 
COMMON  /BCUT/XB  # YB 
; COMMON  /OCUT/XBO ,Y8C 

COMMON  /DGAS/DZ.G.RO.GMl  ,GM1H.  GMIG . GGM  1,  GP1  • GPGM 
C OMMON/I DEBUG/I DBUGt  80 ) 

REAL *8  Nil «N21 *N31 *N1 2 * N22 • N32 • N 1 ■ N2 . N3 
COMMON  /POI  NT/X1  .VI  .Pi  .R1  .01  ,T1  .THl.PSlt BE1  .X2.Y2.P2. 
4R2.Q2.T2 « TH2  »PS2  *BE2  *P2C .T20.PSQ .DEL (4 ) *PI(4)*THl ( 41 » 

*P 11(4) .BE 1(4) . X I ( 4 ) • YI ( 4 ) .OEI (4}.DX(4)*DY(4). DC C 4 ) ■ 

*PP< 5.4) .T12.P12.B12 . T2 1(4) ,P2 I (4 ) .82 1 ( 4 >. D2I (4 ) 

COMMON  /FLDBD  Y/S I NT 12 .COST  1 2 .SINP1 2. CDS PI 2*  S INB12. 

4CQSB 1 2 «S I NO 21 (4 ) .COSD2I (4 » .S INT2 I (4 ) . COST  2 U 4 ) • 

4S1NP2I (4) .COSP2I (4) .Si NB2 1(4) . COS  82  I (4) . OX ON (4  I. 
40Y0M4)  «DZDN(4)  .OTOXC4)  .DTDY  ( 4 ) .OTOZC4)  .OPOX(4), 

♦DPDYt  4) . DPDZC 4) «OER(2 .5 ) • AL( 3 .4 ) • CR(4 ) . AM ( 3. 3) • CM ( 3. 4 ) 

1 FORMAT ( • 0 FAILED  IN  BODY  AT  POI NT ■ • (5 . 5X t * Zl- * . 

4 1PE12.5.5X*  *DZ=*.E12 .5//5X • »X1 =• .El 2 .5.5X. • Yl  = » . El 2. 5. 

♦ 5X.  *Plss».El  2.5  .5X.»T1  = »,E12.5  ,5X  . * PS  1 = « * El  2 • 5/  5X*  • X2=  * 
4.E12.5.5X. * Y2—  • .E12.5.5X* ■P2=*  « El 2.S .5X • • T2*1 , El 2. 5. 5X 
*• *PS2—  • »E 1 2 • 5/45X . ' P2  0= * .E12.5.4X. * T 2 0= * . E12.S. 4X. 
*«PS20=»  *El  2 * 5//) 

X1  = X(I  ) 

Y1=Y( Z ) 

Pl-PC I ) 

R !=RHO( I ) 

Q1=Q(I> 

THl=sTHETA(  l ) 

PS1  =PSI  ( I ) 

T1=P1/(R0*R1  ) 

SI  NT1— DSI  N(  TH1  ) 

COSTl-OCOSITHl ) 

SI  NP  I=DSI  N<  PS1  ) 

COSP1 =DCOS( PS) > 

8E l—DARSI N( D SORT ( G 4P1/R1 ) / 01 ) 

P2-P1 

R2=Rl 

Q2=Q1 

TH2=TH1 

PS2=PS 1 

T2  = Tl 

BE2=BE1 

CALL  F INDBPI XI *Y1 .XBO.YSO. IBP1 ) 

CALL  BDYFITCAl »B1 .Cl .XBD.YBO. IBP1.IT1) 
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CAUL  BDYF  I TI  A2  «B2  *C2  .XB  *YB*  IBP1  . IT2) 

CALL  NCRMALIXl *Y1 • ZE • A1 *B1 *C 1 t Nl 1 * N21 • N3 1 • IT  1.  IBP 1 ) 
CALL  NOR MALI XI *V1 • ZE+OZ • A2 *02 *C2* NI2 . N22. N32. IT2. IBP1 > 
OEL(  l)=DARCOSI-Nl  1 *SI  NT1  *COS PI  *N21 *COST 1-N31 *S  INTI* 
*SINP1  > 

IFfN31.LT.  0.001  OELI  1 1=6  *28 31 8 5 30 7 DO— DEL  III 
DEH  2>=0ELI1)-1 .57079632700 

1FIDELI  2J.LT.O.OO)  OELI 2) ®6 • 203 1853Q7 DO  + DEL (2 ) 

DELC  3 ) =DEL< 1 1 + 1 . 57079632700 
IF (OELI 3 >.GE.6« 203 16530700) 

♦DEL ( 3 ) *DELt  3 1—6.2831 053  0700 
DO  21  XT®I.4 
THI < IT) = THl 
PIK  IT)=PS! 

BEK  IT)=BE1 
DE1( [ T) =DEL( IT) 

21  CONTINUE 
IBP2=IBPI 
V2=Y1 

DO  25  I T=1  *25 
T12=0.  500*1 TH1+TH2) 

P 1 2—0. 5D  0*1 PS1*PS2I 
8 1 2* 0.500*1 BE I+BE2) 

SI NT 12*0 SI N( T12) 

C0ST1  2=QCOS(T12) 

5INP 12=0 SI N(P12> 

CO  SP 1 2=DC0St  PI  2 ) 

SI N8 1 2*0 SI N(812 I 
COSB 1 2=DC0S( 81 21 
Nl  =0*  500 *(  Nl  1 *N12  I 
N2=0* 500*1 N2 !♦ N22 ) 

N3=  0. 50  0*1  N31  + N32) 

CC  1=N3*C0ST1  2+C0SP1  2-N1  *C0ST1 2*S  1 NP1  2 
CC  2=N1 *SI NT1 2— N2*C0ST1 2 *C0SP1 2 

CC3=( N3*SI NT12— N2*COST12*SINPl2 I *02*CC1*X 1 +CC2*Y1 

Y20®Y2 

IBP® IBP2 

X20=X2 

CALL  SOL  VBP( X2 *Y2  *CC1  *CC2 .CC3. A2 . 82. C2. IT2.X20.Y20> 
CALL  CKB0YPIX2.Y2fXB.Y8.IBP2) 

IF( IBP.NE.IBP2)  CALL  BOYF I T ( A2 ■ B2 • C2.XB. YB. IBP 2. IT2) 
CALL  NORMAL  I X2.Y2.  ZE+OZ.  A2.B2  .C2  • Nl  2.  N22,  N32  • IT2.  IBP2) 
SINT2=OSIN(  TH2) 

CO  ST2=DCQSI TH2 ) 

SI  NP  2®0  SI  NIPS2) 

COSP2=OCOS< PS2) 

DELI  1) ®DARCOSC-NI2*SlNT2*COSP2*N22*COST2“N32*S INT2* 
♦SINP2) 

IF ( N32.LT.O.OO)  OELI 1 ) =6 .28318S307O0-0ELI 1 ) 

DELI  21 *DELll >-1.57079632700 
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lF(OEL(2) .LT. 0.001  OEL ( 2 1 =6 . 283 1 8530700 +DEL (2 > 

OELC  3 ) =OEL(  1 1 +1 « 57079632700 
IF(DEL<  3) .GT. 6. 28318530700) 

*DEL{  3)  =DEL(  31 -6.2831  8530700 
DO  23  J-I  .3 

T2IC JI=0.5D0*<TH2+THI< Jll 
P2IC  J) -0.500*1  PS2*PI  1 I Jl) 

B2K  Jl=0* 500 *( BE2 ♦BEl ( Jl > 

021  C J»=0.  5DO*COEL<  J)  fOEKJ)  ) 

02  I ( J) -021 ( J )-3. 1 41 592654  00 

IF ( D 21 ( J) .LT.O.OO ) D2I ( J) -6. 2831 85307D0+D2I ( J ) 

SINT2I  ( J)  -0S1NIT21  ( J)  ) 

C0ST2I  ( J)=DC0SIT2I  I Jl  ) 

SINP2I  I J1=DSIN(P2I  C J)  J 
C0SP2I ( Jl »DC0S(P2 1(J)  1 
SINB2I  ( J)=OSIN(B2I(  J)  ) 

COSB2I ( J) =0005(02  It J>  1 
St NO  21 ( J)=DSI NC02I (Jl) 

COSD2E ( J)=DC0S402I< J) ) 

F1=SIN82I  ( J14C0502K  Jl 
F2-S1NB2I ! J) 4SIND2I ( J) 

F3=C0SB2I  ( J14C0ST2K  Jl-Fl  4SINT2M  Jl 
OLl  J) *DZ/f  F34C0SP2 1 ( J)  — F24SINP2  1 ( J > 1 
0X<  J)  = SINT2I J Jl 4COSB2 1 ( J ) +F1 4C0ST2 I ( J 1 
DY( J ) — F3 *S I NP2I C JJ+F24C0SP2I ( Jl 
XI  C J)  = X2-DX(  J)  *DL<  J) 

YK  Jl=Y2-0YI  J)  «DL(  Jl 

CALL  CETPTCXKJl.YK  JI.PPII.J1) 

RA0=0SQRTCG*PP<1 .J1/PPC2. Jl I/PP13.J1 
1FC (RAO. GE. 0.00) • AND. ( RAO.LE. 1 .00 ) 1 GO  TO  22 
WRI TEC  6.2)1  ,J.P1 *R1 .Q1.TH1 .PSl.eEl.P2.R2.Q2.TH2.PS2. 
4BE2.P20.R2a.Q20.T2  0.PSO.OELU  1 .DEL  (2  1 .OEL 13  1 *N1 1.N21. 
4N31.N12.N22.N32.N1 . N2 . N3 • OLC J) • OX C J 1 *0Y I Jl . OEI ( J 1 • 

♦PPt 1. Jl »PP(2. Jl .PP(3 . Jl *PP(4. Jl .PPI5.J) .THt I Jl .PI ICJ I. 
*QEU  J) 

2 FORMAT {*1  FAILED  IN  BODY  AT  I=* . I3.5X* • J=* • 13// 
4IP6E14.5/6E1 4.5/SE14.5/3E14.5/9E14.5/4E14.S/8E14.5I 
STOP 

22  CONTI  NUE 

CALL  DELTAFt  J.DELl  J}  .DEI  ( J)  1 
1F(DA8S(0EI( Jl ) *LE. I .0-81  DEI ( J)=O.DO 
CALL  GETOERt XI ( Jl «Y1 ( Jl .DERI 
PI  C J1=PPC 1 « Jl 
THIC  J»*PP«4.J» 

PU  t J)=PPI5*J) 

BE I ( Jl =DARSIN(DSQRT<G*PP< 1 ■ J1/PPC2 . J 1 1/PPC3. J I I 
T2I(  J)=0.5DQ*(TH2+THI(  J)) 

P2I<  J)=0.500*CPS2+PI  I<  Jl  1 
B2I ( Jl *0.50041 8E2+8EI I Jl I 
D2U  J)  =0.500* (OELC  Jl+DEH  J)  I 
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02  U J)=02l< Jl— 3, 141 592 6 5 A DO 

IF(02I( J) .LT.O.DO)  021 C J) =6.2831 8S307DO*D2I<J I 
SINT2K  J>=OSIN(T2I<  J»  ) 

COST2I4 J)=DCGS4T21< Jl  I 
5INP2I f JJ  =OSI  N4P2II  Jl  ) 

C0SP2K  Jl=DCOS(P2l4  J)  ! 

SINB2I  I J)=D5JMB2It  Jl  I 
COS32K  J)=DCGS4B2IC  J)  ) 

SIND 21 ( J 1*0 SI NCD2IIJ1 ) 

CO SO 21 4 J I =DCOS( 02  I 4 J)  ) 

DTOXl  J I =D£R (1*41 
OPOXC J»=DERf |,S) 

OTOYC J)=DER42.41 
DPOYI J)=DER42,5) 

Fl=COST2I C J> 4SIND2I 4 J» 

F2=SINT2If  J)«S1ND2](  J» 

F3=S!NB2I ( J| *DLf  J) 

F*=COST2I 4 J) +C0SD2 1 4 J) 

0XDN4  J l'*—F  I 

OYON4  J1=F2*SINP2I  ( J) ♦COSP2 I 4 J I+COSD2  1 4 <J) 

0ZDN4  J l=F2*COSP2 1 C JI-SINP2I 4 JI*C0S02 14 Jl 

F5=XIf  JI-X2 

F6=YI4J)~Y2 

DTD24  J)~(0T0X4  Jl  ♦FS+-OTOY 4 J > *F6 1/DZ 
OPOZf Jl=f DP0X4 J)*F5*DPDY4 J)*F6)/DZ 
AL4  l.J)=COSB2If JI/CSINB2II J)*PPf 2,J)*PPf3, J>**2) 

AL4  2,  J)=COS02I  4 Jl  -F3*SI ND2  I C J) *DZDNI  J l/DZ 
AL( 3*J1=F1*F3*F»*DZDNC Jl/DZ 

CRC Jl=F3*f  SIND2I 4 Jl *<DTDX ( J 1 *DXDNC  J I FDTDY 4 J I *DYDN 4 J I * 
•DTDZt J ) ♦0Z0N4 Jl I — F444DP0X 4 JI*OXDN<  J 1 *DPDY  4 J ) *OY ON 4 J »♦ 
♦DPDZf  J I *DZON(  J)I1  +ALI  1»J1*PI(J1  FALf  2 • J 1 4THI  (JI+ALI3.J) 
**PIIC Jl 
23  CONTINUE 

A1=AL(2.1 >*ALC1 ,21-ALf 1»1I*AL(2*2) 

B 1=A1_{  3«  1 1 *AL<  1,2  l-ALf 1 ,1  l*AL(3,2l 
C1=CR<  1 1 *At_(  1 ,21-ALI  1,1)  *CR42) 

AMI  2,1  l-=TH2 
AMf  3,1  )=PS2 

CAUL  NEWRAPfAMf2.il «AM(3  *11 • A1 • Bl *C l • N1 2 * N22.N32 I 
A1=ALI 2*2 1 *AL( l ,3 1 — AUl 1,2) *ALf2,3 1 
B1=AL(3«2)*ALC1,3)-ALI1 ,2>*AL<3,3) 

C 1 =CR 4 2 1 *AL4 1 ,3 1- ALC 1 ,2 1 *CR 43  I 
AMI  2 *2 )=TH2 
AMf  3.2I-PS2 

CALL  NEWRAPf  AM4  2 ,2) * AM43.2 1 , A1 *B1 • C 1 • N1 2.N22 ,M32 1 
Al=ALl2.3)«ALf l.ll-ALfl »3)*AL42,1 I 
B1-AL4 3,3)*AU( 1,1 I-AL4 1 *3) *ALf3 ,11 
Cl=CR4 3>*AL<1 ,1)-ALC1 ,3)«CRfl) 

AMf  2 ,3)=TH2 
AMC  3 ,3)=PS2 
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CALL  NE WRAP ( AM  f 2 *3) .AN ( 3 .3  > . A 1 . B1  ♦ C 1 , N 1 2 . N22. N32  > 

AM<  X.1)=<CR(1  )-AL(2*l)*AM(2*l )-AL(3«l)*AM<3*  1)  )/AL(l*l> 
AM< I ,2>=ICR(2)-ALC 2,2)*AM(2,2>-AL(3.2)*AM<3.2> )/AL< 1*2) 
AM< 1 .3l=ICRf3J-AL<2.3>*AMI2»3)-ALC3.3>*AH|3.3> )/AL(l.  3) 
P2>0.500*(  AMI  l *1  ) + AMI1  <3)1 
TH2=0.5D0*f  AM<2.1  ) + AM<2.3>  » 

PS  2=0.  500*C  AMI  3*1  ) +AMI3  *3  ) ) 

T2=TT0D*I  P2/PTG0 (1 ) ) **GH1 G 
R2-P2/CRGAT2) 

a2=DSQRT<2.00*GGNl*RO*(TTOO-T21> 

IFC 1T.EQ.1)  GG  TO  24 

IFCDABSU  P2-P20I/P2I  .LE.  1.0-7)  GO  TO  26 
IFIKT.EQ.25)  GG  TO  25 

24  CONTINUE . 

P2D=P2 

R 20— R2 
Q 20-02 
T20=TH2 
PSO=PS2 

BE2— DARSI NIDSQRTIG4P2/R2I/Q2) 

25  CONTINUE 

IFf DABSC|P2-P20>/P2I.LE*l*D-5)  GO  TO  26 
WRI TE( 6 .1 )I *ZE *0Z  .XI *Y1 .PI *TH1 •PS1.X2.Y2, P2.TH2.PS 2. 
♦P20.T20.PS0 

IF  ( 0ABS4 ( P2— P20I/P2 ) »GT • 1 .0—4}  STOP 

26  CONTINUE 
XNI I ) = X2 
YNI I )=Y2 
PNC  I )=P2 
RN( I ) ~R2 
ONI  1 ) =Q2 
TNI  I 1 = TH2 
P SNI  I ) -PS2 
RETURN 
END 
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SUBROUTINE  CARDIN 
DIMENSION  A(20) 

1 FORMAT! 20A41 

2 FORMAT!  5X  *2  0A4 ) 

3 FORMAT!  *1  * *38X. 'INPUT  CARDS*} 

4 FORMAT!  • 0 * .4  X .9!  *0  • 1 .1 0 ! • 1 • ) .10  < • 2 * 1 • 10  ! * 3*  I » 101*  4*  ) , 
*10!  *5*}  .1  0!  *6  *1  .10!  *7*  > .*8  */5X.Bf*  1234567890*  1/1 

IP=5Q 
WRITE!  6.3) 

WRI TE!6«4) 

20  CONTINUE 
READI5.1 ,ENO*99)A 

IF! IP.NE.O)  GO  TO  21 

WRITEI6.41 

WRITE!  6*3) 

WRITE! 6,4) 

1Pa50 

21  CONTINUE 
WR] TEC  6*2 1 A 
WRI  TE ( 8*1  ) A 
IPalP-l 

GO  TO  20 
99  CONTINUE 
REWIND  8 
WRITE! 6*41 
RETURN 
END  ' 
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SUBROUTINE  CUT (ZB) 

COMMON  /BCUT/X8(50> •YB<50> 

COMMON  /’GEOM/ZCI  50  ) *XC  C 50  *501  * YCC50  • SO  ) , NS  * NP 
REAL *4  A ( 3 • 3) 

z-ze 

DO  21  1=2, NS 

IFCZB.LE,ZC( f) ) GO  TO  22 
2t  CONTINUE 
I=NS-1 

22  CONTINUE 

IFI  I #EQ,NS)  I =NS— 1 
IP=I+l 
I M=I —1 
ZO=ZC  ( I M) 

Z1=ZCC I ) 

Z2=ZC< IP) 

Z 10=1.0/C Zl-ZO) 

Z 20=  1 ■ 0/<  Z2— ZOI 
Z2l=1.0/(Z2-Zl» 

ZPl=ZO+Zl 

At  3, 1 >=Z2l*<  Zi 0-220) 

At  3,2)=-Z21*ZlO 
A ( 3,3)=  Z21*Z20 
A(  2,l)=-Z10-ZPI*A(3,l ) 

A(2,2)=  Zl 0— ZP1 *At  3*2) 

AC  2,3)=  — ZP1 *A (3,3) 

At  1 , 1 )=— ZO*  t At  2,1  ) + At  3 « I )*ZO)«-1*0 
At  1 , 2)=— Z0*1  AC  2«2)+At3,2)  *Z0) 

AC  1 , 3 ) =— ZOAt  AC  2,3)  + AC3»3)  *Z0) 

DO  23  J=1 , NP 

B=A( 1,1| *XCC J *1 M) FAC  1 ,2) *XCt J,I >+A< 1 ,3)*XCC J, IP) 
C=A(2.1)*XCI  J,1  M)+Af  2,2>*XCt  J. I >+A<2 .3)*XCt J . I P) 
D=A<  3,1>*XC<  J.IM)+A(3,2>  *XCC  J.I  )+A(3,3)*XC<  J.IP» 
XBI J)=B+CCFD*Z) *Z  t 

B=At I • I ) *VCI J*IMI«-A(1 ,2) *YCt  J, I )*AC1 »3I*Y  Ct J, IP) 
C=A(  2,1)*YC(  J,  t H)  4-AC2.2)  *YCt  J, I > +AC 2, 3) *YC t J,IP) 
0=A(3,  1)  *YCC  J,  I MM- AC  3,2)  *YCC  J,I  )*A(3,3 >*YC( J.  IP) 
YBt J)=B+CCFO*Z| *Z 

23  CONTINUE 

RETURN  S" 

ENO 


\ 
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SUBROUTINE  OEt-TAF  I J .DEL  t DELlI 
IMPLICIT  REAL *81 A-H.O-ZI 

COMMON  /F LOBOV /DUMMY ( 1 4 ) *S i NT  I ( 4 ) .COST 1(4 I.SINPII 4I« 
4C0SPI 141 ,SI N8I (4) .CCSBI (4) 

A— SI NT 1 1 J) 4COSPI I J>  *SI NBl ( J> 

B = SINPI( J)*SINBIC J| 

C=-0. 5DO*COSTI  ( J)  4C0SPI  ( Jl  *(SINBI  ( J)  t-COSBI  ( J > J 
SI  ND=D SI NC DELI 
CO  SD=DCO  S ( DEL } 

OELI=DEL 
ASBS=A**2+B**2 
RAD=AS8S-C**2 
IF1RAO.LT. O.OO)  RETURN 
RAO-OSQRT(RAO) /AS8S 
SINI=B*C/AS8S 
SINI  l-SINt>A*RAO 
SINI2-SINI-A4RA0 

IF ( OABSI SI ND— Si NI 2 » • LT . DABS IS  I ND— S INI 1 ) ) GO  TO  21 
SI NI =Sl NI 1 

IFCDABSI  SINI  I.LE.1.D0I  CO  TO  22 

21  CONTINUE 

S I NI =SI NI 2 

IFC OABSI SINI I.GT. I.OOI  SINl=SCNll 
IF(  DAB  SI  SINI  I.GT.  1.001  RETURN 

22  CONTINUE 
COSI=A*C/ASBS 
COSI l=COSI*B*RAD 
COSI 2=C0SI-B*RAD 

IF  1 OAB  St  COSD— COSI  2 ) . LT ■ DABS  ( COSO— COS  III}  GO  TO  23 
COSI^COSI  1 

IF I OABSI COSI ) .LE.l .DO)  GO  TO  24 

23  CONTINUE 
COSI»CGSI2 

IFI OABSI COSI I.GT.l.DO)  COSI=COSlf 
TF( OABSI COSI I.GT. 1.00)  RETURN 

24  CONTtNUE 
0ELI=0ATAN2( SINI .COSI » 

IF < DEL I .LT .0.00)  OELI-6. 2031853O7DO+OELI 

RETURN  / 

END 
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SUBROUTINE  01  ST  I DELS  ) 

COMMON  /GAS/G.RO .GM1 .GM1 H. GM1G. GGM1 • GP 1 • GPGM 
COMMON /VALUE /ZE* XI 1000) .V(1000) .PI 1000 >• RH0< 1000 > . 
♦0( 1000) • THETAI 1 000 1 .PSI (1000) .N 
I NTEGER *2  NE(8«1000) 

COMMON  /NEAR/NE.DMIN 

OMIN=1.E10 

QMI N— 1 .El  0 

AMA  X-0 . 0 

DO  22  1-1  . N 

A— P( I ) /RH0( I ) 

IF(A.GT.AMAX)  AMAX-A 
1FIQII I.LT.OMIN)  QMI N— Q( I ) 

00  21  J«1  .8 

K-=NE  { J * I ) 

D = t XI  I )-X(K})  **2+( Y(  I l-Y(K)  1 **2 
IF  ( D.LT.OMI  N>  DMI  N=D 

21  CONTINUE 

22  CONTINUE 

AMA  X=SQR  T ( G *A MA XI 
OELS-OMI NOSORT (DMI Nt/AMAX 
OELS-SQRT(DM(N*ABS( ( QMI N/AMAX ) **2- 1 .0) ) 

RETURN 

END 
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SUBROUTINE  FIELD! I I 
IMPLICIT  REAL *8 1 A-H.O-Z) 

REAL *4  ZD.XNI1 000) tYNClOOO) .PN( 1000) *RNC 1000 )• ONI  1000 )• 
•PSNI 1000) .ZE.Xf  1 000) .Y( 1000 I.PI1000) «RHQ( 1000)  • 

• TNI 1000) .PSI(IOOO) .Q! 1 000) . THET A! 1000 i 
COMMON  /NVALU/ZD.XN.YN.PN.RN.QN.TN.PSN 
COMMON  /VALUE/ZE  *X  «Y  .P.RHQ.Q.THET  A.PSI * N 
COMMON  /DGAS/DZ  »G *R0 «GM1 • GM1H. GM1G . GGM l • GP 1 . GP GM 
C OMMON/D STAG /PTODI 1000) .TTOD,  ATOD 
COMMON/IDE8UG/IDBUGC80 ) 

COMMON  /P0INT/X1 *Y1 .PI «R1 .01 .T 1 « THl . PS l • BE1 * X2 * Y2. P2. 
•R2.Q2.T2.TH2.PS2.BE2 . P20.T20.PS0.DEL  (4  ) .PI  ( 4 I.THI  (4  ). 
*P1H  4)  .BE  1 1 4).  XI  ! 4)  .YU  4)  .DEI!4).DX  ( 4 ) • DY  ( 4 > «DL  ( 4 1 . 

*PP|  5.4)  * T 1 2 .PI  2 *B  1 2 «T2  1(4)  *P2  M4)*B2I14}«  D2I  (4  ) 

COMMON  /FLDBDY/SI NT1 2 . COST1 2 . S I NPl 2. COS PI 2.S IN B12. 
•COSB12.SIND2I  14)  .COSD2  I <4  ) . S INT2  1 !4  ) . COST  2 1 1 4 ) , 

*S INP2I (4) .C0SP2I <4)  . SIN82 1(4) .COSB21 14) .OXONI4 1 * 

*DYDN<  4)  .DZDN(4)  .0TDX(4)  . OTDY  14  ) .DTDZI4  ) .DPDXC4  ). 
*DP0Y(4)  .DPDZf  4)  .OERC  2 .5  ) * AL( 3 .4  ) .CR 14 ) • AMI 3. 3)  .CM 1 3 • 4 ) 
l FORMAT! * 0 FAILED  IN  FIELD  AT  POI NT • . 15 .5X. • Zl*» . 

4 1PE 12.5.5X. • DZ= * .El  2 .S//5  X . * XI =* .El  2 . 5 .SX* • Yl» « .El  2.5. 

• 5X.  aPl=>  .E12.5.5X.  *T1^*  .El  2.5  «5X.«PSls*  .El  2.5/5X.  • X2>* 

• *E1 2*5 *5X . • Y2  — • .El  2.5  «SX i 'P2*1 .El 2.5.5X. • T2=* . E12.5. 5X 
*.  *PS2—  * .E12.5V45X.  • P2  0=*  .El  2 . 5 ,4X  • *T  2 0=  • . E 1 2 .5,  4X . 
**PS20= • .E12.5//I 

X1=X(  I ) 

Vi=Y(I) 

DELI  11=0. DO 

IF( XI .NE.O.DO.OR.V1. NE.O.DO)  DELI  1 ) =DAT AN2 { Y 1 • XI ) 

IF (DELI  1 ) .LT.O.DO)  DELI  1 ) =6 .2831 8530 TOO +DEL! 1 ) 

DEtm^OELUl 

DO  20  I T=2 .4 

DEL!  IT)=DEL! IT-1  I ♦ 1 • 570796327D0 
IF (DEL! IT). GE.6.2831 85307D0) 

•DEL! IT)=DEL( IT)-8. 28318530700 
DEI!  IT»=OEL!  IT) 

20  CONTINUE 
P1=P!II 
R 1 =RHO  ( 1 ) 

Q 1=011) 

TH1=THETA(I > 

PSl=PSIII) 

T 1 =P  1 /{  RO  *R  1 ) 

P2=P  1 

R2=Rl 

Q2=Q1 

TH2=TMl 

PS2=PS1 

T2  = Tl 

BE1=DARSINID50RTIG*P1/R1)/0I) 
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BE2=BE1 
00  21  1T=1.4 
TH  1 1 1 T)=TH1 
PII(IT)=PS1 
BEK  IT  ) =BE  1 
21  CONTINUE 

DO  26  ITsl.2S 
T12=0.  500*1  TH1  + TH2) 

PI  2—0*  SDO  + I PS1 *PS2 I 
8 12=0.500*1  BE1  +BE2  I 
SI NT  1 2=0  SI NCT12I 
COST12=DCOS( T12) 

SI NP 1 2=0 SI NIP12) 

C0SP1 2=DC0S(  P12I 
SI NB l 2=DSI NIBI2I 
C0SB1 2=DC0S(B1 21 

X2=X1  + SINT1  2*DZ/I  COST124COSP12) 

V2=V  1 + SI  NP1 2 *OZ/COSP1  2 
00  22  J=1 .A 

T2II  J)  =0.50  0*1  TH2*THHJI  ) 

P21 ( J)=0.5D0*CPS2+PII ( Jl ) 

B2IC J)=0.5D0*f  BE2+BEM J)  ) 

D2I  C J ) =0.500*1  DELI J)+DE 1 1 J)  ) 

SINT2I I J) =D  5INCT2I1J) ) 

COST2XC JI=0CCSCT2IC J) ) 

SINP2I  ( J)=DSIN(P2H  Jl  » 

C0SP2I  1 J)=DC0S(P2  If  J>J 
SINB2I I J)=DSIN(B2I ( J)  I 
CO  SB 2 1 I J) =DCOS (B2IIJI) 

StND2I ( J)=OSI N(02I( Jl  ) 

C0S02I1 J)=DCOS(D2I4 Jl  ) 

F 1— S I NB2I ( Jl *C0SD2 I ( Jl 
F 2= SI NS 2 1 ( J| *SIND2I < Jl 
F3=COS92I ( J) 4C0ST2I < J)-F1 *S INT2I ( J ) 

J DL( J>=DZ/<F3*C0SP2I I Jl — F2*SJ NP2I ( Jl I 
DXI  J)  = SINT2I<  J)*C0SB2I(  J) +F1  *C0ST2  I ( J) 

OYI J )— F3*SI NP2I < J)+F2*C0S P2 I ( J ) 

XII  J)  = X2-DX(  J)  »DH  J) 

VI  ( J)  = Y2-DY<  Jl  *DL<  Jl 

CALL  GETPTIXI  CJI.YM  J)  .PPI1.J11 

I FI OAB  SIPPt  A, J) ) *LT • 1 *0—8 1 PPI4. J)=0.00 

IF( 0ABSIPP1 5 , J) ) *LT. 1.0-8)  PP(S • J)=0.00 

CALL  OELTAF I J.DELI J) «DEI( Jl ) 

IFIOABSf 0E1I J) ).LE. 1.0-8)  DEI < J 1=0. DO 
CALL  GETOERIXI  I Jl.YK  J), DERI 
Pt< J)=PP( l.JI 
THI I J) =PP ( 4 « J) 

PIK  J)=PPI5.JI 

BE  I ( J) =OARSI N{ DSORT IG+PPC 1 . JI/PPC2 • J ) 1/PPI 3. J 1 1 
T2I  { J)  =0.500*1  TH2FTHI  1 Jl  ) 
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P2I  C J > =0.50  0*1  PS2  + PI  I ( J)  1 
B 21  ( J)  =0.  500*1  BE2+8EI  I J>1 
0 21 ( Jl =0. 5D0*(DEL( J)  +DE I ( J) I 
SCNT21  ( J)  =0  SI  N<  T2UJ)  I 
COST2K  J)=DC0SIT2I<  J)  ) 

SINP2I < J) =DSI N«  P2I I Jl  > 

C0SP2I  < J)=OCQS(P21(J)) 

S1NB2I  ( J1=DSINIB2  1 1 Jl  1 
CO  S3  21  I J)=DCOS(B2  f(J)l 
SIND2M  J) =D SIN (0211  J)  1 
COSO 21  C J)  =DC05  (02 1 4 J)  I 
DTDXC J1=0ER( 1 *4} 

DPDX( J)=DER( I ,51 
DTDYI  J ) =DER<  2,4) 

OPOY( J)=DER<2,51 
Fl=COST2I( J)*SIN02I ( J1 
F2=SINT2I ( J) 4SIND2I < J) 

F3=SIN02I I J) *DL( J) 

F4=C0ST2I I J>  4C0SD21 < J> 

DXDNI  J)=-:F1 

OYONC J1=F2*SINP2I C J1+C0SP2I I Jl *COS 02  I < J 1 
DZONI Jl=F2*COSP2I C J1-SINP2I C J)*COSD2Il J) 

F5=XJ(  J1-X2 
F6=YII JI-Y2 

OTDZI J)=(DTDX( J> *F5*0TDY( Jl*F61/OZ 
OPO  Z( J ) — I OPDX ( J) *F5fOPOY( J) 4F61/DZ 
AL(  1 , J1=C  0SB2I ( J! /( SI N82I I Jl *PP ( 2 • J ) *PP4  3, J 1442) 

AL ( 2,J1=C0S02I I J)-F3*SIND2! t Jl*DZONI JI/OZ 
ALC3,  Jl=FH-F3*F4*DZDNt  Jl/OZ 

CR ( Jl=F3*( SI  NO 21 ( J 1 *IDTOX I Jl *DXON< Jl+DTDY ( J1*DY0N(J 1* 
♦DTDZI  J>  *DZDN(  J)  )-F4*|DP0X  t J1*DXDN(  JH-OPOY  ( J)*DYDN(  J >♦ 
*OPOZt  J)  *OZONl  J)  > 1+ALU  • JI*PI<  J1  + ALI2,  J 1*TH 1 1 Jl  *AL  I 3.  J 1 
**P  K 1 I J 1 

22  CONTINUE 
00  24  M=1  ,4 
L=M 

DO  23  J=1 ,3 
AM|  J,1  )=AL(  I ,L) 

AMI J.2)=ALI2 ,L) 

AMI J,3)=AL(3.L) 

CM<  J , Ml -=C R (LI 
L=L*1 

1FIL.GT.41  L=1 

23  CONTINUE 

CALL  SIMQf AM.CMC1 .Ml ,3.11 

24  CONTINUE 

P2  =0. 250041 C M| 1 .11 +CM( 1,21 +CMI 1.3) +CMI 1*411 
TH2=0« 25D  04ICMI2.il +CM< 2,21 +CMI 2 ,3 ) FCMl 2. 41 1 
PS  2=0*  2 50  0*(CM(3,1)  *-CM{  3,21  +CMI  3,31  4- CM  (3*411 
T2=TT00*C  P2/P  TOO I { 1)**GM1G 
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R2=P2/CR0*T2> 

02=0  SQRTC 2*00 *GGM1 *R0+<  TTQ0—T2 J » 

IF ( IT* EQ* 1)  GG  TO  25 

IF ( DABS! ( P2  — P2 0} /P2 ) *LE* I • 0—71  GO  TO  27 
IF( IT.EQ.25I  GO  TO  26 

25  CONTINUE 
P20=P2 
R20=R2 
020=02 
T20=TH2 
PS0=PS2 

8E2=DARSIN(0SQRT(G*P2/R2)/Q2) 

26  CONTINUE 

IF ( DABS!  i P2—P20) /P2  ) *t_E*  t *0— 5>  GO  TO  27 
WRI TEC 6*1)1 .ZE.DZ.X1 ,V1 .PI *TH1 • PS1 *X2* Y2. P2# TH2.PS2. 
♦P20.T20.PS0 

IF { DAB  SC ( P2— P20 )/P2 ) *GT* 1 • 0— *|  STOP 

27  CONTINUE 
XNCII=X2 
VNC  I l=Y2 
PNC  I )=P2 
RNC  I J=R2 
QN( I ) — 02 
TNC I )=TH2 
PSNC I >=PS2 
RETURN 
END 
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SUBROUTINE  FINDCI  *L1*L3«NEI 
INTEGER *2  NEC  1 1 
REAL*4  DISTC3) 

COMMON/VALUE/ZE .X  C 1000 1 . Y ( 1000 ) 

L2=L1+L3 

J=0 

OO  22  L-Ll *L2 

0=1  XI  I 1-XtLl  1 **2«-<YCt  > -YCLI » *42 

IFCJ.EQ.3)  GO  TO  21 

J=  J4-1 

NEC  J) =L 

DISTC  J1=0 

KF( J.NE.3)  GO  TO  22 
CALL  S0RT1OIST.NE.3) 

GO  TO  22 

21  CONTINUE 

• IF*  O.GE.DISTO)  ) GO  TO  22 
OISTC  31=0 
. NE(3)=L 

CALL  SORT CO I ST  #NE  *3) 

22  CONTINUE 
RETURN 
END 
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SUBROUTINE  FINDBPf XI *Y1 »XB.YB . {POINT ) 

COMMON  /6EOM/ZCI50I . XC ( 50 .50 ) • YC (SO • 50 ) *NS «NP 
DIMENSION  XB(50) «YB(50) 

20  CONTINUE 
IPOINT=l  . 

DMI N=( Xl-XBC I > ) **2  + C Y1 -VB(I ) ) **2 
DO  21  J—2 , NP 

D=»(  Xl-XBC  J)  ) **2+CYl-YBC  J>  > **2 
IFC  O *GE  «D  MI N ) GO  TO  21 
IPOINT-J 
D MI  N=D 

21  CONTINUE 
RETURN 

ENTRY  CKBDYPIX1 »Y1 .XB.YBt IPOINT) 

J-IPOI NT 

DMI N=C Xl-XBC J)) **2  + C Y1 + Y8C  J) 1 **2 
J=J-1 

IF(J.EQ.O)  J=NP 

0=(  Xl-XBC  JM  **2*{Y1-YB{  J>  ) **2 

IF(D.LT.DMIN)  GO  TO  20 

J— I POI NT* 1 , 

IFC  J»GT«NP)  J=1 

D=(  Xl-XBC  J)  t **2*(Y1-VB|  J>  1**2 

IF(D.LT.DMIN)  GO  TO  20 

RETURN 

END 
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SUBROUTINE  FI  TCI) 

COMMON/VALUE/ZE  .X(  1000  ) *YI  1 000  > tPtlOOO  ) *RHO<  1000)  • 
*Q< 1000) .THETAC 1000) .PSl (1000) .N 
INTEGER42  N£(8.1000) 

COMMON  /NEAR  /NE 

REAL *8  R .8(12.12) .C( 1 2 .5 ) , V AL ( 5 1 . DER < 2 . 5 ) . OX. OY, FAC 

OO  22  J=1  .12 

1F( J.GT.9)  GO  TO  21 

8(  1 . J)=1.00 

K=1 

ZF(J.NE.I)  K*NE( J— 1 * I ) 

8( 2 . J ) = X( K) 

8( 3. J1=Y(K> 

GO  TO  22 

21  CONTINUE 

0( l • J) =0.D0 
. 8( 2. J) =0.00 
B( 3. J) =0.00 

22  CONTINUE 

OO  24  Jjs4.ll 
13—  J J 
KL=  J— 1 
Kl=l 

IF(J.NE.l)  KlsNE(KL.I) 

DO  23  K=l  »K1_ 

K 2=1 

IF(K.NE.l)  K2=NE( K— 1 • 1 ) 

R=IXCK11-XC K21 )**2M Y(K1 l-Y(K2> ) 4*2 
8 ( J J • K ) =R  *0  LOG  I R ) 

23  CONTINUE 
8(JJ.J)=0.D0 

24  CONTINUE 
8(12*1 1=0. DO 
DO  26  J=1  .11 
K=1 3- J 
KL=K-1 

DO  25  L=1 .KL 
0( K . 1 3-L)=8( L • J) 

25  CONTINUE 

26  CONTINUE 

DO  29  J-1,12 
IFfJ.GT.3)  GO  TO  28 
DO  27  K=l * 5 
C( J.KlsO.OO 

27  CONTINUE 
GO  TO  29 

28  CONTINUE 
KL=1 3- J 
K=I 

IF(KL.NE.l)  K=NE ( KL-T  • 1 ) 
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CC  J.I |=pf KJ 
CC J.2l=RHO( K) 

C<  J.3) -Q( K) 

C< J,4)=THETA«K) 

C( J*5)=PSI( K) 

29  CONTINUE 

CALL  SX  MQf  B.C.12.51 
RETURN 

ENTRY  GETPTf XI  »Yl  «VAL) 

DO  3 0 J*1  .5 

VAL.C  Jt^Cfl2»J)+CflliJ)»Xl  + CC10»JI*Yl 

30  CONTI  NUE 

DO  32  J-l *9 
K=l 

IFfJ.NE.lt  K-NE(J-l.l) 

R=<  XI- X( K) ) **2+< Yl-Yt K1 > **2 
IFCR .EQ. 0.00)  GO  TO  32 
DO  31  L— 1 .5 

VALC  L)  *VALCUtC(  J .LI  *R*DLOG(  Rl 

31  CONTINUE 

32  CONTINUE 
RETURN 

ENTRY  GETDERCX1  .X2.DER) 

00  33  J=1 .5 
DERf I.Jt— Cf 11 • J) 

OER ( 2.J)  =C(  1 0 . J) 

33  CONTINUE 

DO  35  J-l  .9 
K=I 

IF(J.NE.l)  K=NEtJ-l.I) 

DX=X1-X(K> 

0Y=*Y1-VIK> 

R=DX**2fOY*+2 
IFfR.EQ.O.DOt  GO  TO  35 
DO  34  L=1 .5 

FAC=2.D0*C( J *L) 4( 1 .DO+DLOGf R| ) 

DERf 1.L1=DER<1 .L1+FAC40X 
DER(2.L)=DER(2.L)+FAC*DY 

34  CONTINUE 

35  CONTINUE 
RETURN 
END 
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SUBROUTINE  I DENT  (II 
IFd.EQ.2l  <50  TO  21 
CALL  CALCMP415. 0.2. 0.10.0) 

CALL  SYMBOL  40. 5.7.5 .0.  84  . * FOLLOW  I NG  PLOTS  FOR«.0.0.19> 
CALL  SYMBOL! 1.2  .5. 5 .0.98 . *W.  C.  ARMSTRONG*  .0.0. 15) 

CALL  SYMBOL!  5.5.3.5.0.B4.*  CSB— EAO*  .0.0.7) 

CALL  CLASS!  1 .3.0.0) 

CALL  CALCMP!  O.S.O.S.O  .3) 

RETURN 
21  CONTINUE 

CALL  CALCMP! 0.0. 0.0. 0.3} 

CALL  CLASS! 4.3.0.0) 

CALL  SYMBOL ! .00 *7. 5.0. 84.*  PRECEED ING  PLOTS  FOR*.0..20) 
CALL  SYMBOL! 1.2.5. 5.0.98. •».  C.  ARMSTRONG* . 0 .0 . 15) 

CALL  SYMBOL! 5 . 5 .3 . 5 .0. 84 . * CSB— EAO* .0 .0.7) 

CALL  CALCMP! 0.0 .0.0. 9999 .21 

RETURN 

END 
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SUBROUTINE  INGEOM 
COMMON  /TI TLE/HEADI20I 

COMMON  /GEOM/ZC ISO ) .XC(50.50> ,YC(50  .50J.NS.NP 
COMMON/I  DEBUG/I QSUG! 80) 

COMMON/VALUE/X(50.5f»  .YI50.50  ) *5  (50  ) • IT  (SO  ) • IV  ( SO  ) 

1 FORMAT! *1 *.10X.20A4//20X. ‘GEOMETRY  INPUT*. 10X* 

♦ * NO.  OF  STATIONS  ='.I3/>  j 

2 FORMAT! ' 0 Z - * ,1 PEI  3 • 5 *5X . • I V =*,13.5X,*IT  =».I2// 

♦ 8X, « X*  *1 4 X . *Y  •/> 

3 FORMAT! 1P2E1 5.5) 

4 FORMAT!  * 1 * . 1 OX .20A4//20X. * CONVERTED  GEOMETRY*//) 

5 FORMAT!2!S«EiO.O) 

6 FORMAT! 8E 10 • 0 ) 

C IV=  NO.  OF  POINTS  AT  THIS  STATION 

C I T=  1 - ALL  POINTS  GIVEN  IN  RECTANGULAR  COORDINATES 

C 2 - ALL  POINTS  GIVEN  IN  POLAR  COORDINATES 

C 3 - ONE  POINT  GIVEN  IN  POLAR  COORDINATES 

READ ( 8 .5 ) NS 
IF(  I DRUG!  I I • NE . 0 ) 

♦ WRITE! 6.1 ) HEAD  * NS 
00  21  1=1 ,NS 

RE  AO !8.5)IV!t) .IT!I).ZC!I) 

IF! IDBUG!  D.NE.O) 

♦ WRITE!  6.2) ZC! I). IV! I). IT!!) 

N=I V! I ) 

IF I IT! I ) «EQ. 3)  N=l 

READ! 8,6 > I X! J. 1 1 , Y! J.I ) • J=l ,N) 

IF!  IDBUG!  D.NE.O) 

♦WRITE! 6.3) t X! J.I 1 «Y! J* I ) . J=1 «N) 

21  CONTINUE 

DO  26  1=1  .NS 

IF! IDBUG!  2 ) • NE . 0 ) 

♦WRITE! 6.4) HEAD 
IF ( IOB  UG C 2 ) • NE . 0 ) 

♦ WR I TE  (6*2)  ZC  ( I )*tvm.IT(() 

N=I V! I ) 

IF(  IT!  n.EQ.l)  GO  TO  23 

IFCITCI) .EQ.2)  GO  TO  23 

0T=360.0/N 

Y(  1 .1  )=0.0 

DO  22  J=2,N 

XI  J.I  )=X(  1,1) 

Y!  J .1  )=!  J-t  ) WOT 

22  CONTINUE 

23  CONTINUE 

DO  24  J=1 *N 

XX=XI  J.I  )4CQS!0.  017453294YIJ.il) 

Y(  J • I ) =X(  J.  t ) 4SI N!  0.  01  7453294Y ( J.  1 ) ) 

X(J.I)=XX 

24  CONTINUE 
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25  CONTINUE 

IF  f I DBUG(  2) *N£«0) 

*WRlTE(6i3)( X ( J • I) ,YI J.I I ,J=1 ,N» 

26  CONTINUE 
NP»IVlll 

DO  27  I =2 »NS 

IF( NP»GE  *1  VC  I > 1 GO  TO  27 
NPalVI I > 

27  CONTINUE 
NPM=NP-1 
FAC*1 ,0/NRM 
DO  33  lsl,NS 
N=I VI I) 

Sf  1)=0.0 
DO  28  J=2,N 

SI J)=SI J-l) +SQRTI < X( J, I ) —X ( J— 1 * I ) ) **2* 
*1  VI  J.tl-YI  J-l. I)  1**2) 

28  CONT I NIC 
OS=S!N»*FAC 
XCt 1,1 )=X(1,I) 

VCI  1,1 I=Y<1  ,I> 

XC(NP,I)=X(N,I  > i 

YC<  NP,I)=Y(N*I  } 

KL=2 

on  32  J=2 ,NPM 
3P=OS*<  J-l  ) 

DO  29  K=KL,N 
IF(SP-S<Kn31  ,30,29 

29  CONTINUE 
K=N 

30  CONTINUE 
xcu,n=x<K,i) 

YC  I J * I )=>Y(K,I  ) 

KL=K 

GO  TO  32 

31  CONTINUE 

RAT— I SP— SI K— 1 I |/<S<K>-S{K-1»I 

XCI J*I 1-XIK-I  ,1 )♦  I X(K«1 1— XIK— 1 ,11  I ♦RAT 

YCI J,I1=Y|K-1,II*<VIK,II-Y(K-1 ,1>>*RAT 

KU=K 

32  CONTINUE 

33  CONTINUE 
RETURN 
END 
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SUBROUTINE  INPUT 
COMMON  /TI TLE/HEAOf 20) 

COMMON/VALUE /ZE  *X(1000) *YC 10001 *P< 1000  I .RHO* 1000)  • 

♦Q<  1 000)  * THETA  < I 000)  *PSI  UOOO)  »N 
COMMON  /STAG/PTGt 1000) .TTO. ATOT* ZO . 2MAX .NB( 200 ) * 
♦NBDC200)  • JMAX,  JPT  .1  B.  I PLOT  , IT  1 1 • IT  12 
COMMON  /GAS/G*RO,GM1  ,GM1 H*  GM1 G, GGM1  , GP1  • GPGM 
COMMON/* DEBUG/I  DOUG (80 ) 

INTEGERS  NB.NHO 

COMMON/NVALU/XR  *20) «YR*20) ,PR<20 ) *THR ( 20 ) • PSR* 20). 
*QR(20) »PTR*20) 

1 FORMAT* M •«10X«20A4//20X, 'STARTING  PLANE  lNPUT*//5X, 

♦ •GAMMA  = • ,F5,3  ,5X,*PT0  - * *1  PEI  2 .5 ,5X , • TTO  =»,E12.5, 
*5X.*R  =*,E12.5> 

2 FORMAT  C * 0 ZO= • *1 PEI  2 ,5 ,5X , • ZMAX=* , El  2 , 5 , 5X, • IT YPE= • * 

♦ I2.5X,  «IMV=*  ,I2,SX,*NPTS=».IS,SX,*NRAYS=*  • I3.5X* 

♦ • JMA  X=  **I5,5X, • JPT=*  * 1 5/A AX *• IPLOT=» , 12, 4X* • IT  11=  *.  12. 
*5X, *IT12=« ,12 ,7X,'INEIGH=* . I 2//8X , • X» *14X, »Y* , 1*X. • P* , 

♦ 14X.  *Q  1 *1  2X,  • THETA*  ,1 1 X . • PSI  * « 1 2X , • PT 0* // ) 

3 FORMAT* 1 P7E1 5,5) 

♦ FORMAT*  7E 10,0,15) 

5 FORMAT*  2E  10,0*1015) 

6 FORMAT* »1 • *10X,20AA//20X. 'STARTING  PLANE  VALUES*//8X. 

♦ •X'flAX*  • V • *I4X,  'Pf*l3X,  *RHO*  *13X*«Q*  , i 3X  • • THT  • « 1 2X, 
♦•PSI ■ ,X2X* •PTO»//*XPaElS.5l ) 

7 FORMAT* «0  TOO  MANY  BODY  POINTS*) 

C I M V=0  0*1)  =Q  THETA*  l )=THET$  PS  I*  I )=PS  I 

C IMV=1  Q*I)=M  THETA* I )=THETA  PSI<I)=PSI 

C IMV=2  Q ( I ) =W  THET  A*  I )=U  PS  I < I)=V 

C ITYPE=0  ONE  RAY  GIVEN  IN  POLAR  COORDINATES 

C ITYPE=1  PLANE  GIVEN  IN  POLAR  COORDINATES 

C ITYPE=2  PLANE  GIVEN  IN  RECTANGULAR  COORDINATES 

READ*  8,4)G*PT0 ,TTC,RO 
WRITE* fi.l)HEAD,G,PT0,TT0,90 

READ*  8 ,5  ) ZO.ZMAX,!  TYPE  . I MY  .NPTS  . NR  AYS  , JMAX,  JPT  * IPLOT. 

♦ IT11  ,IT12  * I NEIGH 

WRITE  (6,2)  ZO  *7  MAX,  I TYPE,  I. MV, NPTS,  NR  AYS*  JMAX,  JPT,  IPLOT, 

♦ IT1  1 * I T 1 2 ,1  NE  IGH 
N=NPTS 

I B=Q 

DO  20  1=1 ,N 

READ*  H, 4) X(  I ) , Y<  r ) ,P(1  > ,0(1 ) , THET A* I ) , PS  I * ! 1 ,PTO*  * ) ,NBP 
IF ( PTDCI ) ,EQ, 0, 0)  PTO(I)=PT0 
I F ( IMV.FQ.2)  GO  TO  IP 
THETA* I ) =0,01 7453294THETACI ) 

PSI  { I )=0.  01  7A5329+PSI  *1) 

19  CONTINUE 

TF(  NEJP.FO.  0)  GC  TO  20 
T3=IO+ 1 

IF*  IF,  GO  TO  99 
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NS{ TB)=I 

20  CONTINUE! 

IF<  IDOUG!  3I.NE.0) 

*w^iTrc6.3)(x(n.Y(i).p(n  ,o(  i > »thetmi  j «psic  i>, 
*PTO(  I > .1=1  *N) 

1FC ITYPE.EQ.l)  GO  TO  24 
IF! ITYPE.E0*2»  GO  TO  26 
DO  21  1*1  *N 
XR(  1 > = X<  I ) 

YR« I l=Y( 1 ) 

PR(I)=P!I> 

QR  ( I )=Q!  I ) 

THftf I >=THETA< II 
PSR!  I l=PSI<  1 1 
PTR!I  »*PTO!IJ 

21  CONTINUE 

ie*i 

J=0 

N8I=M8< 1 ) 

DO  23  1=1 .N 
J=J41 

X( J) =XR t I ) 

Yf  J)*YR(I > 

P( J)=PR(I ) 

Q<  J)=GR( t ) 

THETA | J)»THR!IJ 
PSH  J)=PSR!I  K 
PTOC  JI=PTR(I  )' 

IF(I.EQ.l)  GO  TO  23 

NPTS=3 • 1 41 593/ARSI Nl 0.5/FLOAT ! 1-1 1 I+0.5 

DT=360.,0/NPTS 

NPT  S=NPTS— 1 

IF(NBI.EQ.I)  NB( I8)=J 

DO  22  K=1 .NPTS 

J=J+1 

X(J)=XRII) 

Y(  J)*YRtl  )*FLOAT< K)*DT 
P( J)=PR(I I 
Q( J)=QR< I ) 

THETA!  J)*THR!  1 I 
PS1  ( JI=PSR«I  I 
PT01 JI*PTR!I» 

IF(NBI.NE.I)  GO  TO  22 
IB=IB+1 

IF! IB.GT.200)  GO  TO  99 
NB!1B)*J 

22  CONTINUE 

23  CONTINUE 
N=J 

24  CONTINUE 
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00  25  I=i .N 

ST*SI Nl  0. 01 745329 *Y f I ) ) 

CT=COS<  0.  01 745329*Yl  I ) I 
YC  II=XC1  I *ST 
X(  I ) = X ( I I *CT 

TANP=SQRT ITANI  THE  TA< I ) ) **2+T ANIPS I C I ) 1**21 
PSI C I I =ATA  Nl  TANPYSTI 
THETA!  t ) =ATAN( TANP*CT*COS<PSl ( I ) ) ) 

25  CONTINUE 

26  CONTINUE 
ATOT=SQRT(G  *RO*TTGJ 
GM 1=G-  \ . 0 
GM1G-GM1  /G 
GP1*G+1.0 

GMIH— 0«5*GM1 

GPGM=GPI/GM1 

GGM1=G/GM1 

XNB=0.0 

YN9  = 0.  0 

DO  261  1*1  .IB 

J*NS ( t ) 

XN8=XN8*X(  J1 
VNB=YN8+YlJ) 

261  CONTINUE 

XN8  = XN3/IB 
YNB*YNB/I  B 
DO  27  1=1, IB 
J=N0 1 I ) 

RHOd  )=ATAN2(  Y(  J)-YNB.X(  Jl-XNB) 

IF |RHO( I I.LT.0.0)  RHOI I )=6 .2831 85+RHGl I } 

NBOII >=J 

27  CONTINUE 

CALL  SORT! RHG.NBO.IB) 

NB( I B* 1 ) =0 

NBO< IB*1 >*0 

DO  30  1=1. N 

IF ( IMV.EQ. 1 ) GO  TO  28 

V=Q< I > 

IF  C IMV.EQ. 2)  V=SQRT(V**2+THETA|I)**2+PSI(II**2) 

SM={ V/ATOT1 **2 

SM=SQRTC SM/I 1 . 0— GMl  H*S  N) 1 

IF < IMV.EQ. 0)  GO  TO  29 

PSI  II  )=ATAN2CPSI|I>.Q(in 

THETAC  I )=ATAN2(  THETA  It  » .01 1 >/COS  I PS  I < I 111 
Q ( 1 l=V 
GO  TO  29 

28  CONTINUE 
SM=QII) 

Q ( 1 ) = SM*A  TOT /SORT  1 1 .0  + GM1  H*SM**2  I 

29  CONTINUE 
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SMS=SM**2 

SMFAC=1.0+GM1H*SMS 
T=TTO/SMFAC 
RHQC I )=P( I)/(RO*TI 
30  CONTINUE 

IF(  IDBUGI  4 ) • NE • 0 ) 

*URI  TE(6»6)HEAD«(X(I  ) ,Y C I > *P< I) *RHO< I I • Q < I I • T HET A(  I ). 
*PSK  I ) *PTO<  I ) .1=1  *N1 

IFll  ITYPE.EQ.O)  .AND.  (INEIGH.EQ.OM  INEIGH=1 
IF( INEIGH.EQ.O)  CALL  NEIGH 
IF(  INEIGH.EQ.M  CALL  RNEIG 
IF( 1NE IGH.EQ.2)  CALL  READNE 
RETURN 
09  CONTINUE 
VR1TE<«.7> 

STOP 

END 
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SUBROUTINE  NEIGH 

COMMON/VALUE/ZE.XI1000I.YC1000I  *PC  1000  I *RHDC  1000)  • 
*0(1000) , THETAI 1 000 ) .PS1 (lOOO).N 

integer *2  nei&.iooo) 

COMMON  /NEAR/NE ,DMIN 
COMMON/] OE8UG /I D8UGC80 ) 

DIMENSION  LIST!. 19)  ,L0NG(19) 

DATA  L I ST/1 *2*  8 *20 • 39*64*95  *133 *177* 227 *28 A *34 7*4 16* 
*491 .573,661 .755.856.963/ 

DATA  LONG/0.5 .11 • 18 .24 ,30 .37 .43 .49 .56. 62. 68* 74 .61 .87, 
*93.100,1  06,0/ 

1 FORMAT I • 1 ■ *2 OX, * NEIGHBORS  1 »5X « * DM I N =«,1PE12.4 
*//6X , ‘BASE • ,9X,«1 * *9X« • 2 • ,9X*'3« ,9X,"4« , 
*9X.,5,,9X,*6,.9X.*7,.9X.,8*//) 

2 FORMAT!  91 1 0) 

DIMENSION  DXM(B) 

DM  I N=l  «E  1 0 

DO  23  is]  ,N 
IC-0 

DO  22  J-  l.N 

IFI  I.EQ. J)  GO  TO  22 

D=(X(I  J-Xt  J)l  **2  *■(  Y(  I ) —V  ( J)  ) **2 

IF(  IC.E0.8)  GO  TO  21 

IC-IC+1 

NE  ( I C , I ) 3 J 

DXMI  IC)=D 

IFIIC.EQ.8)  CALL  SORT t DXM ,NE ( 1 • I ) ,0 ) 

GO  TO  22 

21  CONTINUE 

IF(D.GE.DXM(8)  ) GO  TO  22 
NE<  8,1  ) =J 
DXM(  8)=D 

CALL  SORT (OXM,NE( 1 ,1 ) «B ) 

22  CONTINUE 

IFIDXMI 1 I .LT.OMIN)  DMIN=DXM(1) 

23  CONTINUE 
GO  TO  29 
ENTRY  RNE1G 
J-l 

DO  28  1=1, N 

IF(  I .GE  .LIST!  J+l)  ) J=J+1 

IFI J.NE. 1 ) GO  TO  25 
DO  24  K=l  .6 
NEIK.I )=K4l 

24  CONTINUE 
NE( 7.1 )=11 
NEC  8.1 1=17 
GO  TO  28 

25  CONTINUE 
LL=I — 1 
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IFCLL.LT. LISTC  JM  LL— LI  ST<  J*1  ) — 1 
LH^IM 

IFCLH.GE. LISTC  J+l))  LH=LISTCJJ 

NEC  1 • t l=LL 

NE  C 2 • I I =LH 

IFCJ.NE.2)  GO  TO  26 

LL=LL-1 

IFfLL.LT. LISTC JJ ) LL»LI ST C J*1 I -I 
LM=LH+ 1 

IFCLH.GE .LISTC J+i >)  LM=LIST(J) 

NEC  3.1 1 -LL 
NEC  4*1  ) 3*1 
NEC  5.1 I=LH 
K=J+  1 
GO  TO  27 

26  CONTINUE 

CALL  FINDC I .LISTC J-l) .LONGC J-l I . NEC 3. 1 » I 
K— J-fr  1 

IFC  N.LT.LISTCKU  K=J-2 

27  CONTINUE 

CALL  FINDC I .LISTC  K) .LGNGC  K) • NEC 6 . 1 ) ) 

28  CONTINUE 
GO  TO  29 
ENTRY  READNE 

READC  8 .21  CK.CNEC  J.KI  . J=*l  .8)  .I=l.N> 

29  CONTINUE 

IFI  COBUGC  5)  . NE  .01 
♦ MR I TEC  6.1  IOMIN 
IFCIOBUGC  5I.NE.0) 

A MR X TEC  6.21  C I *C  NEC  J.II.J=1.8) *1=1. N> 
return 
END 
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SUBROUTI NE  NEWRAPI  T«P*A»8»C*N1*N2» N3) 

IMPLICIT  REAL*8(A-H,0-ZI 
R£AL*8  Nl*N2*N3*J 

1 FORMAT  < ^ 0 FAILED  TO  CONVERGE  IN  NEHRAP  T*«.1PE12.5. 
*5X»  1 TN=  ■ .E12.5.6X*  1 P=*  ,E1 2*5 *5X PN= • , El 2.S/5X , «TO»*  • 
+E12.5.5X,  *PO*«  .E12.5/6X  ••  A- • ,E1 2.5 ,6X. • B=  • . E12  .5.  6 X. 
*'C*» .E12.S.5X. *N1 «E1 2.5 »SX * »N2»* •E12«5*  5X. * N3** • 
♦E12.5J 
TO— T 
PO=P 

DO  22  1=1*50 
SINT=OSIN«T» 

COST=DCOS< T> 

SINP=OSIN<Pt 

COSPsDCOSIPl 

F-A*T+B*P-C 

G=Nl*COST*COSP+N2*SI  NT+N3*COST*5  INP 

DGT=-N1*S1  NTACOSP4-N2  *COST -N3*S1NT*SINP 

DGP=— N1 «COST*SI NP+N3*C0ST*C0SP 

J=A*DGP-B*DGT 

TN=T-<F*DGP-G*B1/J 

PN=P+(F4DGT-G*A)/J 

IF(DAQS((  T-TN)  I.GT.L.D-61  GO  TO  21 
IF<DABSC<P-PN»  ).LE. 1.0-6)  GO  TO  23 

21  CONTINUE 
IFCI.EQ.50)  GO  TO  22 
T=0.5D0*I  T+TNJ 
P-0.5D0*{  P+PN) 

22  CONTINUE 

MR  I TEC6*l}T,TN*P*PN*Ta*P0* A.B.C.Nl *N2*N3 

CALL  ERRWCA 

STOP 

23  CONTINUE 
T*TN 
P=PN 
RETURN 
END 
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SUBROUTINE  NORMAL CX »Y * Z* A , 8. C .N1 .N2.N3. IT • IB) 
IMPLICIT  REAL*8(A— H*0— 2) 

REAL  *4  ZC  • XC  • YC 

COMMON  /GEOM/ZCC  50} .XCI50.50I * YC ( 50 • 50 1 . NS . NP 

REAL48  N1 .N2.N3 

COMMON/ IDE BUG /I 0BU6( 80 I 

REALMS  F(3»5)  »FA(5)  »FB(5)  *FCI5I 

1 FORMAT! *0  NORMAL  X** • 1 PEI 2.5.SX. • Y=» ,E1 2. 5.SX. » Z=» • 
*E12.5.5X«2HA=.E12.5 ,5X. »B=* .El 2.5.5X. • C=» .El  2.5) 

2 FORMAT!  8X • ' T1  • .1 2X • • T2  ■ .1 2X • * A2 • • 1 2X  • ■ B2 • .12X**C2»  • 
*12X**N1  * .1 2 X.  1 N2  * • 12X.*N3*/IP8E14.5) 

IF ( IOBUG(  11  1 • NE  • 0 ) 

+URITEI 6.1 ) X.Y.Z.A.B.C 
NSM=NS— 1 
00  21  I *2 iNSM 
IF (Z.LE.ZCt 1)1  GO  TO  22 

21  CONTI  NIX* 

I -NSM 

22  CONTINUE 

CALL  BDVFI T{ A1  .B1.C1  *XC< 1 • t-1 ) • YC€ 1 • 1-1 ) . I B • IT  1 ) 
CALL  BDYFITf  A2,B2*C2.XC(1  .1  >,YC<l.I  1.I8.IT2) 
CALL  B0YFIT(A3.B3.C3*XC(l . I +1 ) . YC < 1 . I +1 ) . IB. IT  3) 

IF ( I T2.EQ. 4 ) GC  TO  23 

IF( ( I Tl.LE. IT2) .AND* (I T3.LE.IT2) ) GO  TO  24 
GO  TO  27 

23  CONTINUE 

IF( C I T1 .NE. 1T2) .OR. (IT3.NE. IT21J  GO  TO  27 

24  CONTINUE 

FAC  1=1 «D0/(  ZC(  I l-ZCCI-1)) 

FAC 2=1 *O0/( ZC ( I+I )— ZC( I— 1 ) ) . 

FAC  3=1  .0  0/1  ZC(  14-1  ) — ZC(  I )) 

FAC4=ZC( 1-1 ) +ZC ( I ) 

A2=CA2-A1)*FACI 
B2=( B2—S1 ) *F AC 1 
C2=( C2-CI 1 4FAC1 
A3— ( A3— A1 ) *FAC2 
B3=< 03— Q 1 ) 4FAC2 
C3=CC3-CI )*FAC2 
A3=CA3-A2I*FAC3 
B 3=( B 3— B2 ) *F  AC3 
C3=(C3-C2»4FAC3 
A 2=  A 2- A 3 *F  AC  4 
B2=B2— B34FAC4 
C2=C2— C3*FAC4 

Al=At-(A2*A3*ZC< 1-1 } )*ZCC I-ll 
B1=B1-(B2*B3*ZC(I-1I 142CII-1) 

C1=C1-CC2+C3*ZC[I-1> l*ZC( 1-11 
AX=AIF( A2+A34Z) *Z 
8V=BU(02fB3*Z)  42 
IF ( IT2.NE.4)  GO  TO  25 
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N2=-2.D0*t X-Al-f A2*A3*Z)*Z> 

N2=— DABSt  N2 } *0 S I GNt 1 • DO  .X-AX ) 

N3=-2.00*( Y-Bl-tB2*B3*Z> »Z> 

N3-— OABSC  N3 1 *DSIGN t 1 *00  »Y— BY ) 

Nl=-N24t  A2+2.D0*A3*ZI-  N3*l  B2+2 • D0*B3*Z1 +C2»2 .D0*C3*Z 
GO  TO  26 

25  CONTINUE 

N2=-DABS<  AX) 4DS I GNI I .00 .X-AX ) 

N3=— OABSC  BY) 4DS1GNI 1 .DO .Y-BY) 

Nl=-l  A2»2.D0*A3*Z)4X-(B2+2.D0*B3*Z)*Y+CC2+2.D0*C3*Z) 

26  CONTINUE 

D=DSQRTt  Nl**24-N2**24N3*42) 

N1=NI/D 
N2=N2/D 
N3— N3/D 

IF C IDBUGt  m • NE • 0) 

♦■RITEt  6.21 T1 .T2.A2 .82.C2.N1 . N2.N3 
RETURN 
2T  CONTINUE 

IF C 1T1 .EQ.4)  GO  TO  29 

Ft l . I ) =0. DO 

Ft  1 »2)=A1 

Ft  I *3 ) =0 .DO 

Ft  I * 4 ) =B I 

FtI.5)=Cl 

GO  TO  30 

29  CONTINUE 

Ft  1.11=1.00 

F< 1 ,2>=-2.00*Al 

Ft 1.3)=1. DO 

Ft  1 .4)=— 2.D04B1 

Ft  1 • 5)  =C1-A  1442-SI  442 

30  CONTINUE 

IF t f T2.EQ.4)  GO  TO  31 

Ft  2.1 )=0.D0 

Ft  2.2) =A 2 

Ft  2.31=0. DO 

Ft2*4)=B2 

Ft  2 . 5)  =C2 

GO  TO  32 

31  CONTINUE 

Ft  2.11=1.00 

Ft2,2)=-2.D0*A2 

F(2.3l=l.D0 

Ft  2.4)=-2.D0*B2 

Ft  2.5)=C2-A2*42-B2**2 

32  CONTINUE 

IFf I T3.EQ.4)  GO  TO  33 
Ft  3.1 )=0.D0 
Ft  3.21 =A 3 


59 


AEDC-TR-78-68 


F(3. 31=0.00 
F(3.41=03 
F(3.51=C3 
GO  TO  34 

33  CONTINUE 

FC 3.11=1.00 

F( 3*2)=— 2.D04A3 

F( 3.31=1. DO 

F ( 3* 4 ) =—2.00*03 

F( 3 .51  =C3— A3**2— B3**2 

34  CONTINUE 
0Z21=Z2-Z1 
0Z31=Z3— Z1 
DZ32=Z3-Z2 
Z12=Z1«-Z2 
OO  35  K=1 . 5 

FB(Kl=(F(2  *K1 — F ( 1 * K1 1/0Z21 
FC(K)=((F(3.K)-F( l.K} 1/DZ31 -FB< K > 1/DZ32 
FBI K1=FB(K)-FC(K) *212 
FA(K1=F( 1 »K1  — IFBCK1 +FC C KJ  *Z1 1 *Z1 

35  CONTINUE 

N2=2.D0* (FA(  1 1 C FB ( 1 1 *FC(  1 1 *Z1  *Z  1 ♦X4FA  ( 2 l+CFBI  2 1+FC (21  *Z  1*Z 
N3=2.D0*(FA(3)*(Fa(3»*FC(3)*Z)*Zl*Y4FA(414(FB(4)4FC(41*Z>*Z 
N1 =— (FBI 1142  »O0*FC{t 1 *Z1*X**2-(FB(2)*2.D0*FC(2)*Z1*X 

* -CF0(3|42.OO*FC(31*Z1*V**2-(FBC4142.OO*FCC41*Z1*Y 

* 4(FB(5I42.00*FC(51*Z1 
AX=FAllH-(FB(i  H-FC(l)*Z>*Z 
IF ( AX. EQ. 0.001  AX=— 2.00 
AX=-2.D0*(FA(2)+<F8(21 *FC(2 1 *Z)*Z)/AX 
BV=FA(31*(FB(3)*FC(31*Z1*Z 

1F( BY. EQ. 0.001  BY=-2.D0 

B Y=— 2.00*( F A(  414-  (FBI  41  +FC  (41  *Z1  *Z1/BV 

N2=-DABS(N2l*OSlGN<l .OO.X-AXl 

N3=-DABS(N3t  *OS IGN(  1 • 00  .Y-BY 1 

GO  TO  26 

END 
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SUBROUTINE  OUTPUT  1 1 T) 

COMMON  /TITLE/HEADC20) 

COMMON/ VALUE /ZE • X ( 1 000 )«Y!1000) • PI 1000 ) .ft  HO  I 1000)  • 

TOC  1000) • THETA < 1000)  ,PSI (1000  > *N 
COMMON  /STAG/PTQC  1 000  I .TTO. ATOT.ZO. ZMAX.NBC200 )• 

♦ NBO!  200) • JMAX, JPT.lB.l PLOT, IT11. If 12 
COMMON  /G AS/G  ,RG .GM1  .GM1 H ,GM1G .GGM1 . GP1 • GPGM 
INTEGERS  NB.NBC 

1 FORMAT! *1 • .20X.20A4//30X, *ZE  »• * 1PE12.5, 15X. 'STEP  NO.* 
♦.I 5//2X,  *NO.  • .BX.'X*  .14X.  *Y*  .1 4X  • • P*  . 1 3X  • * RHO*  • 13X*  *Q* 
♦*12X. 'THETA*  .1 1 X • *PSI  • .13X.*M«  *8X  • • BOOT  * //  ) 

2 FORMAT  (1X.I4.1  P8E1  5.5  .1 X * A4 ) 

3 FORMAT! ' 1 * t20X .20A4//10X, * BODY  POINTS  AT  ZE  =• • 1PE12.5 

♦ ■1  5 X , *STEP  NO.  **I5//2X,'N0.  *.8X.*X*,14X,*V*  . 14X.'P*. 

♦ 13X.  *RH0*  .13X.*Q*,1  2X,  • THETA*  .11X.  *PSI*  . 1 3X.  • M • //  ) 

4 FORMA T( 1 X • I 4 ,1 P8E1 5 .5 ) 

DATA  NES/*  YES*/, NO/*  •/ 

DIMENSION  XPC1 0021 .YPC1 002) 

OATA  XO/O • 0/, OX/25. 0/*YG/ 10 .0/.DY/5.0/ 

IFC  I PLOT «EO.  0)  GO  TO  21 
F=q.o 

. DO  20  I— 1 »N 
xpc  n^F 

YPC I )=PC I ) 

F=F+1. 0 

20  CONTINUE 
XP(  N-t-1  ) = XO 
XP(N*2)=DX 
YPC  N+ 1 ) —VC 
YPC  N«-2)=DY 

CALL  A XI  SC  0.0, 0.0. 'POINT* ,-5 ,1 4 . 0 .0 • 0 . XO, DX . 1 0 .0 ) 

CALL  AXIS(0.0»0,0,*P'*1  .1  0.0  .90.0. YO.DY,  10.0) 

CALL  LINECXP.YP.N.l  .0.1  ) 

CALL  SYMBOL! 8.0  .8* 0.0.15. 'STEP  • .0.0*5) 

CALL  NUMBER! 999 .0,999.0 .0.15 .FLOAT (IT) ,0 .0.-1 ) 

CALL  CALC MPCO.O.O.O .0 ,2) 

21  CONTINUE 

IF{  ITH.EQ.O)  GO'  TO  22 

WRITE!  11 ) ZE ,N , I B. NB  * NBO, ! X f J } ,Y ( J) ,P ( J ) , RHO  ( J ) , OC  J ) » 
♦THETA! J) .PSICJY.PTC(J) . J=1 ,N1 

22  CONTINUE 

IFC IT12.EQ.0)  GO  TO  24 
WRITE! 12)ZE.N,18. NB . NBO 
DO  23  J=1 .IB 
K=NBQ( J) 

WRITE! 12) X(K) ,Y(K) ,P(K) «RHO(K),Q(Kl .THET A (K ) ,PS I C K ) 

23  CONTINUE 

24  CONTINUE 

IF(JPT.EQ.I)  GO  TO  25 

IF ! < IT.EQ. JMAX) .OR. !ZE.EQ.ZMAX) ) GO  TO  25 
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IF  C MOO ( ITt JPT ) • NE  *0 ) RETURN 

25  CONTINUE 
L-i 

DO  28  J=1 »N*50 
N2«J+49 

IFCN2.GT.N)  N2=N 
WRITE! 6.1 ) HE AO • ZE • IT 

00  27  K=J.N2 

SN=Of  Kl /SORT CG*PC  K) /RHOC Kl 1 

1 N=NO 

IF  { NBC L ) * NE • K ) GO  TC  26 
I N-NE S 
l_=L+ 1 

26  CONTINUE 

WRITEC6.2IK.Xt  Kl  .YIK)  ,P(K) «RHO(K)  . Q< K I «T HET A (K  I • 
*PSI  ( K ) * SM  *1  N 

27  CONTINUE 

28  CONTINUE 

WR I TE ( 6 • 3 I HEAD . ZE  , I T 
DO  29  J=l  .IB 
KsNBOl J) 

SH-Q{K)/SQRT(G*PtK>/RHO(K»  I 

WRITEC6.41K.XCK)  .VC  KJ  ,PCK)  .RHQCK)  ,Q  I Kl  .THET  ACK  I. 
♦PSI'C  K)  *5M 

29  CONTINUE 
RETURN 
END 
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SUBROUTINE  S I MQ  ( A ,R  , M.  N) 

IMPLICIT  REAL *8 ( A— H .O— Z ) 

DIMENSION  All)  .R(I) 

OATA  EPS/I.  D-5/ 

111  FORMAT! *0  SINGULAR  MATRIX  IN  SIMQM 
IF(M)23»23.1 

C SEARCH  FQR  GREATEST  ELEMENT  IN  MATRIX  A 

1 IER=0 
PIV=0.00 

mm=m*m 
nm=n*m 

DO  3 L=1 .MM 
T0=OABSCACL1  ) 

|F( TB-PI V)3.3.2 

2 PI  V=TB 
I=L 

3 CONTINUE 
TOL=EPS*PI  V 

A!  1 1 IS  PIVOT  ELEMENT.  Pf  V CONTAINS  /M  II/. 

START  ELIMINATION  LOOP 
LST=1 

DO  17  K-l  .M 
TEST  ON  SINGULARITY 
IF(PIV)23«23.4 

4 IF(IER>7,5.7 

5 IF ( PI  V— TOL) 6*6.7 

6 IER=K-1 

7 Pi VI^l.OO/A! I 1 
J=C 1-1 I/M 
1=1- J+M-K 
J=J+ 1—  K 

I+K  IS  ROW-INDEX.  J+K  COLUMN-INDEX  OF  PIVOT  ELEMENT 
PIVOT  ROW  REDUCTION  ANO  ROW  INTERCHANGE  IN  R. 

DO  S L=K«NM,M 
LL=L*I 

TB=PIVI *R ( LLI 
R(LL)=R(L> 

8 RfL)=TB 

C IS  ELIMINATION  TERMINATED 

IF t K—  Ml 9 .18.18 

C COLUMN  INTERCHANGE  IN  MATRIX  A 

9 LENO  = LST«-M-K 
1F< JII2.12.10 

10  II=J*M 

OO  11  L=LST  .LEND 
TB-A  I L ) 

LL=L4II 

A!L)=A<LL) 

11  A! LL)=TB 

C ROW  INTERCHANGE  AND  PIVOT  ROW  REDUCTION  IN  MATRIX  A 
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12  00  13  L=LST *MM  « M 
LL=LFI 

T0=PI VI *A(LL> 

A(LL>=A(LJ 

13  AtL)  = TB 

C SAVE  COLUMN  INTERCHANGE  INFORMATION 

A( LSTI 

C ELEMENT  REDUCTION  AND  NEXT  PIVOT  SEARCH 

PIV=0.00 
LST=LST+l 
J=0 

00  16  I I=LST  »LEND 
PI  VI  =— AC  II) 

IST=I I+M 

J=JF1 

DO  IS  L— I ST • MM • M 
LL=L- J 

A ( L 1 =A (L1+PIVI ♦A(LLI 
TB-OAB  S(  ACL)  ) 

IF  ( TB-PI  VI  IS.  15  *14 
1 4 PI  V^TB 
I=L 

15  CONTINUE 

DO  16  L=K*NM'N 
LL=L+J 

16  R(LL)=RILL> +P1 VI *R(L> 

17  LST=L5T+M 
END  OF  ELIMINATION  LOOP 

SACK  SUBSTITUTION  AND  BACK  INTERCHANGE 

18  !F(M-1)23«22.1 9 

19  CST=MM+M 
LST»M*1 
DO  21  I —2 «M 

1 I=>LST-1 
I ST=I ST-LST 
L=IST-M 
L=A(L)fOt SDO 
00  21  J=I I * NM • M 
TB=R{ J 1 
LL-J 

DO  20  K— I ST  » MM  «M 
LL=LL*1 

20  T0  = TB— AC  K1  ♦R(LL) 

K-J+L 
R{ J ) — R ( K) 

21  RCK)=TB 

22  RETURN 

C ERROR  RETURN 

23  IER=-1 

MR  I TE ( 6 • 1 1 1 ) 

STOP 

END 
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SUBROUTINE  SOL VBP< X «Y . A .B.C * A1 • B1 .CL • IT • X 1 • Y1 > 

IMPLICIT  REAL*8< A-H.O-Z) 

1 FORMAT  C * 0 NEGATIVE  RADICAL  IN  SOLVBP* ) 

GO  TO (21  .22  .23.24)  .IT 

21  CONTINUE 
X=C  I 

Y=(C-A*X)/B 

RETURN 

22  CONTINUE 
Y=Cl 

X=!C-B*Y)/A 

RETURN 

23  CONTINUE 

Y=(  A*C1-C*A1  ) /(A*B1-8*AI  ) 

X=( C-B+Y) /A 
RETURN 

24  CONTINUE 

IFI  DABSt  BI.GT.DABSt  A)  » GO  TO  26 

DI=l.O0+(B/A)**2 

D2=C/A-A1 

D3»(D2**2*B1  4*2— Cl  l/DI 
D2=( B*D2/A+B1 ) /D1 
R AD=D2#*2— 03 

IF(RAD.GE.O.OO)  GO  TO  25 
URKTEI6.1) 

WRITE(6.2)X.Y.A.B.C.A1.81.C1 .IT .XI .VI *D1 *02. 03 

2 FORMAT!  *0  X= * .1  PEI  2 .5 .5X • *V= * .El 2 .5/3X . • A=* *E1 2 .5* 5X. 
**B=»,E12.5  »5X,  •C=*.E12.5/2X.'Al=«.E12.5.4X.*  B=*  .El 2.5. 
*4X.  »CI  = « ,E12*5/2X.MT=*  *I2.14X*  *Xl=*  .E12.5.4X.  »V1*». 
♦E12.S/2X. ■D1=».E12,5,4X.'D2»'.E12.5.4X.*D3='.E12.5) 

CALL  ERR MCA 
STOP 

25  CONTINUE 
RAO-D SORT (RAD) 

V=O2«-RA0 

IFI OABS( Yl— Y)  • GE .DABS! VI— D2+RA0) ) V =02 -RAD 

X=ic-B*V)/A 

RETURN 

26  CONTINUE 
01-t.D0fr(A/B)4*2 
D2=C/B-Bl 

03*( D2**2*A 1 **2-Cl ) /01 
0 2=1 A 402 /B+A 1 1/01 
RA0=D2 **2-03 

IF{ RAD.GE.O.OO)  GO  TO  27  ^ 

WRITE! 6*1) 

WRIT€!6.2)X.Y.A.B.C.A1  .81  .Cl  .IT, XI. Yl  .0t.D2.03 

CALL  ERRWCA 

STOP 

27  CONTINUE 
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RAO  sO  SORT  (RAD) 

X=D2+RAD 

IF ( 0A8S( XI— X) * GT *0ABS(X1 — D2+RAO) ) X =02 -RAD 

Y=(C-A*X)/8 

RETURN 

END 
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SUBROUTINE  SORTCON*NE»N) 
DIMENSION  DNINI 
INTEGER *2  NEC  N) 

Nl  — t 
N2=N 

20  CONTINUE 
DMIN^DNC  Nl» 

DMA  X=DNC  N2) 

IFCDMIN.LE.DMAX)  GO  TO  21 
D-DNC  NI J 
DNCN1 )— DN( N2> 

DN{ N2) — D 
M-NEIN1) 

NEI  NI  ) *NEC  N2» 

NEC  N2  J *M 
GO  TO  20 

21  CONTINUE 
N1=N1+1 
N2=N2- 1 

IF I N 1 • GT • N2 1 GO  TO  24 
DO  23  I — N 1 • N2 

1FCDNC t ) .GE.DMIN)  GO  TO  22 

DN(Nl-I)-DN<n 

DN<  I )=DMI  N 

DMI  N=DNf  Nl-1  I 

M-NE C I » 

NEC  I I =NE ( Nl  — 1 ) 

NEC  Nl-I  >=M 
GO  TO  23 

22  CONTINUE 

IF(DN( n .LE.DMAXl  GO  TO  23 
DNCN2+M-ONC  I) 

0N<  I > -D  MAX 
DMA  X=DNC  N2+1I 
M=NEC I ) 

NEC  11=NE(  N2+1  > 

NECN2+n»M 

23  CONTINUE 

GO  TO  20  . 

24  CONTINUE 
RETURN 
END 
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SUBROUTINE  THREED 

COMMON/VALUE/ZE  .XllOOO)  »Y(  10001  .PI  1000 » .RHOI 1000)  • 
*0( 1000) • THETA I 1 000 ) .PSI (10001.N 
COMMON  /GAS/G.R0.GM1  »GM1  H. GMl  G.  GGM1  . GPt  . GPGN 
COMMON/N V ALU/D Z *XN( 1000) «YN( 1000 ) *PN 1 1000 ) »«N« 1000). 
*QN<  1000) .TNI  10001  .PSNI1000) 

COMMON  /BCUT/X0C5O) .YBI50) 

COMMON  /OCUT/X8OI50) .YBOISO ) 

COMMON  /STAG/PTOI 1000)  .TTO. ATOT • 20. ZMAX . NBC  200 )• 
*NBO( 200) • JMAX. JPT  .IB 
INTEGER*2  NB  ,NBO 

REAL  *8  ZD.DG .DRC.DGM1 .OGM1H.OGM1 G. DGGM1 . DGPl.OGPGM 
COMMON  /DG A S/ZD . DG .DRO.DGMl ,DGM1 H.DGM1 G. DGGM1 . DGP l * 
♦DGPGM 

REAL*S  PTDO.TTOD.ATOD 

COMMON  /O STAG /PTOD 1 1 000 ) .TTOO.ATOO 

OQ  19  1=1. N 

PTODI I )=PTCC  I ) 

19  CONTINUE 
TTOD=TTO 
A TOD -ATOT 
OG--G 
OSO=RO 

OGM l=OG— 1 .DO 
DGM1H=0.5DO*OGM1 
DGM1G=DGM1/DG 
DGGM1  = 1 .D0/DGM1  G 
DGP 1=DG* 1.00 
DGPGM=DGP1/DGMI 
CALL  CUTIZO) 

ZE=ZO 

DO  20  1=1. SO 
XBOI I )=XBC1 1 
YBOI I )=YB|I ) 

20  CONTINUE 

DO  26  J=1 . JMAX 
CALL  OUTPUT! J-ll 
CALL  OIST(DS) 

DZ=  0.  5*0 S 

IFKZE+DZI.LE.ZMAX)  GO  TO  21  ’ 

IF(ZE.GE.ZMAX)  RETURN 
DZ=Z  MA  X-ZE 

21  CONTINUE 
ZO=DZ 

CALL  CUTtZE+OZ) 

* = 1 

DO  23  1=1. N 
CALL  PIT! I ) 

IFINB(K) »EQ* I)  GO  TO  22 
CALL  FIELD! 1) 
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CO  TO  23 

22  CONTINUE 
CALL  BGDYC I I 
K=KM 

23  CONTINUE 
ZE-ZE+OZ 

DO  24  1=1  * N 
X(  I >=XN<I  I 
YC I I =YN{ I ) 

PI I I =PNt I ) 

RHO! I ) =RN<I > 

QI  I )=QN(  I J 
THETAC I >=TN( I ) 
PSII  I )=PSN!  I I 

24  CONTINUE 

DO  25  1=1 «SO 
XSO<  I )=XB(  I ) 

Y0Q! I 1 =Y8 ( I 1 

25  CONTINUE 

26  CONTINUE 

CALL  OUTPUT! JM AX) 

RETURN 

END 
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APPENDIX  B 
EXAMPLE  PROBELM 

THE  SAMPLE  PROBLEM  IS  FOR  M = 4.1  AXISYMMETRIC  NOZZLE  WITH  GAMMA 
= 1.2.  THE  NOZZLE  AXIS  IS  ALIGNED  WITH  THE  z AXIS.  THERE  ARE  48  STATIONS 
FOR  THE  BODY  GEOMETRY.  THE  STARTING  PLANE  INPUT  HAS  ONE  RAY  WITH  11 
POINTS.  THE  COMPLETE  STARTING  PLANE  IS  GENERATED  FROM  THESE  11  POINTS  TO 
A TOTAL  OF  346  NEARLY  EQUALLY  SPACED  POINTS.  THE  ARRANGEMENT  OF  THESE 
POINTS  IS  SHOWN  BELOW. 
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INPUT  CMOS 


OOflOAOOOO  111111111  1 2222»222233333333334l444«M«45fS95SSa9S66666666ftA7T7T7Trrm 
1234567fl9OI?14967A9A12349eTe9Ot234567B9At2345frre9O12345A70f9l23496789OI234MrAM 


•a 


36 

3 

0*205 

1*01532 

0*0 

36 

3 

0*25 

1.02299 

0.0 

36 

3 

0*30 

1 *0334* 

0*0 

36  . 

3 

0*35 

1*04594 

0.0 

36 

3 

0*40 

t. 06051 

0.0 

36 

3 

0*45 

1*07710 

0*0 

36 

9 

0*50 

1.09561 

0.0 

36 

3 

0 *ss 

1*11506 

0*0 

36 

3 

0*60 

1*13762 

0*0 

36 

3 

0*05 

1*16050 

0*0 

36 

3 

0*70 

1.10442 

0*0 

36 

3 

0*00 

1*23355 

0*0 

36 

3 

0*9(1 

1*20315 

0.0 

56 

3 

1.00 

1*33273 

0*0 

36 

3 

1 .10 

1*30226 

0*0 

36 

3 

1 .20 

1*43171 

0.0 

36 

3 

1.30 

1*40105 

0*0 

36 

3 

1.40 

1*53420 

0.0 

36 

1*50 

I *57936 

0.0 

36 

1 .60 

1 .62031 

0.0 

36 

1*70 

1*67711 

0*0 

36 

1.00 

1.72574 

0.0 

36 

1 *93260 

1 *76990 

0*0 

36 

2*09370 

OOOOOOOOO  I 1 I 1 1 1 1 1 1 1 2222222222333  33333334444444444555S5559S566666666G6777  777TTT7Q 
127456769012145670901 2345670901 £34 56769 01 23456 709012345670961234567690 123456 7090 
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INPUT  CARDS 


bO 


00000600011 11  lit  II  12222222222  J3J3333333444444444453S3353B55666666G6667T77777T778 
123456789012345678901234567890!  234*078901 2345670901 2345678901234567090 123456  7890 


1*00711 

0*0 

36 

3 

2*28244 

1 •95576 

0*0 

36 

3 

7*50080 

2*05662 

0*0 

36 

a 

2*75167 

2*16697 

0*0 

36 

a 

3*04016 

2*29026 

0*0 

26 

3 

2*37131 

2*42576 

0*0 

36 

3 

3*75138 

2*57362 

0*0 

36 

3 

4*18794 

2*73284 

0.0 

36 

3 

4*68976 

2*90607 

0*0 

36 

3 

5*26701 

3*06965 

0.0 

36 

3 

5*93200 

3*26331 

0.0 

36 

3 

6*69770 

3*46569 

0*0 

36 

a 

7*57930 

3.69340 

0.0 

36 

3 

8*59135 

3*90291 

0*0 

36 

3 

9*74737 

4.10686 

0.0 

36 

3 

11*0567 

4*3690! 

0*0 

36 

3 

12*5192 

4*48414 

0*0 

36 

3 

14*1196 

4*63912 

0*0 

36 

a 

15*8202 

4.76413 

0*0 

36 

3 

17*5562 

4,85610 

0*0 

36 

3 

19*2369 

4*91634 

0*0 

36 

3 

20*7800 

4*95043 

0*0 

36 

3 

22*0350 

4*96647 

0*0 

36 

3 

23*0095 

4*97238 

0*0 

36 

3 

24.0 

4.97416 

0*0. 

1*24 

100*0 

100*0  i?io«0 


000 00009 01 1 II  lit  1 1 l 22222 222223 33 33 333? 344444 4444455555555556666666666777777 7T7T6 
1234567890 123456789012345670901 234SA7890  1 2345678901 23456789012345676901234567690 
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INPUT  CHOPS 


00*0  #0 000  I 1 1 1 1 1 11 1 1 322222222233?  3333333*4444444*40590000599660066666677777777776 
12343670901 234S07B90I 2343670901 2345670901 23466709012345670901234367090 1234967090 


0*205 

24.0 

0 

I II  36 

100  10O 

0*0 

0.0 

59.035 

1.0105 

0*0  o.o 

0*0 

0 

0.10153 

0.0 

54*007 

1 *01 291 

0*22466  0.0 

0.0 

0 

0.Z0  306 

0.0 

54.439 

1.02022 

9*47544  9.0 

0*0 

0 

0.30460 

0.0 

53*676 

1.03271 

0. 77961  0*0 

0*0 

0 

0.40613 

0.0 

52*571 

1.05091 

1.16097  0.0 

0.0 

0 

Q. 00 766 

0*0 

51*079 

1.07571 

1.00158  0.0 

0*0 

0 

0*69919 

0.0 

49. 127 

I.I0B59 

2-36959  0.0 

0.0 

0 

0.71072 

0.0 

45.513 

1*15103 

5*30375  0*0 

0-0 

Q 

0-01226 

0.0 

43*985 

1.20005 

4.50134  O.O 

0 *0 

0 

0^91379 

0.0 

39. 245 

1.20600 

0*32552  0.0 

0*0 

0 

1.01532 

0.0 

39.953 

K .39203 

0*57014  0.0 

0*0 

1 

0000060001 tlltllll 1 2222222222333333333344444444443990059036666666666677777777770 
12345670901234567000123496769012343670901 2343676901 2340678901234067896 1234367690 
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MACH  * NOZZLE  <G«1,24> 


STARTING  PLANE  INPUT 
GANNA  *1,240  PTO  = 1.00000E  02 

ZD=  2,050  00C-0!  ZM A X=  2,4Q0Q0E  01 


X 


T 


P 


T TO  - 1.00000E  02 


1 • 71 600 E 09 


I Type  * 0 inv-  1 

I PLOT  = o ITU-  0 


NPTS*  II  MOATS-  » 

IT 12-  0 INEIGH-  0 


O 


THETA  PS1 


JM A X—  1 00  JPT- 

PTO 


100 
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MACH  4 NOZZLE 


(Gal. 241 


ZC  - 2*051006-01  STEP  NO.  D 


NO. 

X 

V 

P 

X) 

X 

o 

O 

THETA 

Pfil 

II 

1 

0.0 

0*0 

5.50350c 

01 

3*600156-04 

4,399526 

02 

0*0 

0*0 

1*010506 

00 

2 

1,01  5 JOE- 01 

0*0 

5*400706 

01 

3*592346-04 

4.40887E 

02 

3*92105E-03 

0*0 

1*012916 

09 

9 

3.0765IE-02 

0*792 746- 02 

5.400706 

91 

3*592346-04 

4.408B7E 

02 

1-940536-03 

3.39573E-03 

1*012916 

00 

4 

-5*078406-02 

0*792766-02 

5*400706 

01 

3*  592 34 E— 04 

4.40BB7E 

02 

-1*960526-03 

3*398746-03 

1*012916 

00 

9 

-1. 01530E-01 

1*606796-07 

5-4  00706 

01 

3 *592346— 04 

• *400876 

02 

-3.921  05E-03 

6*201556-09 

l«0l£9ie 

00 

6 

-3. 6765  IE-02 

— 0*792746—  02 

5.40070E 

01 

5*992346-04 

4*408876 

02 

-1.98053E-03 

-3*395736-03 

1 • 01 291 E 

00 

7 

5.07448E-O2 

-8*792  766-02 

5*480706 

01 

3.59234E-04 

4*4088  7E 

02 

1 .980 5 2 E— 03 

-3. 395746-03 

1 *012916 

00 

a 

2,03060E~01 

0*0 

5.44390E 

01 

3 *56  BATE —04 

4.43717E 

02 

B.29798E-03 

0*0 

1.02022E 

08 

9 

t.rssssE-oi 

1.91&306-01 

5.44390F 

01 

3*560676-04 

••4371 7E 

02 

7*1 882 3E— 03 

4.1490SE-03 

1 *020226 

00 

ID 

U41S30E-61 

I *750556-01 

5*443906 

01 

3*560676-04 

4*4371 7E 

02  . 

4*  14094C— 03 

7*186296-03 

1*020226 

08 

11 

2,074036-07 

2*030606-01 

5*  443  90E 

01 

3,560676*04 

4.43717E 

oz 

1 .05185E-08 

8*297906-03 

1.02022E 

oo 

12 

-1.015306-01 

1 *750556-01 

S.44390E 

0 I 

'3  *568676—04 

4*4371 7E 

02 

—4  » I4894E— 03 

7*186306-03 

1 *02022E 

00 

13 

-l *750556-01 

1*015306—01 

5*443906 

01 

3*568676-04 

4.43717E 

02 

-7*188236-03 

4* 149076—03 

1 *02 022 E 

00 

I* 

-2.0306OE-01 

3*211596-07 

5.443906 

01 

3.560676-04 

4*437176 

02 

-8.29790E-O3 

1.312446-00 

1 *02  0226 

00 

15 

-1 *755556-0 & 

- 1*015306-01 

S *443906 

01 

3*560676-04 

4*4371 7E 

02 

-7,188246-03 

-4*149056-03 

1 *020226 

00 

ID 

-i *oi33oe*oi 

-1.750 536-01 

5 *443906 

01 

3*  56 56 7 E -04 

4*4371 TE 

02 

-4*140966-03 

-7*106296-03 

1 *020226 

00 

IT 

-3*049166-07 

-2*030606-01 

5*443906 

01 

3*568676-04 

4*4371 7E 

02 

-1.57293E-88 

-8*297986-03 

t *020226 

oo 

IS 

1 .015306-01 

- 1 *750556-01 

5.44390E 

01 

3*  56B6/FH>4 

4.437176 

02 

4*140936-03 

-7*186306-03 

1 .020226 

00 

19 

1 . TS6S5E— 01 

- i«oiG3oe-oi 

5*4439 OF 

01 

3.56067E-O4 

4*437176 

02 

7*186236-93 

-4*149076-03 

1 *020226 

00 

20 

3*04fipO€-01 

0*0 

5*367606 

01 

3*528285-04 

4.4B53BE 

02 

1*  361026-02 

0*0 

1 *03271 E 

00 

21 

2 * 8S096E— 01 

9*890336-02 

5*367606 

01 

3*  52 02 BE— 04 

4.4053OE 

02 

1 *287206-02 

4 *419476— 03 

1 *032716 

00 

22 

2.403726-01 

1*070896-0! 

5.367606 

01 

3*  520  2 8E— 04 

4*405366 

02 

1 *074036-02 

0*359896-03 

1 *03271 E 

oo 

29 

1 « 6660  IE— 01 

2*5500 16-01 

5*3676  06 

01 

3. 529206-04 

4.48536E 

02 

7*44394E— 03 

1*139426-02 

1 *09271  E 

00 

24 

7.477546-02 

2*952  796-  01 

5*367606 

91 

3*528286-04 

4.4B536E 

02 

3.3 4104E-03 

1*319306-02 

1 *03271 E 

00 

2S 

-2*9193 16-02 

3.035606—01 

5*367606 

01 

3*528286—04 

4,485306 

02 

-1*123076-03 

1*356386-02 

1 *032716 

00 

26 

-1*223366-01 

' 2*709456-01 

5*3676  OE 

01 

3 * 520 28 E— 04 

4*4853 6E 

02 

-5*466996-03 

1*246406-02 

1*032716 

00 

27 

— 2*06  2046*01 

2*241 02E-01 

5*3676  OE 

01 

3*520286-04 

4.4S336E 

02 

-9.217706-03 

1 *001376-02 

1*032716 

00 

2a 

-2  .078076-01 

1.44974E-01 

5*  367606 

01 

3.520  28E-O4 

4*485 36E 

02 

-1*196976-02 

6,470106-03 

1*032716 

00 

29 

“3 • 004456—0 1 

9.01366E-02 

5*3670 OE 

01 

3 *528286-04 

4*405166 

02 

-1*342466-02 

2*240356-03 

1*032716 

00 

30 

-7*004466-01 

-5*013456-02 

5*367606 

01 

3*528286-04 

4*405366 

02 

-l. 342466-02 

-2*240266-03 

1,032716 

10 

31 

— 2*6700 0E— 01 

-1*449726-01 

5*5676 OE 

01 

3*  520 20 E— 04 

4*405366 

02 

-1.196906-02 

-6  *470  016-03 

1 .032716 

00 

32 

-2*069016-01 

— 2*  24 10 IE- 01 

5*367606 

01 

3*528206-04 

4*48  5366 

02 

—9  *21 7446—03 

-1*001366-82 

1*032716 

00 

33 

-t  *223586-01 

-2.709446-01 

5.3676QE 

01 

3*528206-04 

4*405366 

02 

-5*467086-03 

-1*248406-02 

1 *0327!  E 

00 

34 

-2.51SS9E-32 

-3*035596-01 

5.36760E 

01 

3*  52826 E— 04 

4*405366 

02 

— 1 *123976—03 

-1*356376-02 

1*032716 

00 

35 

7 *477256—02 

-2*992006-01 

5*  36  7606 

01 

3.92B26E-04 

4*405366 

02 

3*340916-03 

-1*319386-82 

1*032  716 

00 

36 

? . 66S99E-01 

-2.550026-01 

5*367606 

01 

3,520206-04 

4.40S36E 

02 

7 *443066—03 

-1.13943E-02 

1*032716 

00 

37 

2*  49371E-01 

- 1. 870916-01 

5.36760E 

Ol 

3*520206-04 

4*405386 

02 

L.07402E-02 

-0*359906—03 

1 « 03271 E 

00 

as 

2 *060956—01 

-9*890606-02 

5 *367606 

01 

3*520266-04 

4*405366 

02 

1*207276-02 

-4.419606-03 

1 • 03271 E 

00 

39 

4.061306-01 

0*0 

5.2S7K0E 

01 

3. 469 59 E- 04 

4*555236 

02 

2*040246-02 

0*0 

1 *050916 

00 

40 

3 *933716—01 

1*010005-01 

5*  25710C 

01 

3*469596-04 

4*555236 

02 

1*9761 36-02 

5.07461E-03 

1 « 05091 6 

00 

41 

3.55B94E-01 

2*956546-01 

S.Z5TI0E 

01 

3.46959E-04 

4.5S523E 

02 

1 *707846-02 

9.82995E-03 

1*054916 

00 

42 

2,960666-01 

2,700156-  01 

5*2571  OE 

01 

3.46959E-04 

4*555236 

02 

1 *487226-02 

1*396746-02 

1*050916 

00 

49 

2*176156-01 

3.429076-01 

5*  2571 OE 

01 

3.469S9E-04 

4*555236 

02 

1*093166-02 

1*722706-02 

1*05091 E 

00 

44 

1*255026-01 

3*842526-02 

5*2571 OE 

01 

3*469596-04 

4.55523E 

02 

6*304306-03 

1,940406-02 

1 *090916 

00 

45 

7.5S016E— 02 

4*053296-01 

5.2571  OE 

01 

3*469596—04 

4*555236 

02 

1*281016-03 

2*0382 IE— 0? 

1*050916 

00 

46 

-7.61006E-02 

3*984366-01 

5.25710E 

01 

3*  48959 E— 04 

4*555236 

02 

-3*022746-03 

2*004116-02 

1*05091 E 

00 

47 

* 1 • 7 29 2 IE— 01 

3*674786-01 

8*2571 OE 

01 

3*469596-04 

4*595236 

02 

-8 *686386- 03 

1 *048 1 IE-02 

1 *050916 

00 

40 

-2  *560776— 01 

'•129296-01 

9*2571 Oe 

01 

3*469596-04 

4.55523E 

02 

-1*300446-02 

1*572126-02 

1*050916 

00 

49 

-3.285666-01 

2*3  07186-  01 

5*2571 OE 

01 

3*469596-04 

4.55523E 

02 

-1 *650546-02 

1*149336-02 

1*050916 

00 

50 

— 3*  7741 OE— 01 

1.49507E-01 

5.2571 OE 

01 

3 * 489 59 E— 04 

4 .55523E 

02 

-1*896946-02 

7.51184E-03 

1-050916 

00 

BODV 


AEDC-TR-70-6B 


MACH  4 NOZZLE  (01.241 

ZC  - 2.05000E-01  STEP  NO.  0 


NO. 

X 

V 

P 

PHD 

0 

TkETA 

PS! 

M 

51 

-4.029276-01 

5.090226-02 

5*257106 

01 

3.46959E-04 

4.SSS23E 

62 

- 2 *024146— 02 

2*557476-03 

1 *050916 

00 

B2 

-4.02928E-01 

-5*09009E-02 

5*  25710E 

01 

3*469596-04 

4*555236 

02 

-2*024156-02 

— 2*  5674 IE— 03 

1 *05091  E 

00 

53 

-3.776L0C-01 

-1.49S06E-01 

S.2S710E 

01 

3. 46950 E— 04 

4.55S23E 

02 

-1.89695E-02 

-7*511486-03 

1 .050916 

00 

54 

-3.28567E-01 

-2. JB7| 6E-01 

S.2S710E 

0! 

3. 4 69 59 E— 04 

4 .55523E 

02 

-I.65OS5E-02 

— 1 *199336—92 

1 *050916 

00 

55 

-2.589766-02 

— 3*1 2928C— Ot 

5.2STLOE 

01 

3.46959E-04 

4 *555236 

02 

-1.300446-02 

-1*572116-02 

1 *050916 

00 

56 

-1 .72922E-01 

-3.6Z477E-01 

5*2571 06 

01 

3.46939E-04 

4.55523E 

02 

-0*6B644E-OJ 

— 1 *046106—02 

1 *050916 

00 

57 

-7.61023E-02 

-3.96930E-O1 

5*257106 

01 

3*46959E— 04 

4*555236 

02 

-3*022026-03 

-2.004116-02 

1.0&Q91E 

00 

5S 

2*549926-02 

-4.0S729E-01 

5* 23710E 

01 

3.469S9E-04 

4.555236 

02 

1 *280096-03 

— 2 *0362  IE— 02 

1*050916 

00 

59 

1 * 25496E—  0 1 

.-3.H6253E-01 

5*2571  OF 

01 

3*469596-04 

• *535236 

02 

6*304146-03 

-1.94041E-02 

1 .050916 

oo 

66 

?.I7414C-0t 

-3.429006-01 

5*2571 OE 

01 

3.469596-04 

4*5552 3E 

02 

1.093156-02 

- 1.722706-02 

1.050916 

00 

61 

2*960546-01 

- 2*  700 17E-01 

5*  2571 OE 

01 

3 ■ 469 59 E— 04 

4*553236 

02 

1 *487216-02 

-1 *396756-02 

1 *050916 

00 

62 

3 .S5694E-0I 

- 1,95656E-01 

5-25710E 

01 

3*464596-04 

4.555236 

02 

1*787046-02 

-9.83001E-03 

1 .050916 

00 

6 3 

?4  93370E— 01 

-1*010026-01 

5.2571 OE 

01 

3.469B9E-04 

4.555236 

92 

1*976136-02 

— 5 .074C2E-03 

1 .050916 

03 

64 

5.07660E-01 

0*0 

S.10790E 

01 

3*389906-04 

4*649756 

02 

2*934916-02 

0*0 

1.07S71E 

00 

65 

4 • 9T268E-0 1 

1.02191E-01 

5.10790E 

01 

3.38996E-04 

4.649756 

02 

2 *0740 IE— 02 

5.9095SF-03 

1*075 TIE 

00 

66 

4. 0CS1 8E— 01 

?*0019BE-01 

5.10790E 

01 

3*389966—04 

4 .6*9756 

02 

2*697006-02 

I.15T68E-02 

1*075716 

00 

67 

4*1*6096-01 

2*900106-01 

5.10 7906 

01 

3*  38996 E— 04 

4*649756 

02 

2.40875E-02 

1*074946-02 

I .075716 

00 

66 

3.49701E-O1 

•'  3.67946E-0I 

3* 1079 OF 

01 

3.309966-04 

4 *649756 

02 

2*021906-02 

2.12T49E-02 

1.075TIE 

00 

69 

2*  68534E— 01 

4»308£2E-01 

5*  t 0790E 

01  ■ 

3*369966-04 

4*649756 

02 

1 *652  306-02 

2.499096-02 

1 .07571 E 

00 

70 

1 .703136-01 

4*760596-01 

5. 1O790E 

01 

3*389966-04 

4.64975E 

02 

1 *0191  BE— 02 

2 * 752316  — 02 

1*075716 

00 

71 

7.6B746E-02 

S*  O10O0E— 01 

5*107906 

01 

3*  389966—04 

4*649736 

02 

4.44368E-03 

2*9010^-02 

1*075716 

oo 

72 

-2*571 18E-02 

5*070062-01 

5* 1 0790F 

01 

3. 389966-04 

4*649756 

02 

-1 *406256-03 

2*931146-02 

1 .07571 E 

00 

73 

-1*272456-01 

4*91454E-01 

5. IO790E 

01 

3*309966—04 

4.649T5E 

02 

-7*335306-03 

2*841276-02 

1 . 0757! E 

00 

74 

-2*235706-01 

4.0S76OE-D1 

5*1 0790E 

01 

3*389966-04 

4*649756 

02 

— 1 *292366—02 

2.6351 2E-02 

1 * 0757 | E 

00 

75 

— 3#  1074 IE-0  1 

4*014468-01 

S* 1079 OE 

01 

3 • 389  96 E— 04 

4.64975E 

02 

-1 *796316-02 

2.321 1 IE-02 

1*075716 

oo 

70 

- 3*85*19 1E-0  1 

3*306766-01 

5* 107906 

01 

3.38996E-04 

4.649756 

02 

-2.2Z6T5E-02 

1 *912046-02 

1.0T57IE 

00 

77 

—4  *4 30 TOE— 01 

2*403696-01 

5*107906 

Ot 

3 . 389 96 E— 04 

4.649756 

02 

-2*906046-02 

1 *424636-02 

1 . 07571 E 

00 

76 

-4*843782-01 

I . S1976E-01 

5*1 0790E 

01 

3 . 38996  E — 04 

4 .6*9756 

02 

-2.000276-02 

8*708396-03 

1.075716 

00 

70 

-S.Q5055E-01 

S*  13601E-02 

5* I0790E 

91 

3 * 38996E— 0 4 

4.64975E 

02 

-2.919856-02 

2*970106-03 

1 *07571 E 

00 

60 

-S.03055F-01 

-5*135616-02 

5*107906 

01 

3*309946-04 

4 *649756 

02 

— 2*9 1 9 BSE—  02 

-2*969986-03 

1*075716 

00 

61 

—4*  64379E— Oi 

-1* 519746-01 

5*107906 

01 

3. 5094*6-0* 

4.649  756 

02 

—2*800286—02 

-0.70027E-O3 

1 *07571  E 

00 

62 

-4*438716-01 

-2*46367E-01 

5- 10790E 

01 

3. 3B996E-04 

4.64975E 

02 

-2*560046-02 

-1.424026-02 

1 *075  716 

00 

63 

-1.00I92E-O1 

-3.306756-01 

5.107906 

01 

3*  309 96 E— 04 

4.64975E 

02 

-2*22475E-0Z 

-1.912036-02 

1.07571E 

00 

64 

-3*107426-01 

-4.0I445E-01 

5*107906 

01 

3. 3B996E-04 

4.64975E 

02 

-1*796326-02 

—2*321 1 0E-02 

1*075716 

00 

65 

-2.238727-01 

-4.55779E-0I 

5* 1079 OE 

01 

3, 38996E-04 

4.6497SE 

02 

-1*292376-02 

— 2*4351 2E— 02 

1 • 07571 E 

00 

86 

-1.2T247E-01 

-4*914546-01 

5.10790E 

01 

3*309966-04 

4*649756 

02 

-7.J5547E-03 

-2. 04 1276-02 

1*075716 

00 

6? 

— 2*571 44E— 02 

-5.0T008E-01 

5*107906 

01 

3 ■ 30996 E— 04 

4.64975E 

02 

-1 *406400-03 

— 2 *931 1 4£— 02 

1.07571  E 

00 

66 

7*687106-02 

-5*010066-01 

5*  L0790C 

01 

3. 3690C E -94 

4.6497SE 

02 

4.44348E-03 

—2  *901096—02 

1 *075716 

00 

69 

1*765112-01 

- 4*  760006-  01 

5.107906 

ot 

3*369966-04 

4*049706 

02 

1*019176-02 

- 2. 752 32 E — 02 

1*075716 

00 

90 

2.68S3ZE-01 

-4.30824E-01 

5 * 107906 

o: 

3*369966-04 

4 .644756 

02 

1*552296-02 

-2*490906-02 

1*075716 

00 

91 

3.49760E-01 

-3*679406-01 

5* 10790E 

01 

3*  389966-04 

4*049756 

02 

2.02190E-02 

-2*127506-02 

1 *075716 

00 

92 

4*  16667E-01 

-2*9O0K2E-0l 

5.I0790E 

01 

3*  30996 C— 04 

4.649756 

02 

2 *400746—02 

-1*678906-02 

1*075716 

00 

03 

4*  0051 TF— 01 

-2*002006-01 

5*107906 

01 

3. 389966—04 

4*649756 

02 

2.696996-02 

-1*157496-02 

1 • 07571 E 

00 

94 

4*97266E-0l 

-1.02194E-01 

5*107906 

01 

3.3B996E-04 

4*649796 

02 

2.87401C-O2 

-5*909706-03 

1 *075716 

oo 

95 

6.09190E-01 

0*0 

4*912706 

01 

3.2BS08E-O4 

4*773846 

02 

4*135366-02 

0.0 

1.1O099E 

00 

96 

6*00  06  IE— 01 

1*002696-01 

4*912706 

01 

3*205006-04 

4 .773046 

02 

4*076936-02 

6.81035E-03 

1 *108596 

00 

07 

5*761826-01 

1 *976036-01 

4.912706 

OK 

3.205O9E-O4 

4.773046 

02 

3*911186-02 

1*3434^-02 

1*100596 

00 

96 

5. 33 70 7E -01 

2*699426-01 

4*912706 

01 

3.28508E-04 

4 *773046 

02 

3.636716-Q2 

1*949006-02 

1 * 1 08596 

00 

99 

4.A0737E-01 

3*  741  726-01 

4.91270E 

01 

3 • Z6S08E— 04 

4.773046 

02 

3*203036-02 

2.540096-02 

1 *100596 

00 

100 

4*1 25936-01 

4.46196E-01 

4*912706 

01 

3.2650BE-O4 

4*7730 4E 

02 

2 .800376-01 

3*043286-02 

1.1O0S96 

00 

0007 


AEDC-TR-78-68 


MACH  4 MUZZLE  <G»L.2*| 


ZE  » 2,  050906-01  STEP  NQt  Q 


HO. 

X 

r 

P 

rhO 

O 

T«TA 

PSI 

M 

191 

3 . 331 96E— 01 

5.  999936—01 

*.912706 

01 

3. 2a SO 86-0* 

* .7738*6 

02 

2 .201306—02 

3. 4625 BE— 02 

L. 106596 

08 

192 

2. *47096-01 

5.5788 DE— 01 

4.9I2T0E 

01 

?■ 28S0AE— 0* 

4.77384E 

02 

I.66076E-02 

3.787406-02 

1.108596 

80 

199 

1 .495*96-01 

5.  905*96-01 

*•912706 

01 

3. 285606—04 

4.7738*6 

02 

1 .0 1*906— 02 

4.008906—02 

1.106596 

00 

19* 

5,030706-02 

6.  071096-91 

4.91270E 

01 

3*  285086—0* 

*.7738*E 

02 

3,41*026-03 

4.12I26E-Q2 

1.108996 

00 

105 

-5.O3O01E-O2 

6.07109e-01 

*.912706 

01 

3 .285086—0* 

4.773846 

02 

-3.413966-03 

4.12I26E-02 

1 .1O059E 

oo 

IQfi 

-1* *95*69-01 

5.905496—01 

*• 91270G 

Ol 

3.285086-0* 

*•7730 *E 

02 

-I.014B9E-0Z 

4.008966-02 

1.108596 

00 

19T 

-2.**70fl€-01 

5.528896-61 

4.912706 

01 

3. 2 830 8 E— 04 

*•7738 *E 

02 

—1 .66075E— 02 

3. 7674 16—02 

1.10659E 

00 

10B 

— 3 • 331956—01 

5.0999*E—01 

*•912706 

01 

3. 285086-0* 

*•7738*6 

02 

-2,201376-02 

3,462586-02 

1,108996 

00 

too 

-*•125939*01 

4.481  966-  01 

4.912706 

01 

3.285086-0* 

4.7738*6 

02 

-2.Q0037E-02 

3.0*3296-02 

1.LO059E 

oo 

119 

— *• 60736E— 01 

3.7*1 72E-OI 

4.91  2706 

01 

3.28508E-04 

4.7738 *E 

02 

—3.283036—01 

2.5*0906-02 

1.108596 

00 

1 11 

-S.3S760E-OI 

2.899*36-01 

4.912706 

01 

9.285086-0* 

4,7738*6 

02 

-3, 836706-02 

1.969096-02 

1 .108596 

00 

112 

-5.70182E-O1 

1.976956-91 

*•91 2706 

01 

3.205086-04 

*•7738*6 

02 

-3*911186-02 

1.3*3446-02 

1.108596 

DO 

113 

-6.00«BtE-at 

1.002706-01 

*•912706 

01 

)• 285066-6* 

*•7738*6 

02 

-4 .078936- 02 

6.010*26-93 

1.108596 

00 

11* 

—6 • 99 190 E— 01 

9.63492E-97 

*•912706 

01 

. 3.2850BE— 0* 

*•7738 *E 

02 

-4.13538E-02 

6.54*216-06 

1.10899E 

00 

115 

—6 • 008S2E— 01 

-1  •902  666-01 

*•913706 

01 

3.28508E— 0* 

4.773B4E 

02 

-4.07893E-02 

-6,810306-03 

1.I0B54E 

00 

116 

-9.701836-01 

- 1.976026-01 

*•912706 

01 

3.285086-0* 

4.7730*6 

02 

-3.911 166-02 

-1.343*36-02 

1 .108596 

00 

IIT 

-S.3S7O7E-01 

-2.  899*16-01 

4.912706 

01 

3.28508E-0* 

4.773B4E 

02 

— 3. 80071 E— 02 

-1.909006-02 

1.10B59E 

00 

118 

-*•097375^-01 

-3.74!  716-01 

4.912706 

01 

3*  2850 8 C— 04 

4.7738*6 

02 

—3.2830*6—02 

-2.5*0396-02 

1 .100596 

00 

119 

-4.1259*6-01 

-*•*81956-01 

4.912706 

01 

3.28508E-0* 

4.77364E 

02 

-2.800386-02 

-3.0*3286-02 

1.I0099E 

00 

129 

—3  >3319  7E— 0 1 

— 5. 099926—01 

4*91 2706 

01 

3*  285086-0* 

* .7738*6 

02 

-2.261396-02 

-3 .4026 76-02 

I .1 08S9E 

00 

121 

”2 ♦ **7 10E- 01 

-5.578796*01 

4.91270E 

01 

3. 285086 -0* 

* .7738*6 

02 

—1 •660776—02 

-3.767*06-02 

1*108596 

00 

122 

-1 .*95*95-91 

-5.905*96-01 

4.912  706 

01 

3. 28508C— 0* 

*•7738*6 

02 

-1.01*916-02 

-4.008966-02 

1 . 1 0B59E 

00 

129 

-5.031036-02 

-6.071 096-01 

4.912706 

01 

3.285086-0* 

4.7738*6 

02 

-3. *1*2*6-03 

-4.121256-02 

1 .109846 

oo 

12* 

5.030S1E-02 

-6.071  096-01 

4.912706 

01 

3.285086-04 

4.7738*6 

02 

3.413896-03 

-♦.121266-02 

1.106996 

00 

125 

1 . 49S45E-0I 

-5.905*96-01 

4«9l2T0e 

01 

3.285086-0* 

4.7738*6 

02 

1 .014896-02 

-4.008966-02 

1. I 08596 

00 

126 

2*447971-01 

—5. 5 78 8 IE— 01 

4.912706 

01 

3*285086-04 

* . 773B*E 

02 

1.6007*6-02 

-3.787*16-02 

1.106596 

00 

127 

3 * 33192E  — 01 

— 5*099956-01 

4, 912706 

01 

3.2B508E-0* 

4.773B4E 

02 

2.261366-02 

-3.462596-02 

1.I06S9E 

00 

129 

*« 125906-01 

-*.*81986-01 

4.912706 

01 

3.285086-0* 

*•7738*6 

02 

2.800356-02 

-3. 043306-02 

1 .1 00596 

00 

129 

*•007356-0 1 

-3.  7*1736-  01 

4.91270E 

01 

3,289086-0* 

4.7738 *E 

02 

3.26302E-02 

-2,5*0906-02 

1 .10BS9E 

oo 

139 

5.357666-01 

-2. 899**6-01 

4.91270E 

01 

3.2BS08E-0* 

ft .7738*6 

02 

3.636706-02 

- 1.9691 06- 02 

1.1 06596 

00 

131 

5. 76  18 1C— 0 1 

-1.978076-01 

4 #91 2 TOE 

41 

3i  2Q508E-0* 

*•7738*8 

02 

3.911 I7E-02 

- 1.3*3*66-02 

1.I0BS9E 

00 

192 

6.008816-01 

-1.002736-01 

4.912706 

01 

3. 28508E-O* 

4.7798*6 

02 

4.0789 3E— 02 

-6.810016-03 

1 .106596 

00 

193 

7.10720E-01 

O.o 

4. 661 106 

01 

9.1*8836-0* 

*.93*896 

02 

5. 706 | 3E“ 02 

0.0 

1 .151836 

00 

S3* 

7.03486E-01 

1.011*66-91 

4*  661306 

01 

3. 1*8836-0* 

4.93489E 

Q2 

5. 707376-02 

8. 21*976-03 

1.151836 

00 

135 

6.8193  IE-01 

2.002336—01 

*•861306 

01 

3. 14SS3E-0* 

*.934896 

02 

5 .33 23 26-02 

1 .626166-02 

1 .101836 

00 

136 

6.*649*E-0I 

2.952*36-01 

*•001306 

01 

3.I4883E-0* 

4.93489E 

02 

5.2*4356-02 

2.397036-02 

1 .151836 

00 

1 37 

5.070966-01 

3.042**6-01 

*•661306 

01 

3.  1*9036-0* 

*.9  3*89E 

02 

4.8*9996-02 

3.1 I985E-02 

I .151836 

00 

138 

5.371 266-91 

* .554  225- Ot 

*•661306 

01 

3.14B83E-0* 

4.934B9E 

02 

4.356716-02 

3.77840E-O2 

1.151836 

00 

13* 

*•65*236-91 

5.3712&E-Q1 

*•061306 

01 

3. 1*0836-04 

4.4  3*096 

02 

3.T7462E-02 

4,359826-02 

t .151836 

00 

1*0 

3.842**6-01 

9. 978956-01 

4.00130E 

01 

3.1*8636-0* 

4.9  3*89E 

02 

3.116186-02 

4. 052346-02 

1.I9163E 

09 

1*1 

2.952**6-01 

6. 46*936-01 

4- 661  3 06 

01 

3.148B3E-0* 

4.Q34B9E 

02 

Z, 39*2*6-02 

5. 2*6096— 02 

1.151836 

00 

1*2 

2*0023 JE— 01 

6*61 93 16- 0 1 

4.86130C 

01 

3.1*8636-0* 

4.93489E 

02 

1.023O8E-O2 

9.5330*6-02 

1.15183E 

00 

1*3 

1 .011*76-01 

7.034  856-01 

4.66130F 

01 

3. 1 *8636-0* 

4.9  34B9E 

02 

8. 20 16SE-Q  3 

S.707S6E-02 

1 .19183E 

00 

1*4 

9.0*9226-07 

7.107206-01 

4.661906 

01 

3. 140836—04 

4.93489E 

02 

7,309196-08 

5.700136-02 

1.151036 

09 

1*5 

-1.01 1*5E— 01 

7.03*066-01 

*•861306 

01 

3.1*8696-0* 

♦•93*896 

02 

-6.201506-03 

9.707976-02 

1 ,151836 

00 

1*6 

— 2 . 00  23  7E— 0 1 

6*819316-01 

4.661305 

01 

3.14883E-0* 

4.93*896 

02 

-1.62367E-02 

5.533056-92 

1,151036 

00 

1*7 

r2» *524 36—91 

6. *64946-01 

*•661  30c 

01 

3* 1 *863 E—0* 

4.93*896 

02 

“2.39*236—02 

5.2*6066-02 

1.I9183E 

00 

1*9 

-3.0*2*36-01 

9.978966-01 

4.661306 

01 

3.1*8896-6* 

4.93*896 

02 

-3.116176-02 

4.652356-02 

1 ,151636 

09 

1*9 

—* .65*226— Ql 

5.371  27E-01 

*•66130 E 

01 

3.1*8096-0* 

4.9  3*896 

02 

-3.77*816-02 

4.359826-02 

1.1&I83E 

90 

ISO 

-5.371266-91 

4.65*236-01 

4.66130E 

01 

3.1*8036-0* 

*•93*896 

02 

-4 • 356 TIE- 02 

3. 770*16-02 

1,151836 

00 

✓ 


BODY 


89-8i-UX-OQ3V 


MACH  A NOZZLE 

ZE  - 2.0S000E-01  STEP  NO*  0 


NO* 

X 

V 

P 

PHD 

0 

THETA 

PS1 

M 

IS! 

-5.97S95E-0I 

3 •84?45E- 01 

4*6613  OE 

01 

3.14863E-04 

.9J4S9E 

02 

—4. 849986—02 

3.119665-02 

15163E 

00 

152 

-A*46493E-01 

2*952455-01 

4.661 30F 

91 

3. 1 4883 E— 04 

.9 3489E 

02 

-5.244SSE-92 

2*397545-02 

151835 

00 

153 

— 6*81930 E— 01 

?*00234E~01 

4*66! 30E 

01 

3.146835-04 

.93469E 

02 

-5.532316-62 

1.62C 185-02 

151 835 

00 

154 

-7. 0J4B6E— 01 

1. 011478-  01 

4.661 ?oe 

01 

3.14083E-04 

•93469E 

02 

—6*7073  75—  02 

• B.2I504E-D1 

1 51 63E 

00 

1 55 

-7.10720E-0! 

t « 1 2*07E~  06 

4.66130E 

91 

3. 14SO3E-04 

• 9 34B9E 

02 

-5*7661 3E-0 2 

9.1 29805-08 

15183E 

00 

136 

-7.07456E-01 

-1.01  145E-01 

4.66I30G 

01 

3.  14  053E-O4 

,934895 

02 

-5.707375-02 

-B.21486E-03 

IS183E 

00 

157 

-6*6193 IE- 01 

-2.O0232E-01 

4.66130E 

01 

3.I4S83E-04 

-93489E 

02 

-5*532326-02 

-1.62C16C-02 

15163E 

00 

156 

-6.A6494E-0I 

— 2*952435—  01 

4.661305 

0! 

3.14683E-04 

•934B9E 

02 

-S. 244555-02 

-2.39  75  3E-02 

1 SI  835 

00 

l&Q 

-S.970965-01 

-3.84243E-0S 

4.66120E 

01 

3*1 4683E— 04 

•93489E 

02 

-4.64999E-02 

-3.119B4E-02 

151835 

00 

160 

-5.37127E-01 

— 4 *654 21 E— 01 

4*661305 

01 

3.M883E-04 

•934B9E 

02 

-4*356725-02 

-3.776395-02 

I5I83E 

DO 

161 

-4.6S424E-0! 

-5*37| 25E— 01 

4.66130F 

01 

3.148035-04 

•934B9C 

02 

-3.774625-02 

“4 .359815-02 

151836 

00 

162 

-3.R4245E-01 

-5*976956- 01 

4.661 30E 

01 

3.1 4083E— 04 

•9  34695 

02 

' -3*1 16195-02 

—4. 85 23 46— 02 

151636 

00 

163 

-Z.9S24SE-01 

-6.464935-01 

4*6*1305 

01 

3*148635-04 

.93489E 

02 

-2. 39425E-02 

-5.24605E-02 

1 51 03E 

00 

164 

-2.00234E-01 

- A.0193OC-O1 

4 *661  305 

01 

3* I4883E-04 

.93489E 

02 

-1.623695-02 

-5*533046-02 

15183E 

00 

165 

-1 .9! lASE-Ol 

-7.OJ485E-0I 

4.661.305 

01 

3.I4B83E-04 

*9  34695 

02 

-6.201795-03 

—5. 707365-02 

151835 

00 

166 

-4.73621E-06 

— 7* 10720E— 01 

4.66130EL 

01 

3.14863F-04 

•Q3489E 

02 

-3. 840395-07 

— 5. 76C 1JF-02 

15163E 

00 

167 

1 .01142E-01 

-7.03486E-01 

4*661305 

01 

3.1 4883c— 04 

•9  3469E 

02 

0 . 20 1 JOE— 0 3 

-5.70737E-02 

IS1835 

oo 

16  5 

2*09230001 

- 6*  81 9325—01 

4 *661 3 OE 

01 

3.14B93E-04 

.9  3489E 

02 

I.6236SE-02 

-5.53305E-Q2 

151635 

00 

169 

2.9S242E-01 

— 6*  464945—01 

4. 661 3 OF 

01 

3. I486 35-04 

*934895 

02 

2.394225-02 

-5.24606E-O2 

1 SI  835 

00 

176 

84241 E-01 

-5.9769 6F-01 

4.661 305 

01 

3*148835-04 

*934895 

02 

3*1  16156-02 

-4. 852366-92 

161635 

00 

171 

4.65A26E-91 

-3.37129E-01 

4.661305 

01 

3.14863E-04 

,9 3 4895 

02 

3.774796-03 

-4. 399445-02 

1S183E 

00 

172 

5*371252-01 

— 4*654245—  01 

4.6613QE 

01 

3- 14AA3E-04 

.93489E 

02 

4.35670E-02 

-3.77B42E-02 

1S1B3E 

00 

173 

5*  970940  01 

- 3 « B4746E- 01 

4.661 3 QC 

01 

3.148836-04 

*9  34S9E 

02 

* * 8499 fl E- 02 

-3. 119666-02 

15IB3E 

Oft 

174 

6^464922-01 

— 2 * 95 2 4 BE—  9 1 

4.661305 

91 

3.140O3E-04 

.934B9E 

02 

5.244545— OZ 

-2.301575-02 

151835 

00 

175 

6*81930E-01 

— 2«  002 3 6F— 01 

4.66130E 

01 

3.14B83E-04 

.93469E 

02 

6 .5323 IE— 02 

-1.626195-02 

I5163E 

00 

1 76 

7 • 034855-01 

— 1 • 01 1 49E-  01 

4.661 30F 

01 

3. 14883E-04 

•934895 

02 

5.70737E-02 

-8.215195-03 

151635 

00 

177 

8. 12260E-01 

0.0 

4.3J651E 

01 

2.971 7 6£ -04 

5 

•144205 

02 

7.9999 35-02 

0.0 

209055 

00 

178 

8.05S55E-01 

l*01803E-01 

4.33850E 

01 

2.971765-04 

5 

•14420E 

02 

7.93274E-02 

1.00*266-02 

209035 

00 

179 

7.86741E-01 

2.02001E-01 

4.338505 

01 

2*971 76 E -04 

5 

.144205 

02 

7*744205-02 

1*992*95-02 

209055 

00 

ISO 

7*352205-0 1 

2.990 13E-01 

4.338S0E 

01  . 

2.971  76E-904 

5 

.144205 

02 

7.43339E-02 

2.945935-02 

20905E 

00 

151 

7*11769E-r01 

3,913095-01 

4. 33  55  OF 

01 

2.971 76E-04 

5 

.144205 

02 

7.009145-02 

3. 650386-02 

209055 

00 

182 

6*571325-01 

4.77434E-61 

4.  33  85  PE 

01 

2.971765-04 

5 

.144205 

02 

6*466456-02 

4.70C49E-02 

2090SE 

00 

183 

S.92112E-01 

5.560305-01 

4*338505 

01 

Z*97| TOE -04 

5 

•I4420E 

02 

5 .825 876—03 

5.479795-02 

209055 

00 

184 

5* 177S4E-0L 

6.  258575-01 

4.338505 

01 

2.971 765-04 

5 

.144205 

02 

5.093576-02 

6.166315-02 

20905E 

00 

185 

4 * 3B231E-01 

6.65613E-01 

4.338S0E 

01 

2.971  76E-04 

5 

•14420E 

02 

4.281185-02 

6.7BS32C-Q2 

209055 

00 

166 

3.4*8445-0 i 

7.349S4E-01 

4.33050C 

01 

2.971 76E -04 

S 

,|4420E 

02 

3.401545-02 

7*237725-02 

209058 

00 

167 

2*510035-01 

7*  7250SE-  01 

4.33AS0E 

01 

2.  97176E-04 

9 

•144206 

02 

2.468516-02 

7.6061 26—02 

20905E 

00 

166 

1 « 522O3E-0 1 

7.97872E— 01 

4.33850E 

01 

2.9T176E-04 

5 

•14420E 

02 

1.496765-02 

7.85466E-02 

2090SE 

00 

159 

9.1AQ32E-02 

8*106575-01 

4.33650E 

01 

£.971765-04 

5 

•14420E 

02 

5.015*55-03 

7.900ZZE-O2 

209055 

00 

190 

—5* l00l2e-02 

S.10657E-01 

4,336505 

01 

2.971765-04 

5 

•14420E 

02 

—5 .01 3555—03 

7*960226-02 

209055 

00 

191 

— 1 *S220l£— 01 

7. 978 735-01 

4.33650E 

01 

2.971765-04 

S 

.144205 

02 

-1*496745-02 

7.65489E-02 

209055 

00 

192 

— 2 *51001E— 01 

7.72505F-01 

4*338605 

01 

2.971 765-04 

5 

•14420E 

02 

-2.468495-02 

7.606125-02 

20905E 

00 

193 

-3*458435-01 

7.3495SE-01 

4.33850F 

01 

2.971765-04 

5 

.144205 

02 

— 3.401 9ZC—Q2 

7.  23773E-02 

209055 

00 

194 

-4.3523OE-01 

6*858145-01 

4. 338505 

01 

2.971 76E-04 

5 

.144205 

02 

-4.281175-02 

6. 755325-01 

20905E 

DO 

195 

“5*1 779 SC— 01 

6.25658E-01 

4.338S0E 

01 

2.9T176E-04 

5 

* 1 4420E 

02 

-5.Q933AC-02 

6*166315-92 

2090 5E 

00 

196 

—9*921 1 IE— 01 

5.5603  IE— 01 

4.33850E 

01 

2.971 76E-04 

5 

.144205 

02 

—5.825865—02 

5*479815-02 

209055 

00 

197 

-6*571315-01 

4 .774365-01 

4.33850E 

01 

Z. 971 765-04 

3 

•14420E 

02 

-6.466445-02 

4.706475-02 

209055 

00 

196 

-7*117885-01 

3*913105-01 

4*33  06 OE 

01 

2.9  71765-04 

5 

.144206 

02 

-7*005146-02 

3.86840E-02 

209055 

00 

199 

— 7»55220E-O1 

. 2.990 14E-01 

4.338505 

01 

2*971 76E-04 

S 

,144206 

02 

-7.4333*5-02 

2.9469EE— 02 

209055 

00 

200 

-7*867415-01 

2.02002E- 01 

4.336S0E 

01 

2*971 76C-04 

5 

.144205 

02 

-7,744206-02 

1.992515-02 

209055 

00 

BOOT 


9 


89-8£-HJ.-Da3V 


MACH  4 NOZZLE  (G-l.Z*) 


ZE  * 2. 05000E-01  STEP 


NO* 

X 

V 

P 

MHO 

201 

-9 .038556-01 

1.01804E-01 

4. 

338606 

01 

2.971 75E -04 

*02 

-8.1230OE-QI 

1 . 284676—06 

4. 

33850E 

01 

2.9T178E-04 

203 

-0.05ft 556- 01 

- 1.01B02E-01 

4. 

338506 

01 

2.971  TOE -04 

204 

—7.867426—0 1 

-2.020006-01 

4. 

330506 

01 

2,971706-04 

205 

-7, 3322  IE—  01 

-2.99012E-01 

4. 

33850E 

01 

2.97179E-04 

206 

—7 • 1 17902-01 

— 3.9 13 DTE— 01 

4. 

338ft  OE 

01 

2.971706-04 

20T 

-6.571336-01 

-4.77433E-01 

4. 

3305OC 

01 

2.971 766-04 

20ft 

-5.921 I3E-01 

-5.360296-01 

4. 

338506 

01 

2.971706-04 

209 

-5. 177559—01 

-6 .250566- 01 

4. 

330596 

01 

2.9TJT6E-04 

210 

-4.3S23ZE-01 

-6.958136-01 

4. 

338506 

01 

2*971 762-04 

211 

-3.45845E-01 

-7.149S4E-01 

4. 

330506 

01 

2.97176E-04 

212 

-2.510046-01 

— 7 .7250SC— 01 

4. 

338506 

01 

2.971706-04 

213 

-1.5**056-01 

- 7.97BT2E— 01 

4. 

330906 

01 

2.971 76E-04 

214 

-5.10033E-02 

-8. 106376-01 

4. 

338506 

01 

2.971 70E-O4 

215 

S .099836- 02 

-ft. 106576-01 

4, 

32  85  BE 

01 

2.971766-04 

216 

1.921906-01 

-7,970 736 -01 

4. 

338506 

01 

2.971 706-04 

21  * 

2.50997E-9I 

-7.  72507C-01 

4. 

3388 OE 

01 

2.971 766-04 

218 

3.4S941E-01 

-7.349566-01 

4. 

338506 

Ol 

2.471766-04 

219 

4.952206-01 

-6.650156-01 

4. 

338506 

01 

2.971 76E-04 

2*0 

5. 177516—01 

-6.238596-01 

4. 

330506 

01 

2.971766-04 

2*1 

5 .921096—01 

-S .560336- 01 

4. 

33  8506 

01 

2.971766-04 

222 

6 .571 296—01 

-4.774396  - 01 

4. 

330506 

01 

2.971 766-04 

223 

7.  117076-01 

- 3.91 3126-01  - 

4. 

330506 

01 

2.47176E-04 

224 

7.552196-01 

-2.990166-01 

4. 

330506 

01 

2.971 76E -04 

2*5 

7.86  7406-01 

-2.O2OOSE-01 

4. 

330506 

01 

2.971 706-04 

226 

8.0SB546-01 

-i.oieooc-Di 

4. 

338506 

01 

2.971706-04 

22T 

9i  13T90E— 01 

0.0 

3. 

924B06 

01 

2.74O0TE-O4 

22ft 

9.06244E-01 

1 . 005246-01 

3. 

924S0E 

01 

2. 740876-04 

2*9 

ft  v 91673E— Ol 

1.99820E-01 

3. 

92450E 

01 

2 • 740576-04 

230 

8.642766-01 

2.9670TE— 01 

3. 

924506 

Ol 

2.74087E-04 

231 

6.26393E-01 

3.B99846-  01 

3. 

924S0E 

01 

2.740076-04 

232 

7.784766—0 I 

4. 7B527E— 01 

3. 

924506 

01 

2«  74087E-04 

233 

7 .211096-01 

3.61 261ft— 01 

3. 

924506 

01 

2. T40ft 7£— 04 

234 

6.54999E-01 

6.371026-01 

3. 

92450E 

01 

2. 74087E-Q4 

23ft 

5 .009 1BE— 0 1 

7.053696-01 

3. 

924596 

01 

2. 740076-04 

236 

4.99796E-01 

7.049946-01 

a. 

924500. 

01 

2.740876-04 

2 37 

4.12607E-01 

0.  1 93326-01 

3. 

92450E 

01 

2,740076-04 

23ft 

3.204106-01 

6.  3577*6—01 

3. 

924506 

Ol 

2.740076-04 

239 

6.20236-01 

ft. 65S29E-01 

1* 

924506 

01 

2.740876-04 

240 

1.2SB1 3E-01 

9.05129E-01 

3. 

924506 

01 

2.74087E-04 

2*1 

2.SL796E-02 

9. 134436-01 

3. 

924S0E 

01 

2.740076-04 

242 

-7.S4595E-02 

9.10669E-01 

3. 

92450E 

01 

2. 74Q87E-04 

*43 

-U7B183E-01 

ft  .96841 E— 01 

3. 

924S0E 

01 

2. 740876 -04 

2«4 

-2 .727786-01 

B.72126E-01 

3. 

924S0E 

01 

2. 740476 -04 

245 

-3.670646-01 

8. 36 9256- Ot 

1. 

924506 

01 

2. 740876-04 

246 

— 4 . 568946-0 1 

7.91366E— OK 

3. 

924506 

01 

2.74O87E-04 

247 

—5.411 7BE— 01 

7.3630 IE— 01 

3. 

924506 

01 

2. 740076-04 

848 

—6. 1BB92E— 01 

6. 722986—01 

3. 

924806 

01 

2.740876-04 

249 

— 6. 89 094 E— 01 

6.00134E— 01 

3. 

9245 OE 

01 

2. 740876-04 

250 

-T. 509116-01 

3 . 2 06  66  E-  01 

3. 

924B06 

01 

2. 740076-04 

AEDC-TR-78-6B 


MUCH  4 NOZZ1_e  (G«tt24i 


Zfi  • 2*053097-01  STEP  NO.  9 


NO. 

2 

Y 

P 

RHO 

O 

THETA 

P5I 

N 

bgOy 

2SI 

— 9*035 53E— 01 

4.349176-01 

3.924506 

01 

2.74Q07E-O4 

5. 416758 

02 

-9,704908-02 

S.ZT109E-0Z 

1 .2060 OE 

60 

252 

-6.466205-01 

3. *30696-01 

3.92450F 

01 

2, 740678-04 

5. 4L  0758 

02 

-1*022566-01 

4.1 64 DEE— 02 

1,266008 

00 

253 

-0.T931OE-OI 

2. *06*78-01 

3.47450E 

01 

2.74O87E-04 

5.410TS6 

OZ 

— I .062198—01 

3.015428-02 

1.266008 

00 

254 

— 9.01327E— 01 

1. 50406E— 01 

3.424508 

01 

2, 740B7E— 04 

S.41B7SE 

02 

-1 .006096-01 

|,82436E-02 

1.286008 

00 

255 

—9. 1 2402E— 01 

5.03401E-02 

3. 92450* 

01 

7,740878-0* 

5.419758 

OZ 

-t. 102336-01 

6.106676-03 

1 ,286  008 

00 

254 

-9,1  ? *037-0  1 

-5.03372E-02 

3.9Z450E 

01 

Z.74087E-Q4 

5.4167SE 

02 

-I.I0233E-01 

—0.1 06328—03 

1 .206008 

00 

257 

—9.01 32  7E— 01 

-1. 5040 36-01 

3.924308 

01 

2. 74007  E— 04 

5,41 6756 

02 

-1. 000098-01 

-1,0243*8-02 

1.Z06DO8 

00 

256 

-9. 793112-01 

-2i  466446-  01 

3.92450E 

01 

2. 740876-04 

5*418758 

OZ 

-1 *062 | 9 E— 01 

— 3.01B30E-92 

1*266008 

O0 

20* 

-0.46621E-O1 

-3. 4306 76-01 

3.924 50C 

01 

Z.74O0TE-O4 

5,41 07SE 

02 

— 1 .Q2256C— 01 

-4 .169026-08 

I .266008 

00 

260 

-8.03655E-01 

—4,3491 4E— 61 

3.92450E 

01 

2.74Q97E-04 

5.410788 

02 

-9. 705006-02 

-5.271  058-02 

1.266008 

00 

261 

— 7*  50 9336- 01 

-3*206038- 01 

3.924508 

01 

2.740678-04 

9.41675C 

02 

-9 .066528-02 

-6.306028-02 

1.26600E 

00 

262 

-6.990962-01 

-C.00132E-91 

3.92*508 

01 

2. T4O07E-O4 

9.418798 

02 

-0.3101 IE-Q2 

-7.2673 7E-02 

1 *266008 

00 

263 

-6.100948-01 

-6. 722966-01 

9.92*506 

01 

2*  74Q07E— 04 

5.41875E 

02 

— 7.4090 JE— 02 

-6. 13t60e-02 

1.206008 

00 

264 

-5.41190E-01 

-7.36299E-01 

3*924608 

01 

2.74097E-04 

5.41 8798 

02 

-6.524718-02 

- 6.9O039E— 02 

1.266008 

00 

265 

-*.560968-01 

-7.913646-01 

3.9245DE 

01 

2.740078-04 

5.41075E 

02 

-5.51165E-02 

-9. 570 706-02 

1 .206008 

00 

266 

-3.67067E-01 

-6*368238-01 

3.924508 

01 

2.74007E-04 

5.4 187BE 

02 

—4*427238—02 

-1.01 I60E-O1 

1.296008 

00 

26T 

-2. 727028-01 

-0.721 25E-01 

3. 92450E 

01 

2. 74O07E-O4 

5.4107SE 

02 

-3,209568-02 

-1.054058*01 

1 .206 OOC 

00 

266 

-I.75I96E-01 

— 6, 9694 QE- 01 

3.92450E 

01 

Z.74007E-04 

5.41075E 

02 

-2.112408-02 

-1.005648-41 

1,206008 

00 

269 

-7*546588-02 

-9.106686-01 

3.92450E 

01 

Z.74057E-04 

5.410T5E 

OZ 

-9.099I4E-O3 

-1*100278-01 

1.206008 

00 

Z70 

2*51 7598-02 

-9.134426-01 

3 . 9Z450E 

01 

Z.T4067E-04 

5.4107SE 

02 

3,035496-03 

-1.103B9E-01 

1.206008 

00 

271 

1 .25507E-01 

- 9.  051 306-01 

3.92450E 

01 

2.74007E-44 

S.41075E 

02 

1.51 331 E— 02 

-1.49363E-01 

1. 2B600E 

00 

272 

2-2*3197-01 

— B.05029E— 01 

3.924908 

01 

2.74067E-04 

S.4167SE 

02 

2.704908-02 

-1. 070498-01 

1.2S600E 

00 

273 

3.20404E-01 

- 0. 557706—  01 

3.92450E 

01 

2.74O07E-Q4 

5.410758 

02 

3,664128-02 

-1. 034446-01 

1.206068 

00 

274 

4,1 76048-0* 

-0*153348-01 

3.924506 

01 

2.74OB7E-04 

5.418758 

02 

4.976B6E-02 

—9.850748—02 

1 ,286008 

oo 

Z7S 

4*9979  IE— 01 

-7.649976-01 

3.924508 

0! 

2.74087E-04 

5.410758 

02 

6.029708-02 

-9.293666-0 2 

1.206008 

00 

276 

3.009168— Oi 

—7. 0537 IE— 01 

3 .92450E 

01 

2.740978-04 

5.41675E 

02 

T.00491E-02 

— 6 » 53605E- OZ 

1 .206008 

00 

277 

6.5498SE-0I 

-6.371  068-01 

3.9**508 

01 

2, 74 09 TC -04 

6. 410758 

02 

T.90547E-O2 

-7.714358-01 

1.28600E 

oo 

279 

7.21107E-01 

-5.612638-01 

3.924506 

01 

2. 74 6fl 78-04 

S.4I675E 

02 

0.705478-02 

-6*7961 66-02 

1.206OOE 

00 

279 

7,784736-01 

-4.76531E-01 

3. 924506 

01 

2.740078-04 

5.416758 

02 

9 ,399966—02 

-5.798566-02 

1 .206008 

00 

290 

9 , 263918—01 

-3.099868-01 

3.924506 

01 

2.740978-04 

5,41 6758 

02 

9.900478-02 

-4.T2739E-Q2 

1.20600E 

00 

ZB  1 

6.64277E-01 

-2.967128-01 

3.9245DE 

01 

2.740678-04 

5.41675E 

02 

1 .0*3968-01 

-3.507056-02 

1* 296008 

oo 

262 

9. 916725-01 

-1.99835E-01 

3.9Z450E 

01 

Z.74007E-04 

5.41075E 

02 

1.077168-01 

— 2 . 42  37  2 E-  02 

1.206008 

00 

203 

9.O0243E-O1 

-I.005Z9E-01 

3.92450E 

ot 

2 • 74 0 87 E— 04 

5.410758 

02 

1.097208-01 

—1 .219468—02 

2 .266008 

00 

264 

1.01532E  60 

0.0 

3.39530E 

01 

Z. 43670 E-0« 

5.703662 

OZ 

1.914028-01 

0.0 

1*392038 

00 

YES 

265 

1.01027E  00 

U01O93E-0S 

3.395308 

01 

2.4307DE-04 

5.7A300E 

02 

1.5O703E-01 

1.51059C-0Z 

1 .392036 

00 

YES 

266 

9,451996-01 

2.01182**01 

3*  3953 QE 

01 

2.49870E-04 

5.703068 

02 

1 .40  *366-01 

3.02340E-02 

1.392038 

00 

yes 

207 

9. 702I2E— 0 1 

2.99271E-01 

3.395308 

01 

2.430T0E— 04 

5.70300E 

02 

1.446848-01 

4.44604E-02 

1.392038 

00 

ves 

266 

9, 353938— Oi 

3.94365E-01 

3.39530E 

01 

2.436708-04 

5.703068 

0* 

1 .3940 TE— 01 

5.921 768-02 

1*392038 

00 

▼8S 

209 

5 *9 16768- 01 

4.05980E-01 

3 .395308 

01 

2. 430 TOE— 04 

6.7B30BE 

02 

1.329006-01 

7.206708-02 

1.392036 

00 

YES 

290 

6.36B97E— 01 

5.71950E-01 

3.39S30E 

01 

2. 430708-04 

S.7B3BBC 

02 

1.249918-01 

0.576918-02 

1.392038 

00 

YES 

291 

T.T7760E-01 

0.5263SE-01 

3. 3953 OE 

01 

2. 438 TOE— 04 

5.783008 

02 

1,156428-0} 

9.779028-02 

1 ,392038 

00 

YE5 

292 

7 . 06934E- 01 

7. 260338-01 

3.395308 

01 

2,439708-04 

5.78306E 

02 

t, 055*68-01 

1.O0032E-Q1 

1 * 39  20  3 E 

00 

▼ 68 

*93 

6.33442E— 01 

7.93809E-01 

3.39S30E 

01 

2* 43670 E— 04 

5 .763568 

02 

9*421376-02 

L.I07TOE-01 

1.392038 

00 

YES 

294 

5. 50 05 BE— 01 

6.52095E-01 

3.395308 

01 

2. 43870 E— 04 

5.763608 

02 

6.195308-02 

1.275108-01 

1. 39203 E 

00 

YES 

245 

4.63201 E— 01 

9.03S04E-D1 

3, 3953  08 

01 

8.438708*04 

5.703068 

02 

6.086926-02 

1.349968-01 

1 .39* 03E 

00 

YES 

296 

3. 70939E— 0 1 

4.451346-01 

3.395308 

ot 

2.4307OE-O4 

5 .763008 

02 

5*515208-02 

1**1 1398-01 

1.392038 

00 

YES 

297 

2.74991E-01 

9.773T1E-01 

3,395308 

01 

Z. 430  70  E— 04 

5.763B6E 

OZ 

4. 0S7 708- 02 

1.450628-01 

1.392038 

oo 

YES 

298 

t .763IOE-OI 

9.99895E-01 

3.3953  Oe 

01 

Z.4307OE-O4 

9.78300E 

02 

- 2.620306-02 

1.491958-01 

1.392038 

00 

YES 

299 

T.56T59E-02 

1.01Z4BE  00 

3.39S30E 

01 

2.4307OE-O4 

S.70360E 

02 

1.127598-02 

1 .510456—01 

1,392038 

00 

YES 

300 

-2.531092-02 

l.OtSOOE  00 

3.395308 

01 

2.4 30 70E-04 

5.76368E 

02 

— 3.T0134C— 03 

l-i  51  4168  — 01 

1 ■ 392038 

00 

YES 

AEDC-T  R-78-68 


MACH  4 NOZZLE  CG*t*24| 


ZG  - Z*Q5009F-Qi  STEP  MO.  0 


Nd. 

X 

T 

P 

RHO 

Q 

7HETA 

P#I 

N 

0007 

301 

-1,282476-01 

1 *0074  46  OO 

1.395306 

01 

Z *436706-04 

5*783586 

02 

-1*570236-02 

1.503046-01 

I.39203E 

03 

ves 

302 

-Z*ZS9Z9E-0l 

9*098646-01 

3*395306 

01 

2*438706-04 

8*783006 

02 

-3*350096-02 

1.477206-01 

1 .792036 

00 

VE5 

303 

-3.2J3A4E-0I 

9*624506—01 

3*395306 

01 

2.43870E-O4 

5*703606 

02 

-4*007286-02 

1.43886E-01 

1.39203E 

00 

YES 

304 

-4.17566E-0I 

9*254  7IE-01 

3*195306 

01 

2.43B7QE-04 

5.703056 

02 

-6*209616-02 

1.3523TE-01 

1.392036 

00 

VES 

39  5 

-S.076S9E-01 

8*792936-01 

3*395306 

01 

2*430706-04 

3.703006 

02 

—7 *551326—02 

1*314206-01 

1*392036 

00 

YES 

300 

«S,92UfiE-0l 

8*243776-01 

3*  39530E 

01 

2,438706-04 

5 *763006 

02 

-0*019146*02 

1 *232986—01 

l *392036 

00 

YES 

SOT 

-4*718216-01 

7.612696-01 

3*393106 

01 

2.4307OE-O4 

5.783006 

oz 

-I *000046-01 

1*139446-01 

1*392036 

00 

YES 

300 

-7*442016-01 

6.905946-01 

3*395306 

01 

2 . 43870  6-04 

5*703006 

02 

-1*108326-01 

1*034456-01 

1*392036 

00 

YES 

309 

— 8*093446-01 

6.I3O56C-01 

3*3953 OE 

01 

2.43570E-04 

5*703006 

02 

— 1 .205066— 0 1 

0*189986-02 

1*392036 

00 

YES 

3 10 

-8.663626-01 

5. 294206-01 

. 3 * 3953  OE 

01 

2 .43070E-04 

5.703006 

02 

-1.29100E-01 

7.941  996  - 02 

1*392036 

00 

YES 

31  1 

-9,147716-01 

4. 4 053 26- 01 

3*3953 OE 

01 

2. 43870 E -04 

5*783806 

02 

-1,383636-01 

6 ■ A 1 2776—  02 

1 *392036 

00 

ves 

312 

-9.540BBE-01 

3*4 726 IE— 01 

3*  3953 OE 

01 

2*438706-04 

5*753806 

02 

— 1 *422636—01 

8.215976-02 

1.39203E 

00 

YES 

313 

-9.83923E-01 

2*505406-01 

3*3953  06 

01 

2.430706-04 

5*703006 

02 

-1 .467446-01 

3*704536-02 

1*392036 

00 

Y63 

314 

-1.00398E  00 

1*513276-01 

3*395306 

01 

2*430706-04 

3*703006 

02 

-1.497876-01 

2*274406-02 

1 *392036 

00 

YES 

310 

-U01406E  00 

5*061126-02 

3*  3953 OE 

OK 

2*439706-04 

9*783806 

02 

-1*512726-01 

7.808132-03 

1*392036 

00 

VES 

316 

-1.01406E  00 

-5.960706-02 

3*3953 OE 

01 

2. 4ia roe-04 

9.783006 

02 

-1*51 272E-01 

— 7.6075OE— 03 

1 .292036 

00 

YES 

317 

-1.O039BE  00 

- 1 .513236-01 

3*  39530E 

01 

2. 438706-04 

5.T03B0E 

02 

-1.497576-01 

-2*2  74426-02 

1*392036 

00 

YES 

318 

— 9*  93 92 4 E— 01 

-2*505366-01 

3*395306 

01 

2*438 706—04 

9.783086 

02 

-1.467446-01 

-3*784476-02 

1.39203E 

00 

YES 

319 

—9*540096—01 

-3*472576-01 

3.395306 

01 

2.43870E-04 

5*783806 

02 

-1,422646-01 

-5*216516-02 

1*392036 

00 

YCS 

320 

-9.1 47736—01 

-4.49S29E-01 

3*395306 

01 

2 . 438 TOE-54 

5*783006 

02 

-1.36364E-0I 

-0.01272E-O2 

1.392036 

00 

YES 

321 

-0*663646-01 

— 5*  Z9422E— 01 

3*3953 Oe 

01 

2*430706-04 

9*783006 

02 

-1,291066-01 

-7.941946-02 

1*392036 

00 

YES 

322 

-8.093466-01 

— 6.  130526-01 

3.395306 

01 

2.438706—04 

5*783866 

02 

-1*205666-01 

-9*189936-02 

Z .392036 

00 

YES 

323 

— 7. 44284 E— 01 

-A. 9059 IE- 01 

1*395306 

01 

2. 4J8 fOe -04 

5.78 3886 

02 

-1*108326-01 

-I.09445E-01 

1 .39201E 

00 

ves 

324 

-6. 71 824£— 01 

- 7*612666-01 

3*395306 

01 

2*4  38706-04 

5*783856 

02 

-1*000046-01 

- 1*139446-01 

1*392036 

00 

YES 

32S 

-5.9260BE-O1 

-8. 243766-01 

3.3953  OC 

01 

2*438706-04 

5.703086 

02 

-8*019176-02 

-1*232986-01 

1 *392016 

00 

ves 

324 

-5.076626-01 

-«*T9291E-01 

3.39530E 

01 

2.438706-04 

5 .TBJSBE 

02 

-7*851 306-02 

-1  .3142OE-01 

1*392036 

00 

VES 

3*7 

-4.175906-01 

-9.254696-Q1 

3 .3933  OE 

01 

2.43B706-04 

5*783886 

02 

-5*209076-02 

-1*382386-01 

1 *392036 

00 

YES 

328 

- 3*2336 8E-01 

-9*624  496-  01 

3.3QS30E 

01 

2 *43  8 TOE -04 

5.70380E 

02 

-4.007326-02 

-t*43O0«E-Ol 

! .392036 

00 

ves 

329 

-2*259366-01 

-9.89862E— 01 

3.3953 OE 

01 

2 * 43870 e -04 

5. 703686 

02 

-3*350206-02 

-1.47720E-01 

1.392036 

00 

V6S 

330 

-1 *262546-01 

-1.007446  00 

3.395306 

OK 

2.438T0E-P* 

5*703006 

02 

-1*076336-02 

- 1. 50 3 04 E- 01 

1 *392036 

00 

YES 

331 

-2*531706-02 

- 1 .015006  00 

3*39 S3 OF 

01 

2. 43B70E— 04 

5*703006 

02 

— 3 . 76229E-03 

-1.51 4106-01 

1 *392036 

00 

YES 

332 

7*587176-02 

-I.01246E  00 

3,395306 

01 

2.43870E-94 

5.70380E 

OZ 

1.127S3E-02 

- I .810456-01 

1*392036 

00 

YES 

333 

1 *763026-01 

-9*998966-01 

3.39530E 

01 

2.430706-04 

5. 703006 

02 

2.62O26E-02 

-1.491956-01 

1*392036 

00 

VES 

334 

2*749046-01 

-9.77373E-01 

3*395306 

01 

2.4357QE-04 

5*783886 

02 

4*057806-02 

-1.450036-01 

1*392016 

00 

YES 

335 

3*709146-01 

-9*451  366-01 

3.395306 

01 

2.43B706-04 

5*783886 

02 

5.51S1 96-02 

-1*411306-01 

1 « 3920 36 

00 

TCS 

336 

4.63193E-01 

-9*035086-01 

3.395306 

01 

2.435708-04 

5*703886 

02 

8. 88 88 16-0 Z 

— 1 • 3499AE-0 1 

1 . 392036 

00 

VES 

337 

5 .50B53E— 0 1 

-8.528908-01 

3*395306 

01 

2.43B70E-04 

5.783556 

02 

B *195226-02 

— t .275196-01 

1 .39201E 

00 

YES 

330 

6. 330306-01 

-7*938126-01 

3.39530E 

01 

2*438706-04 

5*703606 

02 

9.421316-02 

-1.187706-01 

1*392036 

00 

YES 

339 

7*0893 16-0 1 

-7*  268366-01 

3.395306 

01 

2*438706-04 

5.7B30B6 

02 

1 .055406-01 

-1.OO032E-O1 

1,392036 

00 

VES 

340 

T*T777ee-01 

— 6*  52640c- 01 

3*395306 

01 

2.438T06-O4 

S .789806 

02 

1 • 158428—01 

-9*779706-02 

1*392036 

00 

YES 

341 

8*386936- 01 

“5*  7195SE-01 

3.395306 

01 

2.43870E-04 

5*7030 0E 

02 

1.249906-01 

— 0.57A98E— 02 

1*392036 

00 

YES 

342 

0 .9 16T4E— 0 I 

-4. 005  856-  01 

3.39530E 

01 

2.438706-04 

5 . 7830  BE 

02 

1*329006-01 

— 7*  20A7AE—02 

1 *392036 

00 

ves 

343 

9*355906-01 

-3.943926-01 

9*395306 

01 

2*430706-04 

5.7830 6E 

02 

1 .394076—01 

-5.921506-02 

1*392036 

00 

YES 

344 

9*702106-01 

-2.99276E-01 

3.3953 OE 

01 

2*435706-04 

5.70308E 

02 

1.44C04E-01 

-4*495936-02 

1 *392036 

DO 

YES 

345 

9 . 95 1 08E— 0 1 

711866-01 

3*395306 

11 

2.43070E-04 

5.703086 

OZ 

1.40436E-OI 

-3*023476-02 

1*393036 

00 

YES 

346 

1*010276  00 

-1*0 10976-01 

3.395306 

01 

2,498706—04 

5.783806 

02 

1*507036-01 

—1 .619656-02 

1 *392036 

00 

ves 

AEDC-TH-78-6S 


MACH  4 NOZZLE  (G»1.24) 

BODY  POINTS  AT  ZE  ® 2.05090E-Q1  STEP  NO.  0 


MO. 

X 

Y 

P 

RHO 

0 

THETA 

P5I 

K 

204 

1 .015326  00 

0.0 

3. 3953 OE 

01 

2 » 430TDE— 04 

5,7«3aaF 

02 

1.S1462E-01 

0.0 

1 

.392036 

00 

2BB 

I.0I027E  00 

1.01093E-01 

3.39530E 

01 

2*436702— Q4 

S.7B38SE 

02 

1.5O703E-01 

1.519596-02 

1 .392036 

00 

206 

9.951696-01 

?•  0 1 182E-01 

3.3953 OE 

01 

2.43B70E-Q4 

3.783686 

02 

1.46430E-01 

3.0234QE-02 

.392036 

00 

267 

9.7Q21ZE-01 

7.992 TIE- 01 

3.3933 OE 

Ql 

2 * 436  TOC -04 

5.783B6E 

62 

1. 44684E-01 

4.49584E-92 

* 39203E 

00 

208 

9.355936-01 

3.943 856-01 

3.395302 

01 

2.43870E-04 

3.783806 

02 

1.394676-01 

5.921 766—02 

*192031 

00 

269 

6,9)6766-01 

4.05S60E-01 

3 .393302 

01 

2.436T0G— 04 

5 >783086 

02 

1.329006-01 

7. 266706— 62 

.392036 

00 

290 

5.306976-01 

5.71950E-01 

3. 3953 OE 

01 

2, 4 30 TOE— 04 

5.78J66E 

02 

1.2499  IE-01 

8.5769  IE-02 

. 39293E 

00 

291 

7,777806-01 

6.52635E-01 

3. 29330E 

01 

2. 43070 2-0* 

S.76388C 

02 

1.I&842E-01 

9.779626-02 

.392036 

00 

292 

T.08934E-0I 

7. 26833E— 01 

9.3953 OE 

01 

2.436T0E-04 

5.T6386E 

02 

1.05B46E-0I 

1.08632E-01 

•39203E 

00 

293 

6.31042E-O1 

7.93s  096 -01 

3. 3953 OE 

01 

2 .4  38 70 E— 04 

5.75380E 

02 

9,421376-02 

1.1B770E— 91 

,392036 

00 

294 

5.S06S8E-91 

6.S289SE-01 

3.395302 

01 

2.436706-04 

5 .763666 

02 

8. I9530E-02 

1*275166-01 

.39203E 

00 

290 

4.632012*01 

9.03504E-01 

3. 3953  02 

ol 

2.436702-04 

5.763606 

02 

6.096926-02 

1 .249966— Ol 

.392036 

oo 

296 

3. 709 392— 01 

9,451  342-  01 

3.395306 

01 

2.438706-04 

5.783866 

02 

. 5.5I520E-02 

1.41 L35F— 0 | 

• 39293E 

00 

297 

2.74991E-01 

9.773T1E-01 

3.395302 

01 

2,438  70  E— 04 

5 .763666 

02 

4. Oar TOE- 02 

1 .45BB2E-61 

• 392 03E 

no 

296 

1.763102-01 

9.99095E-O1 

3,395306 

01 

2.438706—04 

S.7630SE 

02 

2,620386-02 

1.491956-0 1 

.392036 

00 

299 

7.567S9E— 02 

1. 012*86  00 

3.395306 

01 

2.43870E-04 

■ 5.783002 

02 

1, 127596-02 

1 .510*56-0 1 

•392036 

00 

366 

-2.53109E-OZ 

I.OI500E  00 

3.3953 OE 

01 

2. 436  TOC— 04 

5.7B360E 

02 

-3.761 39E-03 

I.S1416E-0! 

.392036 

OQ 

30  1 

-1.262472-01 

1 .007442  00 

3.395302 

01 

2.436T0E-04 

5.703802 

02 

— 1 .876236—02 

1.563 04E-01 

.39203E 

00 

302 

-2. 2592 92-01 

9.64664E-01 

3.39530E 

01 

2.430702-04 

5. 7038  BE 

02 

-3.346096-08 

1.477206-01 

*392336 

09 

303 

- 3.23364E-01 

9.624 506-01 

3.39S30E 

01 

2.430702-04 

5,703666 

02 

-4,607266-02 

L .436066-01 

.392036 

00 

304 

-4. 175062-01 

9.25471E-01 

3. 39 5302 

01 

2.4307OE-O4 

5.763002 

02 

-0,2096  IE-02 

I.3973TE— 01 

.392036 

00 

306 

-5.07659E— 01 

8.792932-01 

3.39530E 

01 

2.4307OE-O4 

5.7638 BE 

02 

-7.551  32  E—  02 

1.3142OE-01 

.39Z03E 

90 

306 

-S.9206SE-91 

8,2437 TE-Ol 

9.  3953 OE 

01 

2.43670E-04 

S .703006 

02 

-6.819 I4E— 02 

1 « 23298E— 01 

1.392036 

00 

307 

-6. 716216-01 

7. 61269E— 01 

3.39S30E 

Ol 

2.430702-04 

5.7B368E 

02 

— 1 *000046-0 1 

I.  I3944E-01 

,392036 

00 

306 

—7. 94 26 IE— 01 

6.905942-01 

3.  3933  02 

01 

2.43670C-04 

5.70300E 

02 

— 1 .10832E—01 

1. 034456-01 

.392036 

00 

309 

“8.09344p— 0 1 

6,1 3066E- 01 

3.395306 

01 

2.43670E-04 

5.763062 

02 

* 1 ,205666—01 

9.  18998E-0Z 

•39Z03E 

00 

310 

-6. 663 6 2 F— 01 

5.29425E-0I 

3. 3953 OE 

01 

2 .43670 2-04 

9.783802 

02 

- 1.291 066*01 

7.941 99C-02 

.392036 

00 

311 

*9 .147712—01 

4,4053  2E-01 

3,3953  02 

01 

2.43670E-04 

5.70366E 

02 

-1.36S63E— Ot 

6.612T7E-02 

•3920 36 

oo 

312 

—9 . 540886- 6 1 

3.47261E-02 

3.39B30E 

01 

2,43670E-Q4 

5.793662 

02 

-1.42263E— 01 

5.2156 7E-02 

•392 03C 

00 

313 

— 9 .63923 E— 01 

2.50540E-01 

3. 395306 

01 

2.43B70E-04 

B,7B386E 

02 

-1 .467446-01 

3 ,764536—02 

• 39203E 

00 

314 

-1.003986  00 

1 .5132  7E -01 

3,3  953  OE 

01 

2 » 438 TOE— 04 

S.7S306E 

Of 

*1.497576-01 

2,2  74486-02 

.392 036 

00 

319 

— 1*014066  09 

5.061126*02 

3.39530E 

01 

2.4  38  TOE— 04 

9.76386E 

02 

-1.91272E-01 

7*6061 76—03 

*392036 

00 

316 

- 1 . 01 4060  00 

- 5. 060702—02 

3, 395302 

01 

2. 4 3670 E— 04 

S.76368E 

02 

-1.SI272E-01 

—7*607506—03 

.392  036 

00 

317 

-1.003966  00 

-1.513232-01 

3.395302 

01 

2. 436705-04 

5.78366C 

02 

- 1.497576-01 

-2.27442E-02 

.392036 

00 

316 

-9-63924E-01 

-2.505366-01 

9*  39 530 E 

01 

2.436702-04 

5.78306E 

02 

-L,44T*4E-01 

-3.764476*02 

.392036 

00 

319 

-Q.&40B9E-0I 

-3.479576-01 

3. 39*5306 

01 

2.436702-04 

5.T6306E 

02 

-1.42264E-01 

-9.2155  16*02 

.392036 

00 

320 

—9 . 147T3E— 01 

— 4 .40529C— 01 

3.3953DE 

01 

2. 438 TOE -04 

5.76386E 

02 

-I « 36  364E- 01 

-6,612726-02 

,392036 

00 

321 

-6.66364E-01 

-3.294222-01 

1. 395 302 

01 

2.43S70C-94 

5.783686 

02 

* 1. 291006—01 

-T.941 94E-02 

. 39203E 

00 

322 

-8.09346E-6L 

-6.130522-01 

3.3953 OE 

01 

Z. *30706— 04 

5.78360E 

02 

-1.20S06E-01 

-9.1 0993E— 02 

• 3 92  03E 

00 

323 

—7 . 442846—01 

-6.90391E-01 

3.395306 

01 

2.430702-04 

S. 763862 

02 

-1.I0S32E-OI 

-1 .034466-01 

.39203E 

00 

324 

-6.716242-01 

-7.61266E-01 

3.39530^ 

91 

2.  *30 TOE— 04 

5. 78 38 BE 

02 

-1.000D4E-OI 

-1.13944E-01 

.392036 

00 

32S 

-9. 920866-01 

-6.243762-01 

3.39S30E 

01 

2.438706-04’ 

5.7B386E 

02 

-8.819176-02 

- I .232986-01 

. 39203E 

oo 

326 

—5*076626-0! 

-0.792916-01 

3.395302 

01 

2,430706—04 

5 a 70306E 

02 

—7.551 3 BE— 02 

-1«31*20E-01 

.392036 

00 

327 

-4.175906-01 

-9.25469E— 01 

J.9953QE 

01 

2.43670E-04 

S.763B6E 

02 

-6.20967E-02 

-1.382366-01 

.39203E 

00 

326 

-3.233686-01 

-9.624492-^Ql 

3.39530E 

01 

2. 43 8 TOE— 04 

5.78380E 

02 

—4.607326—02 

-1. 436066-01 

•39203E 

00 

329 

-2.25936E-01 

— 9 • 89862 E- 0 1 

3.39530E 

01 

2. 438706-04 

5.7830 BE 

02 

-3. 356206-02 

-1  .477806-01 

. 392 01 E 

00 

330 

-1  .26  2 546-  01 

-1.00740C  00 

3.395302 

01 

2. 438 7QE-64 

5.76386E 

02 

-1.876336-02 

*1.503046-01 

.392036 

00 

331 

-2.531702  - 02 

-1.01300E  00 

3 .39530E 

01 

2,430706-04 

5*7638 8E 

02 

-3.762396-03 

-1. 814166-01 

.392036 

00 

332 

7.667172-02 

-1.01248E  00 

3.3953 OE 

.01 

2. 430 TOE-04 

S. 783882 

02 

1 .127036-02 

-1  .51 0456-  01 

.392  03E 

00 

333 

1 . 76302E-01 

-9.998962-01 

3.395302 

01 

2.4307OE— 04 

5.763662 

02 

2.020266-02 

— 1 *491926-01 

.39Z03E 

00 

334 

2.74Q64C-01 

-9.773736-01 

3.395306 

01 

2.4  38 70E— 04 

5,783686 

02 

4. 08 76 OE- 02 

-1 .458036-01 

.392036 

00 

335 

3.709342-01 

-9,451362-01 

3*3  953  OE 

01 

2.43870E-04 

5. 783866 

02 

S.S1S13E-02 

- 1.41 136E-01 

,392036 

00 

336 

4*631 936*01 

-9.O33O0E-O1 

3,395306 

ot 

2 .436706-04 

5.76386F 

62 

6 . 6666 IE— 02 

— 1 * 34 99 AE— 01 

. 392036 

00 

A E DC-T  R-78-6B 


937 

5*  S0BS3E-01 

-8.S2890E-01 

330 

6*330 306-61 

-7,93012E-O1 

339 

7.089312-01 

- 7.26036E— 01 

340 

T.77T76C-01 

-e.5204oe-ot 

341 

8.38893E-01 

-9.719SSE-9I 

342 

0.916T4E-01 

-4.BS50SE-01 

343 

9. 35590E-0 1 

- 3.94392E-01 

344 

9. 7021 OE— 01 

-2.99276E-0I 

345 

9*9Siaae-oi 

-2.011  86E-  01 

346 

1 • 0 10  27E  00 

-1*0109  7E-0I 

3«39330E  01 
3*39S30E  01 
3*3993DE  01 
3.39S30E  01 
3.39530C  01 
3*3953 0E  01 
3.39530E  01 
3.39B30E  01 
3*  39530E  01 
3*395392  01 


2»4)b70e-i>4 
2.43e?OE-0* 
2 *43570E— 04 
2.43ST0E-0* 
2.43870E-04 
£»43070E-04 
2t4?»70E-0l 

2*  43879E-04 
*««30T0E-Q4 
2*  433T0E-04. 


9.783AAE  02 
5.TS35QE  02 
5.7B388E  02 
9*T8368E  02 

9.73300c  02 

S.7fl3a8E  02 
5.70388E  02 
5*7Q9aO£  02 
5.T9304E  02 
9 • 70 3S BE  02 


8.195222-02 
9 .421310-02 
I.09S4BE-01 
1.1  504-2C— 01 
1.24990E-01 
U32900E-01 
1.394B7E-01 
1 •44604E— 01 
1 .4O436E-0 I 
1.50703E-0& 


-I.27919E-01 

— 1*1 87  7DE- 01 
-l«0flfl32E-tl 
-9.77970E-02 
-0.S7C98E-02 
-7. 200 TOE-02 
—3.92 tflOE-OZ 
-4.495932-02 
-A»02347E-0Z 

— 1 *51965E— 02 


I.39203E  00 
1.39203E  00 
1*392032  09 
1.J9Z03E  00 
1 « 392036  00 
1 *392036  00 
1 • 39203E  00 
1*392036  00 
U19203E  00 
U39203C  00 


00 

Ul 
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MACH  A NOZZLE  (G-1.24* 

ZE  = Z.40000E  01  STEP  NO*  78 


NO* 

X 

V 

P 

MHO 

0 

THETA 

PSI 

M 

1 

7 

151212-03 

2* 1 1067E-01 

3. 560866-01 

6.16520E-06 

1 .08538E 

03 

2*340088-03 

6*374196-04 

*061408 

00 

2 

4 

49 39 OE— 01 

- 1*091 72E- 03 

3.5321 *E  — 01 

8* 13894E— 06 

1.0BS66E 

03 

-4. 223018—04 

*« .599118-05 

•06453E 

00 

3 

7 

2$  J59E-01 

3 ■ 731  1 7E-01 

3,  52  E3  BE— 0 1' 

6« 1 3 3 67 E— 06 

1.0A572F 

03 

— 2. 627 34 E— 04 

-2.817208—73 

•06SI6E 

00 

4 

-2 

193188-41 

3*  791 33E— 01 

3, 52495E-0I 

6,1 3586E-06 

1 .035696 

03 

?*  I0792E— 03 

-2.586168*03 

• 064908 

00 

5 

-4 

350!  9001 

£.367598-03 

3.522676*01 

6.12367E-06 

1,005808 

03 

4.215968-03 

4.I2436E-04 

.0661 25 

00 

A 

-2 

110528-01 

-3.70560E-01 

3.52267E-01 

6.12567E-QA 

I-9B580E 

03 

3.23I43E-93 

4.478938-03 

.066126 

00 

7 

2 

26I4BE-01 

-3.  75347E-0I 

3. S2474E— 01 

6. L285TE-06 

1*08577E 

03 

2*678198-04 

2*953208-03 

.065 7 7E 

00 

a 

9 

91646E-01 

-1.9S31QE-03  . 

3*300796-01 

5* 82392b -66 

I.08909E 

03 

-1*215338-02 

-3*608688*04 

.103258 

00 

9 

7 

7*9708-01 

4 * 53  797E  — 01 

3«  2 734  BE  —01 

5.77375E-06 

1.089656 

03 

-1.199388-02 

— 7 .172108—03 

• 10462E 

06 

10 

4 

47107E— 01 

7* 694 a BE— 01 

3.281058-01 

S.7B4Q2C-06 

1.0A9S3E 

03 

-7.234458*03 

— 1 *109138-02 

•108258 

00 

11 

-9 

97679F-04 

8*044456-01 

3*2791  IE— 01 

5. 761 76E-06 

1.089568 

03 

-3.08736E-04 

—1 *505958-02 

1 0859E 

00 

12 

-4 

424456-0! 

7*61 3318-01 

3*  29905E-01 

5.8I010E-06 

t *08925E 

03 

6*403238-03 

*1 ,292276-02 

.105008 

00 

13 

-7 

82609E-0  t 

4*31 777E-01 

3.275B68-0I 

5. 77714E-06 

1.0B961E 

03 

1 *257218-02 

-8*540708-03 

,109198 

00 

14 

- a 

as 39  IE— 01 

— 1 * 3600 IE- 02 

3*2846 TE-01 

5.76966E-06 

1*089488 

03 

1.43BS2E— 02 

-1.10577E-O3 

1 0759E 

00 

IS 

-7 

731870-01 

-4*571018-01 

3*  2 73 TIE— 01 

S.  7 74 00 E— 06 

1 .009658 

03 

1*353678-02 

7.21503E-03 

•10958E 

00 

16 

-4 

39390E-01 

-7.694  316-01 

3.27707E-01 

5. 778868-06 

1 *089608 

03 

7*641658-03 

1.330978-02 

• 1 089  7E 

00 

IT 

-a 

67071 e— 03 

- 8*931 79E- 01 

3.274BJE-01 

5*  77567E— 06 
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-6*257396-04 

6.391QZE— 04 

4*150256 

00 

.093186 

03 

-9*386816-04 

5*673486-04 

4*150296 

00 

*093176 

03 

-3*456006-04 

3 *906576-04 

4*150216 

00 

.093206 

03 

3*436386-04 
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4*1 5143E 

00 

•093256 

03 

8*565096-04 

4* 09 23 OF-  04 

4*181066 

00 

*093306 

03 

1,332926-03 

5*858126-04 

4.151706 

00 

.09  3306 

03 

8.958456-04 

6*1 1623E— 04 

4*151 706 

00 

•09329E 

03 

-1*176786-04 

3.28346E-0* 

4*191 5B6 

00 

•093326 

03 

-6*131106-04 

7*71)876-05 

4*151946 

00 
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2.406992 

04 

-4.104792 

04 

1.040922  -Ot 
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03 

310 
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00 
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03 

139 
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03 
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00 
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03 
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00 
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00 

3.039092-01 

6.437072-06 

1 .003462 

03 
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00 

-2.367062 

•0 
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9.436I4E-06  • 

1 .003462 

•3 
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4 .902702 

oo 

-1.O30TQE 

00 
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1.60360E 

03 
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00 

-1.4692TE 

00 
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03 
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4.074362 

00 

•9.639092-01 

3.033162-01 

B.  42042E-O6 
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03 

340 

4.047022 

00 
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3.03 2612-01 
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03 

-lrt3MK-«* 
-iiMasae-Ai 
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-6.699002-04 
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9.8I9HC-04 
T.I474IE-Q4 
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•b4|«7IE-Q5 
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4.IU69E  00 
4.|929Se  00 
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4.I9C14C  00 
4*184962  00 
4.194992  00 
4. 1591 OC  00 
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NOMENCLATURE 


a Speed  of  sound 

h Static  enthalpy 

L Distance  along  bicharacteristic 

L,  M,  N Bicharacteristic  coordinate  system 

M Mach  number 

n^ , n^  Unit  normal  to  body  surface  in  z,  x, 
p Pressure 

q Velocity 

q i q I 


r*  Throat  radius  of  axisymmetric  nozzle 

R Gas  constant 

\ 

T Temperature 

u x component  of  velocity 

v y component  of  velocity 

w z component  of  velocity 

x,  y,  z Cartesian  coordinate  system 


directions 
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0 Mach  angle 

Y Ratio  of  specific  heats 

<5  Parametric  angle  for  bicharacteristic 

0,  \p  Flow  angles  defined  in  Fig.  1 

p Density 

SUBSCRIPTS 

1 i = 1,  2,  ...  refer  to  quantities  at  given  points 

ts  Stagnation  conditions  along  streamline 
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